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The total polycrystalline elastic and inelastic scattering cross sections are computed by 
means of the Born approximation with use of the Fermi (4 function) interaction between 
slow neutrons and bound nuclei. Isotopic disorder and magnetic interaction are aeglected. 
Numerical calculation for iron scattering ‘300-degree’’ neutrons yields an inelastic cross 
section which rises from 0.006 at T7=0°K to 0.192 for a scatterer temperature of 1000°K, in 
units of the “free” nuclear elastic scattering cross section. The total (elastic plus inelastic) 
cross section remains constant for scatterer temperature up to 250°K, but falls off from 0.97 
at T=0°K to 0.90 at T=1000°K. High and low temperature approximate results are also 
computed. Comparison of the results for inelastic scattering with the theory of the diffuse 
scattering of x-rays as developed by Zachariasen is successful in the limiting cases in which 


the two theories overlap. 





I. INTRODUCTION 


F a neutron scattered by polycrystalline 

matter possesses the same amount of kinetic 
energy that it had upon incidence, then it is said 
to have been scattered elastically. If the scatter- 
ing results in a decrease or an increase in the 
neutron’s kinetic energy, then the process is said 
to have been inelastic. In any given experiment, 
elastic and inelastic processes occur simultane- 
ously, so that by mere measurement upon the 
intensity of the transmitted beam one is not 
able to separate the contributions of the two 
types of scattering. Comparison with theory 
thus requires theoretical knowledge of both the 
elastic and inelastic scattering cross sections. 
The present work includes, in addition to a 
presentation of the theory of elastic scattering, 
a theoretical treatment of inelastic processes in 
which energy is given to or absorbed from the 
elastic vibrations of the scatterer by the neutron. 
Loss of energy by the neutron is accompanied by 


the excitation of one or more elastic vibrations— 
the emission of one or more sound quanta, or 
‘“‘phonons;” gain of energy by the neutron is 
accompanied by the quenching of one or more 
elastic vibrations—the absorption of one or more 
phonons. 

Although Halpern, Hamermesh, and Johnson! 
have treated the elastic scattering of slow 
neutrons by polycrystals, the subject is retained 
here since the present derivation is more com- 
plete and somewhat more general. Also, elastic 
scattering is here treated as a special case of 
inelastic scattering—zero-order process.” 

In all the expressions derived herein, there is a 
constant factor left undetermined, which in 
general can be determined only by experiment; 

10. Halpern, M. Hamermesh, and M. H. Johnson, 
Phys. Rev. 59, 981 (1941), henceforth referred to as HHJ. 

? HH] gives no explicit derivation of the crystal temper- 
ature dependence of elastic scattering; instead it borrows 
the exponential factor from the earlier x-ray theory as 


developed by Debye, Waller, et al. This procedure is justi- 
fied explicitly in the present work. . 
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this is the total nuclear elastic scattering cross 
section for the free nucleus. That this quantity 
should differ from the average cross section per 
nucleus for nuclei bound in a crystal lattice 
arises from the fact that nuclei so bound do not 
scatter independently; but interference influ- 
ences completely analogous to the familiar x-ray 
interference effects occur in the scattering of 
neutrons as well. In the present state of the 
understanding of nuclear forces, theoretical 
calculation of the free nuclear cross section is in 
general not at all possible; only the knowledge 
that nuclear forces are of short range (~ 10-" cm) 
is certain in the realm of features relevant to 
the present problem. This fact is expressed in 
the particular analytical form of the neutron- 
nucleus interaction employed. What is here 
sought is a comparison of the average elastic and 
inelastic nuclear cross sections, and the variation 
with scatterer temperature of these quantities 
as well as of their sum. 

For the calculation of all scattering cross 
sections, the present work employs the Born 
approximation with use of the Fermi (6 function) 
interaction between slow neutrons and bound 
nuclei. Attention is apparently limited to spin- 
free nuclei, but a slightly more detailed analysis 
reveals that the results are spin independent; 
for all scattering cross sections (including the 
free nucleus cross section) are multiplied by the 
same spin dependent factor in the case of 
scattering by nuclei of spin different from zero. 

Initially, consideration is restricted to a single 
crystal as scattering body. Since, however, we 
are concerned with the scattering from poly- 
crystalline substances, it is necessary to average 
the single crystal results over all possible crystal 
orientations. The single crystal is assumed suffi- 
ciently large to be approximated, for our pur- 
poses, by an infinite crystal. On the other hand, 
it is assumed that the crystal is sufficiently 
small to neglect the influence of the so-called 
“secondary extinction’’ in the polycrystal. This 
refers to the fact that the incident beam, in 
traversing the scattering specimen, is being 
depleted because of scattering by individual 
crystals, so that crystals farther back in the 
polycrystal are shielded from the full intensity 
of the beam, In HHJ, it is shown that for 
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reasonably small crystals this secondary extinc- 
tion is negligible. 

Throughout, it is assumed that the neutron 
beam is monochromatic: Every incident neutron 
possesses a definite energy, the same for all in 
the beam. By the use of velocity selectors, such 
beams are approximately attainable. 

From the outset, it will be assumed that every 
nucleus of the crystal scatters neutrons in 
precisely the same way. Thus we shall be limited 
to scattering by elements, and only those ele- 
ments in which the abundance of the possible 
rarer isotopes is negligible. In addition, the 
interaction between the magnetic moment of the 
neutron and the inhomogeneous atomic magnetic 
fields in the case of magnetic substances is 
neglected. Actually, its contribution to the 
scattering is always quite small compared with 
the purely nuclear part, if not completely 
negligible. 

The most significant result of the present work 
from the point of view of comparison with 
experiment is embodied in the variation of the 
total (elastic plus inelastic) scattering cross 
section with the temperature of the scattering 
polycrystal. So far as the author is aware, there 
have been no experiments performed which could 
provide the necessary check, so that the test of 
the theory must await more favorable times. 

On the theoretical side, however, there is an 
interesting check provided by a close similarity 
between the inelastic scattering of slow neutrons 
and the diffuse scattering of monochromatic 
x-rays. Except that the scattering of x-rays from 
the single atom is by no means isotropic, the 
only significant point of departure of the two 
problems is that for x-rays the energy of the 
scattered photon may always be set equal to the 
incident energy. That is, the photon energy is 
always large compared with its energy change in 
an inelastic process. For slow neutrons, the 
incident energy and probable energy changes are 
in general of the same order of magnitude. If, 
however, the energy change were negligible for 
neutrons as well, then the results should be 
precisely the same as one would obtain for 
x-rays of equal wave-length assumed to undergo 
isotropic scattering by the individual atom. For 
constant neutron wave-length—hence constant 
momentum—the energy is inversely proportional 
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to the mass, so that in the limit of neutron mass 
equal to zero the equivalence with the case of 
the supposedly isotropic x-rays is complete. 
(Actually the mass needs to be only so small as 
to render the neutron energy large compared 
with energy changes probable.) Modification and 
further development of the work of Zachariasen* 
on diffuse x-ray scattering provide a result 
suitable for the required comparison, as we 
shall see in full detail in ITI. 

Owing to the great complexity of the resulting 
mathematical expression for the inelastic cross 
section, it is difficult to obtain much useful 
information from mere inspection of the formula. 
In the special cases of high and low temperatures 
—compared with both the Debye temperature 
of the scatterer and the neutron ‘‘temperature”’ 
—simplifications occur, however, to render obvi- 
ous a number of interesting qualitative results. 
The work is then completed with a numerical 
computation of the nuclear elastic and inelastic 
cross sections for polycrystalline iron as a 
function of the temperature, for 300-degree 
neutrons. The strict calculations are supple- 
mented by the low and high temperature ap- 
proximate calculations. The results are discussed 
in the final section. 


Il. ANALYTICAL TREATMENT 


The fundamental question is, then: Given a 
crystal in a definite state described by the wave 
function ®, and a neutron incident with a 
definite momentum p, what is the probability 
that, as a result of the interaction between 
crystal and neutron, the neutron will later be 
found to possess a momentum p’, and the crystal 
will be found in the state ®,-? To this end, we 
employ the Born approximation, justified by 
Fermi‘ for the case of slow neutrons scattered 
elastically by bound nuclei. The justification of 
its use for inelastic scattering lies in the circum- 
stance, as will be shown, that the inelastic 
amplitude is in general small compared with the 
elastic. 

It is necessary to know the solutions of the 
Schrédinger equation corresponding to all states 
of the neutron and crystal in which the inter- 
action between them is neglected—that is, the 


*W. H. Zachariasen, Phys. Rev. 57, 597 (1940). 
‘E. Fermi, Ricerca Scient. 7, Part 2, 13 (1936). 
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product of neutron wave functions by crystal « 
wave functions. The neutron functions are 
merely exp (tk-r) for all k, where k= (p/h) is the 
propagation vector of the neutron; r is the 
neutron coordinate vector. 

The treatment of the crystalline potential 
energy by the introduction of normal coordinates 
is carried out in the standard fashion of ex- 
panding the individual nuclear displacements 
from equilibrium in terms of a series of plane 
polarized standing elastic waves, each of which 
is characterized by a propagation vector q and 
one of three mutually orthogonal polarization 
directions (j= 1, 2, 3) corresponding to the given 
q. Most of the details are omitted, since they are 
abundantly covered in the literature.® 

If the position of each nucleus in the crystal 
at any instant is given by the vectors R,, we 
may write: 


R,=e+u,, (1) 


where @ is the equilibrium position of the eth 
nucleus and u, the instantaneous deviation from 
equilibrium. For not too high temperatures, the 
potential energy of the crystal is expressible as 
a quadratic form involving the components of 
all the displacements u,. To obtain the most 
simple expression for this energy, we introduce 
the variables y, and u,, where 


u, = (2/mN)! >, a.(y, CoS @-G.—, Sin @-q,). (2) 


The index s stands for the double index (q, 7), 
and the range of values assumed by s is such 
that for each q the three mutually orthogonal 
directions of a, (unit polarization vector) are 
summed over, and then one-half (determined by 
any plane intersecting the origin) of the q space 
to be defined presently. N is the number of 
scattering centers in the crystal, m the mass at 
each scattering center. 

In terms of the y, and 4u,, the crystalline 
potential energy may be written 


U =i. (1/2)a.?(y.?+u,"). (3) 


Because of the fact that U may be so written, 
the y, and yu,, the “normal coordinates,’’ may be 
regarded as the coordinates of harmonic oscil- 


lators of unit mass and circular frequencies w,. 


5M. Born and K, Sarginson, Proc. Roy. Soc, A179, 69 
(1941), 
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To each s, then, there are both a ‘‘y-oscillator”’ 
and a ‘“‘y-oscillator’’—therefore one of each to 
every (q, 7) combination. 

The definition of the q space over which 
summation is taken is fixed by the usual require- 
ment of periodicity of u, in a “fundamental 
parallelopiped”’ of the crystal. Since all results 
must be independent of the size of this parallelo- 


piped, we take it to be that containing the NV 


lattice points of our crystal, for sake of definite- 
ness. We write: 


@ = 018; + 0282+0383, (4) 


where @1, 2, a3 are the basic vectors of the lattice 
(of dimensions length), and o;, a2, o3 are integers 
peculiar to the eth lattice point. If bi, be, bs are 
the basic vectors of the reciprocal lattice— 
b,= (a Xas)/[ar ° (a2 Xas) |, etc.—we have: 


q=q1b1+9¢2b2+Qsbs. (5) 


Then 
0° Q = 0191+ 0292+ 0393. (6) 


Thus g; must have only such values that if o; 
runs through its “fundamental interval,’ then 
o.qi changes by some integral multiple of 27; 
for only then will u,, given by (2), remain 
unaltered. If the fundamental parallelopiped is 
such that it is traversed exactly once when all 
three of the o; are augmented by the integers g;, 
then we have 


gi=2Qi/z:, (7) 


where Q; is a positive or negative integer. 

Also, q space is limited by the fact that each 
oscillator actually represents a plane polarized 
standing wave of half-interval of periodicity in 
any of the three basic directions a; equal to 
(x/q:). Since it is physically meaningless for the 
half wave-length of a sound wave to be less than 
the distance between successive lattice points in 
the crystal, we must have (x/|q;i|)21, or 
—x£qi<r. By (7), then, the integers Q; are 
limited by the inequalities — (g;/2) < Qi (g;/2), 
whence gq space is limited to gigogs=N points. 
(Unity is neglected compared with g;.) We 
employ this fact later to determine the upper 
limit of the absolute value of q. 

Now, the potential energy U, given by (3), is 
that of a system of 3N uncoupled harmonic 
oscillators of circular frequencies w, and each of 
mass unity. (There are N/2 points of half q 


space, to each q the three values j=1, 2, 3; to 
each (q, 7) combination are both a y and a y— 
that is, 3N terms in (3) altogether.) The wave 
functions are well known to be® Hermite func- 
tions of argument y(w/h)' and u(w/h)* for the 
respective cases. The normalized solution for a 
y-oscillator in the mth energy state is 


n(¥(w/h)*) =[(w/he)?/(2"n!)4] 
Xexp (—wy?/2h)Hn(y(w/h)*), (8) 


where /7, is the mth Hermite polynomial defined 
by the generating function 


exp (—2?+ 22x) = >> (2"/n!)H,,(x). 


The energy of the oscillator in the mth state is 
(n+1/2)hw. 

Thus we may write the total wave function 
describing the thermal agitation of the crystal 
as the product 

®,, = II On, Pn,» 


$182 


where On, = pn, (ys(ws/h) }), rn, = Gn, (He(ws/h) 4), 


Each of the indices s;, sz clearly takes on the 
3N/2 values assumed by s. 
Finally, the total wave function of the system 
neutron plus crystal without interaction is 
WV =exp (ik-r)®, =exp (ik-r) J] dn, On, 


$132 
Because of the interaction between crystal and 
neutron, the state of either or both may be 
changed, so that we denote the wave function 
after interaction by 
Te /. P ‘ 
W’ =exp (ik’-r) IT ng, n;,. 

In the Born approximation, the probability of 
transition from a state V to a state WV’ because 
of the interaction energy V between neutron and 
crystal is proportional to | M|?, where 


u=f f fvsvvdrardy; 


the integrations are carried out over the com- 
plete ranges of the variables; dr=volume ele- 
ment in neutron coordinate space, 


dy=II dys, du=]I dps. 





*Any elementary text on quantum mechanics—eg., 
V. Rojansky, Introductory Quantum Mechanics (Prentice- 
Hall, Inc., New York; 1938), first edition. 
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Fermi has shown‘ that for bound nuclei 
scattering neutrons of wave-length long com- 
pared with nuclear dimensions, the interaction 
energy between the neutron and individual 
nucleus, position R,, is given by 


V,=A,6(r—R,), 


where, for spin-free nuclei (cf. introductory 
remarks that the results also apply to a spin 
dependent interaction) A, is a scalar depending 
upon the type of scattering nucleus. We are 
restricted to scattering by identical nuclei by 
setting A,=A, a constant for all «. 6(r—R,) is a 
three-dimensional Dirac 6 function. For “‘room 
temperature” neutrons, the wave-lengths are of 
the order of atomic dimensions, so that we may 
employ the above form of V,. The total inter- 
action between neutron and crystal is given by: 


=>, A,i(r—R,) = AD, 3(r—R.), 


where the summation is taken over all the nuclei 


of the crystal. 
We write, therefore, since ®,, is a real function, 


u=Af f fex (ir-(k—k’)] 
XX. 8(r—R.) O,O,-dtdydu 


=A Ey ff exp (R-e-K)] 
XII on, on; dye IT bn, dn; dure 








Tar, 2(t) = oe +h) in} exp (tx/a) exp (—x?/a?) Hnia(x/a)Ha(x/a)d(x/a) 
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=A >. exp [io-(k—k’) ] 


x f fexp Liu, - (k—k’) ] II On, bn, dys 
<I On, bn; Cus, 


with the use of (1). The expansion of u, according 
to (2) gives, finally, 
M=A >. exp [ie-(k—k’) ] J] In, i Tnzm, » (9) 


$182 


where 
} es few (iF 7s) Onin de, 


Fu= (2/mN)*(k—k’) -@s; cos o-qQs: 
and 


Fi. = —(2/mN)4(k—k’)-as: sin @-qs:. (10) 


It is necessary, then, to evaluate integrals of 
the form 


Insnald)= J exp (b0/a)dnsn(2/a)a(x/a)dex 
a (t=iaF,’, a=(h/w,)*). (11) 


For n’ has been written m+); i.e., the result 
corresponds to a change from state m to state 
n+n of the given oscillator. 4 may be positive 
or negative, but necessarily an integer (as is m). 
The evaluation follows: 

According to (8), 


= [22+ !(n +2) in} [exp (tx) exp (—x*)Hns(x)Ha(x)dx 


=[22"+°n!(n+ ) in} [exp [ —(x—t/2)*] exp (9/4) Hasa(x)Ha(x)dx 


= [22"+4n!(m-+ A) le }-! exp e/a) f exp (—2?)Hasa(2+1t/2)H,(2+t/2)dz 


exp (1?/4) 





m >> n+xX , 
~ [22+ I-A) Mahe i a ae N , +f exp (—2*)Hn4r-6(2) 1 n_,(2)dz 


exp (17/4) 





n—» _ _ ; 





Po: on . yao nad 
~ Enl(n+a)!Tt ST (ota ( 
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The second-from-last step is obtained through 
the use of the expansion 


Halety)=¥ (") (2y)'Ha-ol2), 
as one may show from the familiar recursion 
formula (d/dx)H,(x) =2nH,_,(x) and use of the 
Taylor expansion of H,(x+~y) about the point x 
in powers of y. Otherwise, the expansion of 
H,(x+y) may be obtained directly through use 
of the generating function for Hermite poly- 
nomials. The next-to-last step in the sequence 
above is given by the orthogonality relation 


f exp (—2?)H,(z)H,(z)dz = 6,82’! x}. 


The summation limits for » (and 6 above) are 
provided naturally by the presence of denomi- 
nator factorials; the result is unchanged if one 
sums from zero to infinity. 

The quantity |M|? (proportional to the 
scattering cross section), for scattering processes 
in which the quantum number of the sth oscil- 
lator changes by the (positive or negative) 
integer \,, is given according to (9) by 


| M|?=|A|? > exp [i(o—o’) -(k—k’) ] 


xII Jn, +A, ,m, Ff +h, Ms (13) 
$182 1 1 1 2 2 2 
where 


oo’ o a’* 
Inia, =I nin Tniny 


If for all s, \,=0, the scattering is elastic; this 
is the case which has been considered by HHJ. 
Positive 4, means loss of energy, of an amount 
\,fw,, by the neutron to the sth oscillator; 
negative \, means the absorption of this amount 
of energy by the neutron from the sth oscillator. 

According to (11) and (10), ¢ is proportional 
to N-!; since N>1, then (for w,#0) «1. Thus 
it is necessary to keep only the lowest powers of 
t in the expression (12) for Ins, .(¢). For A=0, 
it is necessary to expand to the quadratic term 
in ¢; to this approximation, 


In, n(t) =1+(14+2n)(t?/4). (14) 
The term in #, proportional to N~, must be 


kept since the quantity M involves the product 
of approximately 3N such factors; higher order 
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terms involve higher powers of N-' and are 
therefore negligible for N>1. 
For \= +1, only the linear term must be kept: 


Int, n(t)=[(m+1)/2 J; Ins, n(t) =(m/2)4#. (15) 
In the transition probability, proportional to 
| M|?, the lowest appearing power of ¢ due to 
\=+1 for just one oscillator (and \=0 for all 
the others) is the second, which is proportional 
to N-'. Since, however, there are 3N oscillators 
which may so participate, the probability that 
any oscillator so participate is the sum of the 
single probability over all the oscillators. The 
probability is thus rendered finite, no matter how 
large N; the terms in higher powers of ¢ still 
contribute negligibly for large N. If more than 
one oscillator, say B(<JN) in all, undergo changes 
in nm by \=+1, then the factor N-* appears in 


3N 
the quantity | M|?. Summation over the ( P ) 


combinations of 3N oscillators taken 6 at a time 
renders the total probability finite, but the con- 
tribution of higher powers of ¢ remains negligible. 

It is next shown that processes for which 
|| 22 are completely negligible so long as N>1. 
By (12), to the lowest power of ¢: 


Tn+a,n(t) =[n!(n +A)! TAL (a +A) V/A 2-2; 


lana, n(t) = [nm '(n —a) '}-3(n '/a !)¢e2-@/2 
(a= —A>0). 


For || 22, the squares of these quantities are 
proportional to N-*, or higher powers of N~. 
Thus even summation over the 3N oscillators 
still renders the probability negligible for N>1. 

From the expansions (14), (15), it is seen that 
for N>1 the ratio |Zn+1,2/?/|Znn|? is propor- 
tional to |¢|*, or to (h/w,)(F.’)?, by (11). Since 
| In+1,n]2 must be summed over all 3N oscil- 
lators, the ratio of the probability of a first-order 
process to that of a zero-order process is roughly 
proportional to 


L=N(h/w,)(F,*)*=(Ak?/2mw,) 


el 


to order of magnitude, according to (10). Here 
mo is the neutron mass, so that (h?k?/2mp) is the 
incident neutron energy; hw, is of the order of 
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the energy of the sth oscillator. For room 
temperature neutrons the ratio of these energies 
js, on the average, of the order of unity for 
many substances. Hence for such cases, what 
counts in the ratio L is the factor (mo/m), the 
ratio of the neutron mass to the mass of the 
scattering nucleus. Thus, in general, for elements 
of large nuclear mass L is small compared with 
unity; i.e., inelastic processes are much less 
probable than elastic ones. Also, since this same 
ratio L is involved in the comparison of the 
probability of processes involving two oscillators 
with the probability of those involving only one, 
the former are much less likely than the latter 
for (mo/m)<1. It is for this reason that the 
present work shall confine itself to single oscil- 
lator first-order processes only, in addition to 
zero-order processes. For Fe as_ scatterer, 
(mo/m)=1/56, so that the limited consideration 
is well justified in this case. The role of the 
factor, ratio of incident neutron energy to hw,, 
in L is not so clearly defined as a function of 
neutron energy. This is because of crystal 
interference effects, which depend, for processes 
of all orders, quite strongly upon the incident 
energy and, therefore, may play a far greater 
part in determining the relative magnitudes of 
types of scattering than do these simple ratio 
considerations. 

To find the actual probability (eventually the 
cross section) of each of the oscillators performing 
the transition from nm, to m,+A, (A.=+4,/ or 0 
for all s, for our purposes), it is necessary to 
multiply | /|*, (13), by the probability that the 
oscillators originally be in the states n,. Since 
before interaction with the incoming neutron the 
crystal is in thermal equilibrium, the Boltzmann 
probability that the sth oscillator be in the state 
n, (energy =(",+})hw,) is given by (1—x,)x,", 
where 


x,=exp (—hw./koT), x,"*=exp (—n,hw,/koT), 


ko=Boltzmann’s constant, and 7’ =crystal tem- 
perature, Kelvin. Hence the probability that all 
the oscillators initially be in the states m, is 


II (1x, Ja, (1—2,,)ay, —the indices 5}, S2 re- 


ferring, respectively, to the y and 4 oscillators, 


NEUTRONS 7 
as previously. We form, then, the product 


Q=T1 (1— x), (ada, | M2, 


#142 
whence from (13), 


Q=|A|? ¥ exp [i(o—o’)-(k—k’)] I] (1—x,) 


eo’ 8182 


1 oo’ "es: aa’ 


ge (1—x, x. Jn, +A, My, * 
1 1 : 2 2 2 -, | 


Owing to the circumstance that the energy 
levels of the harmonic oscillator are equidistantly 
spaced, an important quantity is the quantity Q 
summed over all m, (zero to infinity) for all s. 
This gives (when properly normalized) the prob- 
ability of oscillator transitions by energies \,/iw,, 
regardless of the initial states m,. Applied to Q, 
this operation gives: 

P=|A|? > exp [i(o—o’) -(k—k’) ] 
XT (1—a)(1—ae) Yarn, Vegas,» (16) 


where 


"Ns oo’ 


Yn. = _ Xs J npthy n,° 


Dropping for the moment excess indices, we 


write 
Y, = on Jat, ny 


with Jay n=Lnprnlb)Inua(w), and t=iaF,’, 
w=—iaF,’, a=(h/w,)'. The calculations of 
Yo, Ya: follow: 

According to (14) 


Jn,n=1+(1+2n)(2+w?)/4, 
so that 


Yo= Dn x*+L[(P+w*)/4][Dn x" +2 Dn nx”) 
= (1—x)-'+[(#+w")/4] 
<((1—x)-!'+2x(1—x)-*], 
Dax*=(1—x)-', (x<1), 
Dn nx" =x(d/dx) >, x" =x(1—x)-*. 
Finally, 

Yo=(1—<x)—! exp [(#+w*)(1+x)/4(1—<x) ], 
since 1+z=exp (z) for z1. Or, if we restore 
the missing indices, 

Y¥s"o=(1—x,)-! exp [—(h/en) 
XK ((Fi")? + (Fi)? } (1 +%.)/4(1 — x) J, 


since 


and 


(17) 
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t and w having been replaced by their respective 
values given directly above. 

In what follows, it is convenient to obtain the 
quantities Yi; as multiples of Yo. To this end, 
it is necessary to expand J,+1,, to powers higher 
than the first of ¢ (or w), as in (15). From (12), 
we get 


Tn+1,n(t) =t exp (/4)[(n+1)/2 }L1+n(#/4) J, 


whence 


In+i,n= (tw/2) exp [(t?+w*)/4 ](n+1) 
x [1+n(2+w?)/4], 
Vi = >on X°S 41, n= (tw/2)(1—x)-? 
Xexp [(#+w?)(1+x)/4(1—x)] 
= (tw/2)(1—x)-'Vo, 


as one finds by applying >>, ?x" =x(d/dx)>>, nx” 
=x(1+x)(1—x)-* (the sums >>, x" and >>, x" 
having been evaluated to obtain Yo), and 
1+2=exp (z) for <1. Finally, 


Yo. =(AF PF." /2w(1—x.) ]¥i0. — (18) 


Similarly, it is found that 

Ys" =[hx. FF," /2w,(1—x.) ]¥eo=%eYer. 

Since x,=exp (—hw,/koT) is less than unity 
for finite T, the relation Y_,=x Y; expresses the 
physically obvious fact that processes in which 
the neutron loses energy are more likely to 
occur than those in which it gains the same 
amount; at J7=0, it is impossible for a neutron 
to absorb energy from a crystal lattice. Only 
for temperatures such that k)T>>/hw, are energy 
gain and loss equally probable. 

For elastic scattering, the quantity P given 
by (16) is denoted by Po; for first-order processes 
in which the fth oscillator undergoes a change 
n; by +1, P is denoted by P’:. According to 
(16), (17), and (18), 


P,=|A|* exp (—2W) 
X |. exp [te-(k—k’) ]|*, (19) 
Pf =|A|?exp (—2W)[h/2w;(1—x,) ] 


X | L. exp [to-(k—k’) ]Fy*|*, (20) 
and 


P!,=x,Py; (21) 
W= Lo [(h/cons) (Fes)*(1 +2081) /4(1 —xs)] 
+ Xa [(h/wea) (Fea) *(1 +22) /4(1 — 2x02) J. 


The removal of exp (—2W) to the left of >, is 
justified directly ; W is actually independent of ¢, 

Since there is a one-to-one correspondence 
between the terms of the two sums of W, this 
quantity may be simplified by dint of the fact 
that ws; =ws2=w, (and therefore xs; =xs2=x,) for 
corresponding terms. From (10), we obtain, for 
corresponding sj, Se, 


(Fi) + (Fix) =(2/mN)[(k—k’) -@, }2. 
Consequently, 
W=(h/2mN) ¥, [(k—k’) a, ]2(1 +24) /we(1 —x,) 


has clearly no dependence on go. 

In the usual fashion, the sum over s is replaced 
by an integral over the half of q space. For 
N>1, this is valid. Also, we employ the Debye 
approximation for the frequency spectrum— 
namely, w=cg, where c is the constant average 
sound velocity. (We adopt the average sound 
velocity—to be characterized in terms of the 
average Debye temperature—instead of keeping 
the individual transverse and longitudinal veloci- 
ties for the sake of great simplification of the 
work involved.) The assumption of a sound 
velocity independent of wave-length actually 
holds in a continuum, therefore for wave-lengths 
long compared with the lattice constant. That 
it breaks down for short wave-lengths (large w) 
is not a serious difficulty so long as we use only 
integral properties of the frequency spectrum. 

With an average sound velocity assumed, the 
values of w, (and hence x,) are the same for the 
three mutually orthogonal directions of the unit 
vectors a, for a given q. Thus we may write 


W = (h/2mN)(k—k’)? b.’ (14+2%,)/w.(1—x,), 


where the prime of >.,’ indicates that only the 
N/2 points of half q space need be summed over 
—the summation over j (indicating the three 
directions of polarization) having been carried 
out for all q. Now >-,’ actually means summation 
over the three sets of discrete variables q:, 2, qs, 
given by (7). To convert this sum into an 
integral, we must first multiply the summand by 
[.N/(2x)* ]dqidqedqs since (7) shows that there 
are N/(2x)* points per unit ‘‘volume”’ of gigs 
space. It is yet more convenient to convert >,’ 
into an integral over the variables gz, gz2, 4s, 
where these are the Cartesian components of q 
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(and of dimensions inverse length). This necessi- 
tates further multiplication by B, where 


dqidq2dq3 = Bdgz,dqz2dq3. (22) 


For regular lattices, the transformation from 
reciprocal lattice to Cartesian components is 
linear, the Jacobian B constant. It will be shown 
in Section VI that B is the volume per nucleus 
in the crystal. Finally, we employ the usual 
artifice of replacing the region of q space by a 
sphere containing the same number—N/2—of 
points. The volume element is thus a hemi- 
spherical shell, or dgz1dgz2dgz3 = 2%q°dq ; the radius 
of the sphere, to be determined directly, is qo. 
To insure that the number of points in the 
hemisphere is V/2, we must have 


go 
(1/2) =[N/(2n)*1B f “2ra*da: 
; " (23) 
go® = (62°/B). 
At length, with application of cg=w, cgo=wo 
(the maximum circular frequency), we find (the 
now unnecessary index s having been dropped) : 


W=(38/4ma)(e—K’)* { [e(1 +x)/(1—x) Jdw 


= (3h?/2mk)@)(k—k’)? 
x {(1/4)+(T/0)?Q(0/T)}, (24) 


where use has been made of the fact that 


x=exp (—hw/koT), and O=(hwo/ko) is the 
Debye temperature of the scattering substance. 
Here, 


Qx(2) = f "{B/Lexp (8)—1]}d8. (25) 


A. Elastic Scattering 


In order to convert Po, (19), to an actual 
cross section (per unit solid angle and per scatter- 
ing nucleus), we must multiply by (m?/4rh‘N), 
where mp is the neutron mass. Thus 


Go(k—k’) = (X/44N) exp (—2W) 
x | 2. exp [to-(k—k’) ]|*, (26) 


is the elastic scattering cross section per nucleus 
per unit solid angle for the single crystal; here 
X=(m,?/rh*)|A\? is the total elastic cross 
section of the free scattering nucleus. 

Our system of basic lattice vectors is neces- 
sarily such that for each integer triple c,o203 


(O$o;<g;) there is one and only one crystal 
lattice point. [The definition of g—>»>1—pre- 
cedes relation (7) directly ; we shall often neglect 
unity in comparison with g;. ] Thus, since 


3 3 
o=) o1a;, (kK—k’) => (ki—k,’)bi, 
i=] i=l 
and : 


o-(k—k’) => oi(ki—k,’), 
t= 
following (4), (5), (6) (where Rk; and &; are 
clearly the reciprocal lattice components of k 
and k’, respectively), 


|1/?=| 5. exp [ie-(k—k’) ]|? 


3 0; 


> exp [io(ki—ki’)]| 


1 o;=0 


o 3 sin? [gi(ki—k,’)/2] 
~ ica sin? [(k;—k/’)/2] 





3 
=(2r)*N J] 5(ki—k,’ —2x7;). (27) 
i=1 
Here 6(k;—k,;’—2x7,) is a Dirac 6 function de- 
fined by the integral property /f(x)8(x—b)dx 
= f(b) for b within the range of integration, =0 
for 6 without the range of integration; 1; is a 
positive or negative integer. The vector 
3 


<=). rb; 
i=l 

therefore connects two points of the reciprocal 
lattice, and is loosely referred to as a “reciprocal 
lattice vector.’’ The validity of replacing the 
sine-squared ratios by 6 functions rests upon the 
facts that g;>1 and that functions by which 
|>01|? is multiplied undergo no violent variation 
in the region in which | >, |? differs appreciably 
from zero; also, only the integral properties of 
|o1|? are required in the present work. 

If instead of using 6 functions whose argu- 
ments are reciprocal lattice components of 
vectors, we introduce 6 functions in which the 
arguments are the Cartesian components of 
these same vectors, we find from the definition 
of 6(x—6) above that 


3 
[I 6(&:—k,’ —2e7,) = B5(k2;— ke, — 22) 


i=l 
X 5(k22— kre — 2wrze)8(kz3—ke3 — 223) 
= Bd(k—k’ —2r-), 
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where B- is the reciprocal of the Jacobian, 
defined by (22), for the transformation from 
reciprocal lattice to Cartesian components. The 
three-dimensional, or vector, 6 function is 
merely an abbreviated way of writing the triple 
product for which it stands. It has the additional 
advantage of independence upon the orientation 
of the Cartesian coordinate system. Thus (26) 
finally becomes 


Go(k —k’), =2°B-1X exp (—2W)6(k—k’ — 292). 


This result expresses for neutrons the familiar 
x-ray Laue conditions for elastic scattering by a 
single crystal—namely, that only for well-defined 
(the better the larger NV) angles of incidence is 
there elastic scattering at all for an incident 
monochromatic beam, and then only at certain 
angles defined with the same degree of sharpness 
as is required of the angles of incidence. 

We now pass to the case of the polycrystal by 
averaging over all crystal orientations. We do 
this by averaging over all directions of +, which 
for any given value has a fixed orientation 
relative to the crystal axes. We choose a polar 
axis along the direction of (k—k’) and denote 
the polar coordinates of + by 1, 6,, ¢,, so that 


E,(k—k’), = (1/47) f _ f Go(k—k’), sin 6,40, 
=(x/2)B-X exp (—2W)S; (28) 


is the polycrystalline elastic cross section corre- 
sponding to the reciprocal lattice vector +. 
Here, clearly, 


2 td 
sf def 6(k—k’ —2r2) sin 0,40, 
0 0 


= (1/4rert) [ f a(r21)8(r) 
Xx 6( | k—k’| — 2123) | cos 6,|—'drz,drz2, 


since the x; and x, components of (k—k’) 
vanish owing to the choice of polar axis, and the 
Jacobian 


D(—2xr2;, —2x722/0,, gr) =42°r* sin 6, cos 6,. 


For a non-vanishing integrand, tz3=7, so that 
cos 6,=1. Hence the polycrystal cross section 
per unit solid angle per scattering nucleus for 
given « is, on substitution of the result for S; 
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into (28), 
Eo(k—k’), = (X/8xBr*) exp (—2W) ; 
x 6(|k—k’| — 2x7) 
= (X/82Br*) exp (—2W,’) 
X 6(2k sin (0/2)—2mr), (29) 


since for k=k’ (elastic scattering), |k—k’| 
=2k sin (6/2); 6 is the angle of scattering. This 
result expresses for neutrons scattered elastically 
the familiar polycrystalline (or powder) scatter- 
ing of x-rays along certain right circular cones 
about the line of incidence—so called ‘‘Debye- 
Scherer rings.’’ The particular cone depends, of 
course, on the magnitude only of «, but we must 
keep in mind that in general there are several 
examples of + in the reciprocal lattice to a single 
absolute value r. Since sin (0/2) <1, only those 
vectors « for which r<(k/r) provide Debye- 
Scherer rings of scattering. 

To obtain the total polycrystal cross section 
per nucleus, we integrate over all 6, ¢—the polar 
angles of neutron scattering. For simplification, 
in the exponent W of (29) we have replaced 
(k—k’) by 2+, its equivalent for permissible 
angles of scattering, and have so replaced W by 
W.*’. Then the only 6, ¢ dependence is in the 
6-function factor. By (24), 


W,’ = (62°h?7?/mko@) 
x { (1/4) +(7/0)?Q1(9/T)}. (30) 

It is here for the first time that we use the 
independence of X on 6, g—i.e., the fact that 
purely nuclear scattering of slow neutrons is 
isotropic from the single isolated scatterer. 

Since Eo(k—k’), is independent of ¢, the total 
cross section for given ¢ is 


E.,.=2x | Eo(k—k’), sin 6d6 
, = (X/4Br?) exp (—2Wo')Ss, 
where 


S= f 5(2k sin (0/2) —2xr) sin 6d0 
0 


= (1/k?) f "(2k sin (6/2) —2er) 


= 2er/k?, X {2k sin (6/2) }d {2k sin (6/2)} 


whence, finally, 
E., = (4X /2Brk*) exp (—2W,’). (31) 
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For simple cubic structure (B=a*, volume of 
unit cell), this reduces to the result of HHJ for 
elastic scattering. 

The total elastic polycrystal cross section per 
nucleus is obtained by summing over all recipro- 
cal lattice vectors * for which 7 < (k/2) is satis- 
fied, together with other conditions on + de- 
pending upon the type of lattice: 


Ea= Dr Ex 
= (4X /2Bk*) >, (1/r) exp (—2Wo"). (32) 


B. Inelastic Scattering 


We pass now to inelastic scattering processes 
of first order; initially, consideration is confined 
to phonon emission—loss of energy by the 
neutron in exciting a sound quantum. From the 
definition (10) of Fy’, we write (20): 


Pf =|A|*{h/4mw,(1 —x,)N J[(k—k’) -a; P 
Xexp (—2W)| >. {exp [ie (k—k’+q,) ] 
+exp Lio: (k—-k’—q,)]}|*. (33) 


The + sign arises from the fact that F;” may 
contain either a cosine or sine factor (according 
to whether f refers to a y- or yu-oscillator). The 
ambiguity may be removed by recollection of the 
limitation to consideration of just one-half of q 
space (determined by any plane through the 
origin). For q;#0 and N>1, only one term of 
the summand in (33) will give a contribution 
appreciably different from zero for suitable k, k’ ; 
each sum is essentially a 6 function of argument 
different from that of the other. Hence, if we 
now open consideration to all of q space, thereby 
departing from the earlier restriction, it is 
permissible to omit one of the terms of the 
summand. That is, the one term provides for the 
“‘non-zeros” of the other by the mere introduc- 
tion of the previously eschewed values of qy of 
opposite sign. 

To convert P;/ to the cross section per unit 
solid angle per scattering nucleus, one multiplies 
by (m?k'/42°kh'N). Thus 


Gi (k—k’) =k’ AX /16rkmw (1 —x,)N?] 
X[(k—k’) -a, ? exp (—2W) 
X | 2. exp [to-(k—k’+q,)]|? (34) 


[where again X =(m,?/rh*)|A|? is the total 
elastic cross section for the free nucleus] is the 
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inelastic scattering cross section per unit solid 
angle per nucleus for single phonon emission, 
corresponding to the fth oscillator. 

From (21), one finds directly that 


G!(k—k’) =x,/Gi/(k—k’). (35) 
Comparison with the treatment given the 
quantity | >-,|?, (27), yields the result: 
G,/(k—k’), =[wk'AX /2kmw,(1—x,) BN] 
X((2re—qy) ay? 

Xexp (—2W,)i(k—k’+q,;—22re), (36) 

where ¢ is a reciprocal lattice vector of the type 

introduced in (27). Since G,;/(k—k’), vanishes 


unless k —k’ = 2912 —q,, this substitution has been 
made in (36); thus W has been replaced by 


W /* = (3h? /2mkoO) (242 —q,)? 
x { (1/4) +(T/0)*0:(9/T) }. 
Now the only dependence upon the angle of 
scattering rests in the 6 function. 
We next compute the total single crystal 
scattering cross section for given k’, but inte- 


grated over all orientations of k’. The quantity 
required is 


2r r 
Girm f ae | G,/(k—k’), sin 0d’ 
0 0 


=[9k’'hX/2kmw;(1—x,) BN] 
X([(24e—qy)-ay P exp (—2W;")Ss, (37) 


where k’, 6’, g’ are the polar coordinates of k’, and 
2r r 
Ss -{ ae 5(k —k’+q,—22) sin 6’dé’. (38) 
0 0 


We choose as polar axis the direction of 
(2%*—q,), so that the x; and xz components of 
this vector difference vanish. Also, we use the 
variables 


kn=k' sin @’ cos ¢’ and kz:=k’' sin 6 sin ¢ 
for the integration, the Jacobian of the trans- 
formation being 


Dk, k'z2/8’, gy’) =k” sin 0’ cos 6’. 
Hence 
S3=(1/k’)? f f 5(k2;—k2)5(k22—k22) 


X 5(k23—k23— | 2ne—q,|)| cos 6’ |—'dkrdkoo. 
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The first two 6 functions yield kz;=k2 and 
kz2=kz2 for a non-vanishing integrand. The con- 


sequent value of k:3 is obtained through the 
condition of energy conservation: 


(h2k!?/2mo) = (h2R?/2mo)—hwy. (39) 


Because of the equality of the first two compo- 
nents of k and k’, this becomes 


keg = + (kz3—2mow,/h)' = +Z, 
and 


|cos 6’| =(|kz3| /k’) =Z/k’. 
Hence, finally, 
S3 = 6(Rzg FZ — | 2re—q,|)/k’Z. (40) 


Next, G“" is averaged over all directions of k 
—a procedure physically equivalent to averaging 
over all crystal orientations, so that the result is 
the total polycrystal cross section for given k’: 


25 © 
Ey (1/4e) [ de | G,"" sin 6,.d6,, 
0 0 


where k, 6%, gx are the polar coordinates of k. 
From (37) and (40), we have 


Ey =[whX/4kmw (1—x,) BN] 
<[(24re—qy)-a, ? exp (—2Wy") Si, 


where, after integration over ¢x, 
Sif 5(R23 FZ — | 2re—q,| )Z-1 sin 6,4 0,.. 
0 


It can be shown that the argument of the 
6 function has the single root 


k cos 0 =(|2re—qy|?+2mow,s/h)/2|2e2—q,| ; 


but since this solution is obtained through 
squaring the term preceded by *, it may pertain 
to either sign. It turns out that for one of the 
signs this solution always holds, for the other it 
is extraneous—which condition depends on 
whether |2r<—q,!? is less than or greater than 
2mw;/h. (In the case of phonon absorption— 
i.e., where w; may be replaced by —wy in the 
conservation of energy and resulting formulae— 
it is for the plus sign alone that the above 
solution—with wy; replaced by —w,;—holds.) 
That this solution gives no new condition is 
ascertained by the fact that it is derivable from 
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the already imposed relation, (k —k’) = (2r<—q,), 
if one uses the fact that 


k-(24e—q,) =k|2x2—qy| cos 6, 


since the polar axis has been taken along the 
direction of (24r<—4q,). 

The ambiguity can be obviated by the change 
of variable: 

2=k cos 0. (k? cos? 6, —2mw,/h)!, 
dz ; Rk? cos 6% sin 0% 
—=-—ksin &+ 
do, (Rk? cos? 6, —2mw s/h)! 

+k sin Ox 
7 (Rk? cos? 0, — 2mw;/h)! 
XL cos 04F (k® cos? 0, — 2mow;s/h)*] 


sk sin 6, 








zk sin 0 
_ (Rk? cos? 0, —2mw s/h)! a 





Thus we have, finally, 


Si=(1/k) { oe | 2xe—qz|)(dz/z) 
=[k|2re—q,|}"', (41) 
so that 
E\’* =[whX /4k?mw(1—x,)BN|2x2—q,| ] 
X[(242e—qy)-a;P exp (—2Wy"). (42) 


Carrying through the same procedure in the 
case of phonon absorption by the neutron, we 
obtain the result expected directly from (35): 


Ev j=x,Ey!". (43) 


Now (42) represents the total polycrystalline 
scattering cross section for excitation of the fth 
oscillator, (43) for absorption of energy from the 
fth oscillator. To obtain a physically measurable 
quantity, it is necessary to sum (42) and (43) 
over all f—for each qy, over the three directions 
of polarization and the two (y and y) oscillators, 
and then over all of q space (the restriction to 
half of q space having been refnoved at (34))— 
and then to add the results. Also, summation 
over suitable « must be carried out. The suitable 
values of « and the corresponding ranges of q 
(the index f is no longer required) are ascertained 
by the conditions 


k—k’=2r7-—q (44) 


and conservation of energy, (39). For phonon 
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absorption, the latter is replaced by setting —w 
for w, so that the general condition for energy 
conservation is 


k? =k? +2mw/h. (45) 


From (44) we find, upon squaring, the condi- 
tion for the cosine of the angle between k and 
(2r2—q) to be less than unity: 


kt —R8 + |2ee—ql? - 
2k | 2r2—q| 
(|2r2e—q| —k)?—k” <0. 





or 


This inequality is satisfied if and only if 
k’—k<& |2re—q| <R+#’, 


a result significant for phonon absorption, where 
k’>k. The condition that the cosine of the angle 
between k’ and (q—22) be less than unity 
yields the corresponding result applicable to 
phonon absorption: The sign of the left-hand 
member of the above inequality is reversed. 
Thus the general restriction for all first-order 
inelastic processes is 


|k—k’| < |2re—q| k+R’. (46) 


The previous introduction of an average sound 
velocity renders summation of E;“* over the 
three directions of polarization quite a simple 
matter. Since for a given q the three ay are unit 
orthogonal vectors, the sum over them gives, 
from (42), 


Ey’? =[whx /2k*mw(1—x)BN] 
X |2xr2e-—q| exp (—2W"), (47) 


where the prime of E,’* indicates that summa- 
tions over the three directions of polarization 
and the two oscillators (y and yu): have been 
carried out for the particular q. 

In summing (47) over all points of q space, 
we again integrate after multiplying the sum- 
mand by the factor 


LN/(2)* }dgqidgedqs 
=[N/(2r)*]B¢@ sin 0,dqd0d¢y., (48) 


whereby the relation (22) has been written in 
polar coordinates. Instead, however, of inte- 
grating over the angle 6,, we introduce the 
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variable \, where 
|2me—q| =2ar+Aq (—1<A<1). (49) 


For r#0 (the case r=0 is handled separately 
later), it can be shown that for simple, body- 
centered and face-centered cubic lattices 247 > qo, 
so that the substitution (49) is valid for these 
structures at least. If the direction of + is taken 
to be that of 6,=0, then (49) becomes on 
squaring 


4n’r’?+¢° —4rrq cos 0,=49'* 7? +4r7rqi+¢_7X’, 
whence 


sin 0,40,=(247r+Aq)dXd/2ar. 


Therefore, finally, the total inelastic cross 
section (for phonon emission) is given by 


Ey’ =(Xh/322°k*mr) 
xf f (2xr+)Aq)? exp (—2W’) 


w(1—x) 
As variable of integration, \ takes on all values 
between —1 and +1 consistent with (46), 
which through (49) now reads: 


|k—k’| <2ertrqc kth’. (50) 


gdgqdnd ¢. 





Also, for given 1, q is limited to those values 
consistent with conservation of energy, the above 
inequalities, and g < go. Integration over g,merely 
yields a factor of 27. 

We employ an approximation to the first mean 
value theorem of the integral calculus by re- 
placing the factor, exp (—2W’), by its value for 
\=0, and thereby remove it to the left of the 
signs of integration; W’ is therefore the same as 
Wo", given by (30). Upon integration over ¢, 
and X, then, 


Ey’ =(Xh/16rk?mr) exp (—2W.") 
qdq 


w(1—<x) 





x f [4etr’h + 2ergr?-+-(1/3) ga"; , (51) 


where Aj, Az are the limits of integration over X. 
They are functions of r and gq, consistent with 
both (50) and —1<A,;SA2< 1. 

In order to evaluate the integral over g, we 
appeal once more to the Debye approximation— 


namely, w=cqg, where c is the constant sound 
velocity. By letting B=hw/koT, whence x 
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tain 
E'x1 =(AXke?T?/164rh?k?mc) exp (—2W%") 


X{4e2etKy +(QerkoT/hc)Ky. 
+(1/3)(keT/he)*Ky }, (52) 
where 
Kt f arn 8" Lexp (8) —1]}-"'+1}de 


honon emission 
(p ) (53) 


K,*= f (A2”—A1”)8’Lexp (8) —1 }-"'dB 


\ (phonon absorption), 


for v=1, 2, 3. 





In order to interpret (50) to obtain the limits 
of integration in (53) as well as the choices of A; 
and 2, we must rewrite the energy conservation 
condition (45) in terms of the Debye approxi- 
mation : 


k’? =k? +2mocq/h. (54) 


The lower limit is in general zero; all limits 
which differ from zero are inversely proportional 
to T. (For, B=hcg/koT, and the limits are 
determined through restrictions upon g alone.) 
If the upper limit for g is just go, then the upper 
limit for B is @/T, since hcgo/ko= 9, the Debye 
temperature of the scatterer. For phonon emis- 
sion, the highest possible value for qg is clearly 
(hk?/2moc), by (54); the corresponding limit for 
B is (h?k?/2mokoT)=(To/T), where T> is the 
neutron temperature. A convenient aid to ob- 
taining the limits and values of \; and Az will be 
indicated in the final section of the paper. 


C. The Case +=0 


We now investigate the possibility of inelastic 
scattering for r=0. For phonon emission, by 
(46) and (54), 


k?—2gk+¢@°< k’? =k? —2mocg/h, 
whence 
hq+2m)(c—v) g 0, 


where v=(hk/my) is the incident neutron ve- 
locity. Hence, for neutrons of velocity less than 
the sound velocity in the scattering substance, 
single phonon emission cannot occur at all for 
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=exp (—/hw/koT) =exp (—8), we at length ob- 7=0. For faster neutrons, (v>c), only oscillators 


for which g < 2mo(v—c)/h may participate. 
For phonon absorption, (46) and (54) give 


k?+2kg+¢@’ > k’? =k? +-2mocq/h, 


whence 

hq > 2mo(c—v). 
These results, obtained in similar fashion by 
Wick,’ are only approximate in that they rest 
upon the Debye approximation for the frequency 


spectrum. 
For r=0, the transformation (49) is neither 


valid nor necessary. We employ (48) directly to 
obtain 


E41=(hX/4k*m) 
x ftx°7 exp (—2W”)/w(1—x) ]g’dq, (55) 


where 
W® = (3h?q?/2mko®) { (1/4) +(T/0)?Q1(9/T)}, 


and the limits are given through |k—k’| <q 
¢k+k’. The evaluation of this integral may be 
carried out numerically for those specific prob- 
lems in which consideration of 7 =0 is required. 

(It turns out that for Fe scattering 300- 
degree neutrons, the conditions necessary are 
not satisfied, so that r=0 does not enter into 
the computation carried out for this work.) 


D. Recapitulation 
The total inelastic cross section is obtained, 


finally, by summing the expression Ein=E\+E-1, 
(52) and (55), over all suitable values of ¢: 


Ein= Xr Ein. (56) 


The suitability of a particular « depends upon 
the consistency of its absolute value r with the 
inequalities (50) for any simultaneously per- 
missible } and g. Also, there are possible re- 
strictions upon the values of + imposed by the 
crystal structure of the scattering substance. 
For example, as we shall see later, a body- 
centered cubic lattice requires that (ar), where 
(a) is the lattice constant, be an even integer. 
For a simple cubic structure there is no such 
restriction. 


7G. C. Wick. Physik. Zeits. 38, 403 (1937). 
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III. COMPARISON WITH X-RAY PROBLEM 


From (54) we see that in the limit mo<(hk?/ 
2cqo) = (h?k®/2ko®)—which is equivalent to 9 
<T, the neutron temperature—we have 
k=k’. In this limit, therefore, we should expect 
that our results should apply to the scattering 
of x-rays as well, provided that isotropy of 
scattering is assumed. Insofar as the mathe- 
matical treatment is concerned, setting k=k’ 
through m)=0 merely means lifting the restric- 
tions on the limits of integration in K,—+v1z., the 
limits are always 0 and (@/T), while \1= —)z 
= —1. This comes about through the bp sgene 
of the condition (50) to read 2xr+ Ag 2k. (It is 
clear, then, that for values of 7 so large that 
2xr+qo>2k, there are still restrictions on the 
values of the upper limits and of A». But it is 
not necessary to consider such special cases here; 
they would of course apply to x-rays as well.) 
For the sake of simplicity, only the term in (52) 
involving K,*! is considered; but the following 
comparison is valid even if the smaller term in 
K;*! is taken into account. (Since 4; = —A.= —1, 
K;*'=0.) 

From (52), (53), and (25), under the conditions 
of the preceding paragraph, 


Ein=E, +E 1 =(XrkeT?7/hek*m) 
Xexp (—2Wo") {0:1(0/T) + (1/4)(0/T)?*}. 


Since there is no comparable formula for x-rays 
readily available, it is necessary to develop 
further the work of Zachariasen,? beginning with 
his Eq. (15), which we rewrite in the notation of 
the present work: 


J2=(X/4r) exp (—2Wo")(242)?(B/m) 


xf" om Q 3} sin® [gi(ki—k,’ +9:)/2] a 
w? It sin? [(k:—k,’ + 9;)/2] (2x) 





Here we have written the angular independent 
(X/4mr) for the non-isotropic Sf*. Also, (k—k’)? 
has been replaced by the equivalent (24*—q)?; 
the subsequent negligence of q compared with 
2re is equivalent in one case to ignoring of the 
K; term, while in exp (—2Wo") it corresponds 
with the use of the mean value theorem in the 
present work to reach the result (51). The 
frequency » is replaced by (w/2r), while the 
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volume element dv is substituted for by 
[4aq*dq/(2xr)*]}—integration over a sphere of 
radius go thus replacing the integral over rec- 
tangular variables. An average sound velocity 
has been adopted, so that the sum >>; y7=1 is 
carried out. In terms of w, 


Q=hwLexp (hw/koT) —1]}-'+ dhw. 


Reducing Zachariasen’s formula to the single 
atom by dividing by N, we obtain 


(J2/N) =(4x°X/m) exp (—2W,") 
xf ‘(Q/w*)8(k —k’ +q—20e)¢¢dg, 
0 


where use has been made of (27) to replace the 
sine-squared ratios. Instead of integrating over 
g immediately, we perform the angle integrations 
over all directions of k and k’ to obtain the total 
polycrystal cross section. These integrations of 
the function 6(k—k’+q— 27) have already been 
carried out above—(38) to (41). The result for 
the x-ray inelastic polycrystal cross section per 
atom (again with negligence of q compared 
with 27) is then 


(J3/N) =(Xa7r/mk?) 
Xexp (—2009) [ '(Q/w*)qtdg. 


Using w=cqg, and introducing B=hcg/koT, we 
finally arrive at the expected result: 


(J3/N) = (X rk? T?*r/hek?m) exp (—2Wo") 
X {Q:1(0/T) +(1/4)(0/T)*}. 


This is precisely the result directly above for 
Ein, and thus provides the desired check with 
the x-ray theory as developed by Zachariasen. 

(Even though, on the one hand, the term in 
Ks; has been dropped and, on the other, q has 
been neglected in comparison with 2, the 
check is still perfect. Detailed comparison with 
the steps leading to (52) readily verifies the 
equivalence of these negligences, and conse- 
quently the exactness of the agreement.) 


IV. RESULTS FOR LOW TEMPERATURES 


For TKO and TKT)y= (h*k?/2mko), both, the 
upper limit of integration in K,, (53), is large 
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compared with unity. Since the term of the 
integrand involving [exp (8)—1]" is close to 
zero for B>>1, the upper limit may be replaced 
by infinity for such small values of 7. For the 
term not involving the exponential, which ap- 
pears in K,! only, the limit must remain finite. 
[ Actually, the criterion for the smallness of T 
varies with 7 since the limits of integration 
depend strongly on this parameter through (50). | 
Also, for T small, it is only the range of small g— 
q<koT/hc—which contributes appreciably to the 
result from the exponential part of the integrand. 
For sufficiently small g, the restrictions upon A 
through (50) drop out, so that \;= —A:=—1. 
Comparison of (53) with (25) shows that for 
this case K;'=20,(©) and K,;'=2Q,(«) 
+f (A2—A,)8dB. For the term in K,! not in- 
volving the exponential, this simplification can- 
,not be made. For 7<0O, W,’ is simplified by 
setting Q,(«)(=2*/6) for Q:,(0/T), (25). 

Because the term of K,! does not involve the 
exponential, it is not possible to perform a 
complete simplification for T small. One can, 
however, glean some general qualitative points 
of interest from the preceding results. Of course, 
for T=0 phonon absorption must vanish, as the 
factor T? in (52) shows, since the integrals K,~ 
converge for T=0. In the case of phonon 
emission, however, since the integral K,' contains 
a term which diverges as 7~’—', inspection of 
(52) reveals a non-vanishing cross section at 
absolute zero. For, although at 7=0 the lattice 
has no energy to give to the neutron, it is still 
possible for an elastic vibration to be excited by 
the impact of the incident neutron. 

For TKO, T&T», both, the most significant 
term, aside from that which is non-vanishing at 
T=0, is proportional to 7*. This quadratic 
term turns out to be directly evaluable, so long 
as the condition 


(39°h?7r?/mko@) <1 (57) 
is met. From the remarks of the preceding 
paragraphs, it is seen that for low temperatures 
Ein = Ey +E_1= Col 4X ho?r/4hek*m) 

X {401( ©) 72+ C1 — Ci(12*h?r?/mkoO) 
X(7/0)?Qi1(%)}, (58) 


where 
C=7 f Os—r)ade 
+(kaT"/2urhe) [ (dst —rs1 648 
+(*/3)(bo/2arhe)? f (str) 84, 


C:= exp ( _ 3h?x*7?/mko@) ’ 


exp (—2W ") has been replaced by the first two 
terms of its expansion. Now, the highest possible 
value of C, is attained for \;=—A:=—1 and 
when the limits of integration are zero and 


(Q0/T). This maximum is then 
Ci = @?{14+4(ko@/2mrhc)?}. 


It can be shown that the second term in brackets 
is at most of the order of magnitude unity; so 
that C,= 0? is a satisfactory order of magnitude 
approximation. The value of the last term in 
brackets in (58) is therefore, to order of magni- 
tude, (129°h?r*/mk,@)Q:1(©)T?. If (57) is ful- 
filled, then this term is negligible compared with 
the first in brackets, and (58) becomes: 


Ein=C2(4X ko? Qi( ©) /hek*m)T?+C;, (59) 
where 
C3 = CiCo(a#Xko?r/4hck?m). 


It is interesting to compare this result with the 
corresponding quantity for elastic scattering. 
Expansion of the exponential in (31) to two 
terms and use of the fact that O= (hcgo/ko) 
= (hc/ko)(6x?/B)*, by (23), gives 


Ex, = Ci—Co( 4X ko? Qi( © )/hetk?m)T? 
(C,=constant). 


With (59), this gives the noteworthy result that 
for low T [and (57) ], the total (elastic plus ine- 
lastic) cross section is temperature independent. 

Now the values of w contributing appreciably 
to the integral Q:( ©) are small—the smaller, the 
smaller T. Since it is this integral which provides 
the term in 7*, we need to consider processes 
involving only low frequency oscillators in the 
consideration of this term. For finite tempera- 
tures these are on the average excited to high 
quantum states to render their energies to the 
order of ko7. For high quantum numbers— 
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whereby the separation of energy levels is small 
compared with ko7—it should be possible to 
give the problem a classical treatment. We 
should expect, therefore, to find a classical de- 
scription of the constancy of total scattering 
with temperature, under the assumption (57), 
for low temperatures. 

Since for low frequency oscillators k—k’ = 27+ 
approximately, 47*h?r? is essentially the square 
of the neutron’s momentum change in being 
scattered. If we consider now a classical collision 
between a neutron and a free scattering nucleus, 
4h??? is therefore the square of the scatterer’s 
momentum change. Also, mkoQ@ is of the order of 
the square of the momentum of the scattering 
nucleus. Condition (57) holds, therefore, so long 
as the motion of the scatterer is hardly disturbed 
by the impact of the neutron. This condition 
obtains if and only if the scatterer mass is large 
compared with that of the neutron. Since we 
have required that 77>, or that the probable 
neutron energy change be small compared with 
the incident neutron energy, it is necessary also 
that the scatterer velocity be small compared 
with the incident particle velocity. Otherwise, 
the neutron would experience a considerable 
“Doppler shift’’ in velocity and therefore undergo 
an appreciable energy change. 

Thus the classical picture is the following: 
Light, rapid particles are scattered by heavy, 
slowly moving centers. All processes, in general, 
involve an interchange, however small, of energy ; 
but owing to the relative mass and velocity 
magnitudes, the scattering (total for all energy 
changes) is just what it would be for neutrons 
scattered by immovable centers at rest. The 
more strongly are (57) and T7<Tp> fulfilled, the 
more closely does the total scattering approach 
the “rigid scattering.” This is equivalent to our 
quantum-mechanical limiting result. 


V. RESULTS FOR HIGH TEMPERATURES 


For T>T>, T>9, both, then the integrals 
are simplified somewhat by dint of the fact that 
8<1 over the range of integration, and exp (8) 
may be replaced by (1+8+ 2/2). Thus we have, 
by (53), 


ns f (de Ay") (81-4 B"/2)d8 


to this approximation and with use of (1+2z)— 
=1—z for 21. As for the factor exp (—2W,’), 
for large T, Q1(0/T)=(0/T)+(1/4)(0/T)?, so 
that, according to (30), 


W.' =(62r°h?7?/mko@) (T/®@). 
It is thus clear that for large T the cross section 


E's varies as (aT +b) exp (—aT), where a, b, a 
are constants. (This fact is hardly evidence that 
the inelastic scattering cross section vanishes at 
very high temperatures; for temperatures too 
high, the entire mathematical treatment breaks 
down, since the assumption of harmonically 
bound nuclei in a crystal is no longer valid.) 

The particular asymptotic expressions for both 
high and low T are employed in the computation 
carried out at the end of this work. 


VI. NUMERICAL EVALUATION 


The remainder of the work is devoted to the 
computation of the inelastic cross section per 
nucleus for 300-degree neutrons scattered by 
polycrystalline iron at temperatures up to 
1000°K. The slight change in lattice constant 
with temperature is neglected in the computation. 

A description of the positions of the cubic 
body-centered lattice points—the structure ex- 
hibited by iron through 1000°K—in terms of the 
basic vectors a, @, a; must be such that the 
position vector @=018:+0282+03a3 of the eth 
lattice point bear only integer components 
71, 2, o3. Also, every integer triple o1, o2, o3 
(0<o:& gi) should denote the position of one 
and only one lattice point. These requirements 
follow from the work of Section II. If ir, ico, ics 
are three mutually orthogonal dimensionless unit 
vectors, then the following triple a:, a2, a; 
satisfies the required conditions for the body- 
centered cubic lattice: 


a3= (a/2) (icy + ize +iz;), 


where (a) is the lattice constant. 
The reciprocal lattice is defined generally by 
the system of basic vectors b;, be, bs, where 


bi=a:Xa;/D, be=a;Xai/D, bs=aiXa:2/D; 
(60) 
D=(ai Xa) *@s. 
For the present choice of ai, a2, a3, we find 
bi =(1/a)(iz:—iz3), be=(1/a)(ir.—izs), 
b3 = (2/a)izs. 


a,= Giz}, ae= ize, 
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Fic. 1. Inelastic scattering cross section vs. temperature 
of scattering crystal. 


A vector * having components 7), 72, 73 in the 
reciprocal lattice system is written: 


c= tibitrebs+rsbs = (r:/a) (iz —izs) 
+ (r2/a) (ize —izs) +(273/a)izs. 


This gives the relation between the Cartesian 
and reciprocal lattice (covariant) components of 
a vector: 


Tz1=7;/a, Tr. = T2/d, tz3 = (1/a)(273—71—72). 


From here, the Jacobian of the transformation 
from the Cartesian to the reciprocal lattice 
system is found to be (2/a*). Consequently, its 
reciprocal, defined by (22), is B=(a*/2) ; whence, 
according to (23), go=(12x*)!/a. 

Also, the magnitude of a vector ¢ in terms of 
its reciprocal lattice components is clearly given 
through 


a?r? = 7;?7+ 72? + (273—71—7T2)?. (61) 


Clearly, for integer triples 71, r2, 73 we must have 
a*7? an even number. Hence in summation over 
2, as in (56), (32), only those + for which a?r? is 
even are summed over. (Physically, this means 
a destructive interference between the two simple 
cubic lattices making up the body-centered 
lattice, for “‘Bragg reflections” at “Miller planes” 
for which a*r?—the sum of the squares of the 
Miller indices—is odd.) 

It is noticed for the special case of a body- 
centered lattice, the Jacobian B is equal to the 
volume per nucleus, here a*/2. It is possible to 
show that this result is perfectly general provided 


ROBERT WEINSTOCK 





that the basic lattice vector system is chosen 
with the “‘one-to-one”’ restriction imposed here, 
In this event, (a:Xa@)-as is the volume per 
nucleus. By direct substitution through (60), we 
find that (b, Xb2)-bs=[ (ai X a2) -as}, the num. 
ber of nuclei per unit volume. Now if we write 
the expansion of b; in terms of the unit Car. 
tesian vectors as b;= )); Bi;iz;, then += >); rb, 
= >i; rBisiz;. Hence r2;= >>; 7:83; the Jacobian 
of this transformation is the determinant whose 
elements are the coefficients 8;;. This is clearly 
equal to (b; Xbz)-bs, from the definition of the 
8:3. Since this determinant is the reciprocal of 
the earlier defined B, we have B=(a;Xae)-a,, 
whence the theorem is proved. 

The numerical work can be greatly simplified 
through consideration of the expression (51) for 
Ey". Even for g=qo, the term (1/3)q@*A\* can be 
shown to be small compared with 47*7*). First, 
since —1<¢A<1, |A|*< |A|. For 70, the smallest 
employed value of r? [by the considerations just 
following (61) above] is (2/a*). Therefore, 
4n°7? > (82?/a*) =79/a*. On the other hand, 
go? = (12?)'/a?=24/a*?, so that (1/3)g?< 8/a’, 
Hence it is quite reasonable to neglect the term 
involving K;*' in (52) compared with that con- 
taining K,+'. The larger r—and the lower T— 
the better the approximation; in the present 
work the error thus involved is much less than 
10 percent, since an appreciable contribution to 
the results comes from values of r? a number of 
times as large as (2/a?). Similarly, the term in- 
volving K;*! is neglected, since for the most 
part (A2?—A,”) =0. Thus simplified, the expression 





















for E+1 becomes 
E's1=X[wa(ar)ko?T?/4hm(ak)**] 
+1 
Xexp (—2W,o")Ki , 







where, for sake of convenience, the reciprocal 
lengths k and 7 have been introduced as multiples 
of (1/a). 

The author has, with great benefit, plotted the 
curves |k-k’| as functions of g—all, by (54), 
parabolic arcs. There is one plot for each of 
phonon emission and phonon absorption, each 
done in units of (1/a). On the same plots were 
drawn the lines 2x7 for all permissible values 
of r. This leads to a simple graphical interpreta- 
tion of the inequalities (50): The double (A, q) 
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integration is performed over that portion of a 
“\—q plane” consisting of the area common to 
the region between the two parabolic arcs 
|kk’| and that between the two straight lines 
Ier+q—and to the left of g=gqo—for each value 
of r. The lines 2mr+g, when ‘“‘integration 
boundaries,”’ are the curves \= +1. If the para- 
bolic arc |k+k’| is an integration boundary, 
then it corresponds to the curve \=(|k+k’| 
—2zxr)/q. Intersections between the straight lines 
and the parabolic arcs provide ready information 
as to the limits of the g integrations (and there- 
fore of those over 8), and as to when one changes 
from one to another of the values of A; and Az. 
The g values of these intersections are obtained 
exactly by the solution of four quadratic equa- 
tions—2mrr+q = | kk’ | —but their meanings are 
most readily interpreted through the graphs. 
Whenever an integration boundary is sup- 
plied by a parabolic arc, the latter has been 
replaced by a suitably chosen tangent straight 
line. Since the integrations are carried out 
numerically, replacing of the parabolic arcs by 
straight line segments serves the great advantage 
of omission of the temperature T from explicit 
appearance in the integrands. Thus only single 
tabulations, good for all 7, need to be carried 
out; the temperature appears only in the limits 
of integration. The only integrals so required are 


Q(z) = f ‘eLexp (8) —1]}"dg, 


Quis)= f tLexp (@)—1}+1}48 
=log [exp (z)—1]; 


in addition are required 


fe {Lexp (8) —1]}-'+1}d8 =Q,(z) +2°/2, 
0 


[Cex @ -17748 = Q4(2)-s 


In the high temperature approximation con- 
sistent with Section V, we have for 1 


Qi(2) =2—(2*/4); Qo=log s+ (2/2) — (2/8). 


The constants from which the computations 
spring follow: h=1.056X10-*’ erg sec., ko=1.38 
X10-'* erg deg.—', mo=1.66X10-** g, m=9.22 
X10-** g, a=2.90K10-§ cm, @=453°K, wo 
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= (ko@/h) =5.93 X10" sec, go=[(122*)!/a] 
=1.69X10® cm, c=(wo/go)=3.50K10® cm 
sec.—!, T)=300°K, v= (2ko 7 o/mo)' = 2.23 X 10° cm 
sec.—!, k=(mov/h) =3.51X10*® cm=, ago=4.91, 
(hk?a/2moc) = 3.24, ak = 10.18. 

The results of the computation may be grasped 
most readily from the accompanying graphs 
(Fig. 1). From 0.6 percent of the free elastic 
cross section, at absolute zero, the inelastic 
cross section rises steadily to 19.2 percent, at 
1000°K. Between 150°K and 400°K, the increase 
of the inelastic cross section is essentially linear, 
but at higher temperatures the slope exhibits a 
steady decrease. 

Use of the low temperature approximation 
given earlier produces the stated parabolic rise 
from T=0; extension of the approximate curve 
to 250°K shows good agreement with the rigor- 
ously computed curve at 150°K, but at higher 
temperatures the deviation becomes appreciable. 
Since the low temperature expression is valid for 
TKT>,—here 300° K—the beginning of deviation 
at 150°K is to be expected. 

The asymptotic expression valid for high 
temperatures gives agreement at 1000°K to 
about the computational accuracy of the work— 
about 1 percent of the result. Although the 
approximation made is valid only for T>0,— 
453°K for Fe—the asymptotic expression gives 
an accurate result even below 400°K; the devia- 
tion at 250°K is only about 8 percent of the 
result. This low temperature agreement is merely 
fortuitous, however; it arises from the compen- 
sating effects of the errors committed in the 
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Fic. 2. Comparison of total and elastic scattering 
cross sections. 
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approximation. Although there is no @ priori 
assurance of the over-all reliability of the high 
temperature expression in general, one is led to 
expect a suitable computational result for any 
given case by means of a plausible interpolation 
between the low and high temperature curves. 
The elastic and total (elastic plus inelastic) 
scattering cross sections are plotted in Fig. 2. 
Since (32°h?r?/mk,@)<K1 here (<=.007a?r?), the 
predicted constancy of the total cross section at 
low temperatures is clearly exhibited. The con- 
stancy is strict up to 150°K; but even up to 
400°K, the entire variation (a monotonic de- 
crease) is no more than a percent. At 1000°K, 
the value has decreased to 90 percent of the 
free elastic cross section, from the maximum 
of 97 percent at absolute zero. Such a variation 
should be readily detected by a reasonably 
sensitive experiment. 


G. GAMOW 
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By use of the method of successive homology transfor- 
mations in conjunction with Bialobjeski-Eddington’s 
formula, it is possible to predict the changes of radii and 
luminosities of massive stars in different stages of their 
hydrogen and contractive evolution. It is shown that at a 
certain stage of contraction the star must reach the 
maximum of its luminosity, and it is calculated that for 
the star masses of 5, 10, 20, and 40 suns the maximum 
luminosities are 1X 10°, 4105, 110°, and 4X10® suns, 
respectively. In the later contractive stages total luminosity 
of the star remains constant whereas its visual luminosity 
rapidly decreases because of the shift of energy into the 
ultraviolet part of the spectrum. It is to be expected that 
for such high values of the (luminosity)/(mass) ratio, 
radiation pressure becomes strong enough to eject stellar 
atmospheres into the surrounding space, and it is shown 
that the ejection will actually take place if the force of 
gravity on the stellar surface will be somewhat reduced by 
the centrifugal force due to axial rotation. It is also shown 
that under the conditions existing in the ejective atmos- 
pheres of Wolf-Rayet stars, radiation pressure is pri- 
marily due to the light scattering by free electrons. In 


1. INTRODUCTION 


T is generally accepted at present that during 
the largest part of their evolutionary life stars 


(Received October 11, 1943) 





discussing the motion of the ejected gases, it may be 
necessary to assume that their original velocity of about 
2000 km/sec. can be considerably reduced by the gravi- 
tational action of the star. If this is the case, gaseous 
envelopes which will form around Wolf-Rayet stars in the 
course of several centuries will possess properties very 
similar to those of the planetary nebulae. This would 
indicate a close evolutionary relationship between these 
two classes of celestial objects. An alternate evolutionary 
road of a massive star consists in the formation of an 
energy-producing shell, which will take place in all cases 
where the convective currents due to axial rotation are not 
fast enough to secure homogeneity of stellar matter. It is 
indicated that the growth of such shells will probably lead 
to the formation of an extensive atmosphere which offers 
certain possibilities for the interpretation of the so-called 
red giant stars as intermediate evolutionary stages between 
the stars of the main sequence and the Wolf-Rayet stars. 
In conclusion the problem of stellar collapse, which is 
expected to take place towards the end of contractive 
evolution, is discussed in some detail. 





receive their energy supply from thermonuclear 


reactions in which hydrogen is being transformed 
into helium through the ‘catalytic’ action of 
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carbon and nitrogen nuclei (C—N cycle).' In the 
case of the sun, for which the central tempera- 
ture, central density, and the hydrogen content 
are determined by the direct numerical integra- 
tion of the fundamental equilibrium equations, 
the calculated rate of energy liberation of C—N 
cycle is in excellent agreement with the observed 
value of luminosity. 

By assuming the sun as the starting point, and 
using the method of homology transformations, 
it is possible to calculate the observable charac- 
teristics of stars with various masses and various 
H contents. The fundamental equations of 
stellar equilibrium can be written in the well- 


known form 2 


dp/dr= —(GM,p/r’), (1a) 
dM,/dr=4rnr*p, (1b) 
pb=(R/uB)eT, (1c) 
L/4ar*= —(co/3kp)(dT*/dr), (k=kop*T-*) (1d) 
L=A:x-ro'?p2T.". (le) 


Here p, p, and T represent the density, total 
pressure, and temperature in stellar interior 
(index ¢ pertaining to the center) ; u, k, and 1—8, 
the mean molecular weight, coefficient of opacity, 
and relative radiation pressure, respectively; x 
the H content, 79 the radius of the small central 
zone where the energy is mainly produced, and 
n the exponent in the formula for the tempera- 
ture dependence of the thermonuclear reaction. 
(For C—N cycle at solar central temperature 
n=18, decreasing as 7.~' in the case of hotter 
stars.) The value of 8 can be assumed to be 
approximately constant through the body of any 
given star, being entirely determined by the 
mass of the star through the well-known Bialob- 
jesky-Eddington equation :* 


[M/M(O©) Ju?=18.0-((1—6)/6*]*. (2) 


Applying the method of homology transforma- 
tions to the equations [(1a)—(le)], we can 


1C, v. Weizsacker, Physik. Zeits. 39, 633 (1938); H. 
Bethe, Phys. Rev. 55, 434 (1939). 

? Compare, for example, G. Gamow, Phys. Rev. 53, 595 
(1938) ; 55, 718 (1939). 

’ This equation was primarily developed for the so-called 
standard model. However, it also gives very good results 
in application to the point-source and contractive models. 
Compare: S. Chandrasekhar, Introduction to the Study of 
Stellar Structure (Chicago University Press, 1939). 
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express the relative changes of stellar radius R, 
luminosity LZ, and central temperature 7. 
through the changes of M, wu, ko, and x. We 
find : 


Rw~ Mete-s—D/ (nt+3a—s+3) (uB) (n—s—4) /(n+3a—s8+3) 
X Roll (nt8a—st3) x 1/ (n+-3a—s+3) (3a) 
L ~~ M1GB+2a)n+8a+s+9) | (n+3a—8+3) 
x (u8) [(44+3a) n+38+12] /(n+3a—s+3) 


x Rg Ot 81 (nt8a—s+3) y—(e— 3a) | (n+ 38a—s+3) | (3b) 
T.~ MuB/R. (3c) 


For temperatures between 1X107°C and 
5 X10’ °C, and densities between 1 and 100 g/cm? 
(i.e., within the range of temperature and density 
variations for the normal stars of the main 
sequence) the changes in opacity can be fairly 
well represented by the exponent values a=0.5 
and s=2.75 (Bethe’s approximation), and our 
expressions become: 


Rw~ M(-3.25)/(a+1.75) (uB) (n—6.75)/(n+1.75) 


K Rol oF 1.78) gol (n+1.78) | (3’a) 


L~ MGnt18.26)1(n+1.75) (113) (5-5n+20.25)/(n+1.75) 


XK bon 8)1(m+1.78) (1.28) +1.78) (3b) 
T.~ MuB/R. (3'c) 


In order to calculate the values of R and L 
for stars of different masses in the early stages of 
their H evolution, we assume that all stars 
possess originally the same hydrogen content 
x=0.35.4 To take into account the variability of 
the exponent m which depends on the central 
temperature, we carry our calculations step by 
step for the star masses M(©),5M(©),10M(O), 
20M(©), and 40M(©), using new values of m in 
each successive step. The values of 8 for the 
masses used are given by Eq. (2) (with ~=1) and 
are, respectively : 1.00, 0.95, 0.85, 0.70, and 0.55. 
As the results of such calculations, we obtain for 
log [R/R(©)]: +0.52, +0.70, +0.86, and +1.02, 
and for log [L/L(©)]: +3.0, +3.9, +4.6 and 
+5.3. 


‘ This value corresponds to the sun which, because of its 
a low luminosity, could not have used but a 
small fraction of its original H content during a few billion 
years of its existence. 




















Re x*x 

+1 4 oo x 
M 
' | e A leg Me» 
” -05 0 +0.5 +1.0 +15 420 +25 
L 
+6) leg Lie) Py +e", 
#4 9 
+9 7 
0 - 
“i 4¢* all Le M 
3 Me 

-2 








“05 0 +05 +10 4-5 +9.0 +25 


Fic. 1. Dependence of radius and luminosity on the 
mass for the stars of the main sequence. Continuous 
curve—theory ; solid dots—observation ; crosses—the data 
for Trumpler stars. 


Thecalculated log [R/R(©) Jandlog[L/L(©) ] 
are plotted in Fig. 1 as the functions of log 
[M/M(©)] (continuous curves); the separate 
black dots represent the observed values for 
individual stars for which reliable data on M, R, 
and Z are at present available.6 A very good 
agreement between the observed and calculated 
values indicates that even for rather massive 
stars, the method of homology transformations 
in conjunction with Bialobjesky-Eddington’s 
relation gives rather precise results. We are thus 
justified in using this method for the calculation 
of evolutionary changes produced in massive 
stars by the change of hydrogen content and the 
subsequent gravitational contraction. 

It must be noticed in the conclusion of this 
section that there is known a special group of 
very luminous stars which do not seem to fit on 
the regular R—M and L—M curves.* Whereas 
the radii and luminosities of these so-called 
“Trumpler stars’’ seem to indicate normal 
masses of the order of 50M(©), mass estimates 


5G. P. Kuiper, Astrophys. J. 88, 472 (1938). 
* In Fig. 1 the positions of Trumpler stars are shown by 
crosses. 


G. GAMOW 


7R. J. Trumpler, Pub. Astronom. Soc. Pac. 37, 249 






carried out by Trumpler’ on the basis of the 
red-shift analysis lead to the values of several 
hundred sun masses. We shall also see in the 
next section that Trumpler’s mass values are 
inconsistent with the alternative hypothesis that 
these stars represent the later stages of con- 
tractive evolution. It seems, therefore, that the 
high mass values claimed for these stars do not 
correspond to reality and are due to the erroneous 
interpretation of the observed shift of spectral 
lines as a purely gravitational red shift. It may 
well be that part of the shift is due to the ordinary 
Doppler effect taking place in the gas masses 
ejected from the surfaces of these stars, and, in 
fact, such ejection can be very well expected for 
the stars of such extremely high luminosity. 










2. LUMINOSITY CHANGES DURING THE 
CONTRACTION 





Since the stars of the main sequence produce 
their energy by transforming hydrogen into 
helium, their H content must gradually decrease 
with the time. The change of hydrogen content 
x with the time is evidently given by: 


Ax = —Lmy/eM - Al, (4) 











where my is the mass of a hydrogen atom, and 
e the energy liberation per proton. For the C—N 
cycle e=1.0X10~-5 erg. Applying (4) to the sun, 
we find that during the two or three billion years 
of its existence it could not have consumed much 
more than about one percent of its original 
H content. Since, however, the ratio L/M in- 
creases rapidly with the mass of the star 
[(L/M)~M?-*], heavier stars must consume 
their hydrogen considerably faster. Thus, for 
example, Sirius A would consume all its hydrogen 
in just about 3X10° years, whereas for the star 
with the mass 10M(©) it would require only 
4X10’ years. This fact leads to the well-known 
difficulty concerning the necessity of the assump- 
tion that many stars have been formed at com- 
paratively recent epochs. On the other hand, we 
are also forced to the conclusion that there must 
exist a considerable number of massive stars 
which have by now completely exhausted their 
H supply and are at present in the state of 
ultimate contraction. 
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We will now try to describe the changes in the 
luminosities and diameters of the massive stars 
evolving from their original positions on the 
main sequence towards the later contractive 
stages, and we shall use for that purpose the set 
of homology transformations derived in the 
previous section. When 35 percent of stellar 
matter is transformed from H into He (the re- 
maining 65 percent being the so-called Russell 
mixture of heavier elements) the mean molecular 
weight u increases from 1.0 to 1.7. The value of 
ko, on the other hand, will remain practically 
unchanged since both H and He contribute but 
very little to the total opacity. Taking m= 15, as 
an average,® we find from (3’) that the radius 
luminosity and central temperature of any given 
star will be subject to the relations 


R~ (uB)9-52x9-06, 
L~ (uB)®x-0-07, 


T-~ (uB)9- 480-96, 


(Sa) 
(Sb) 
(Sc) 











We see that unless x becomes very small as 
compared with its original value, the changes of 
Rand L will be caused entirely by the variation 
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Fig. 2. Construction of evolutionary tracks for massive 
stars in the luminosity radius diagram. 


*For the temperature range involved n changes only 
from 16 to 14, 
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of uw; the representative point of the star in the 
R-—L diagram will move along the line: 
L~R*-10/0.82~ R®,. When most of H is trans- 
formed into He, uw rises from 1.0 to 1.7 and, 
according to (2), the values of 8 (for the same 
mass) must increase correspondingly. For the 
star masses used above, we find the new 6-values 
to be: 1.00, 0.78, 0.62, 0.47, and 0.35. The 
formulas (5a, b) give us now for A log [R/R(Q)]: 
+0.07, +0.05, +0.03, and +0.02, whereas for 
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Fic. 3. Morse’s curves for the coefficient of opacity 
at high densities and temperatures. 
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A log [L/L(©)] we get: +0.85, +0.55, +0.3, 
and +0.2. The central temperatures rise only 
by a few percent remaining generally in the 
region of applicability of Bethe’s approximation 
of the opacity law. In Fig. 2 we have plotted the 
main-sequence positions for different star masses 
in the frame of the R—L diagram; the parts of 
evolutionary tracks corresponding to the H 
consumption are marked as I, and the positions 
obtained for 1. =1.7 are indicated by empty dots. 

When the largest part of H is transformed into 
He, » becomes practically constant, and further 
changes of R and L are given by the x power in 
the expressions (5). The radius begins to decrease 
whereas the luminosity increases inversely with 
the radius. Corresponding parts of evolutionary 
tracks are marked as II in Fig. 2. When x be- 
comes very small, thermonuclear energy pro- 
duction gives place to purely gravitational energy 
liberation, and the star will slowly readjust itself 
from the point source to contractive model. This 
change of the model results in a slight increase 












































of the luminosity, which will be taken into 
account somewhat later. 

Because of the rapid decrease of the radius 
during this stage of evolution, the value of T, 
rises very appreciably (7.~R™), and we are 
soon brought into the region where the ordinary 
opacity formula is no more applicable. We see, 
in fact, from Fig. 3, representing the results of 
Morse’s® opacity calculations for high densities 
and temperatures, that near 7=1X108 °C the 
value of k suddenly becomes constant. This is 
due to the fact that for higher temperatures ab- 
sorption is mainly due to the Compton scatter- 
ing by free electrons. Although this change of 
the opacity law begins near the stellar center 
and then gradually spreads out as the con- 
traction proceeds, we can assume for the purpose 
of our calculation that this change takes place 
in the entire body of the star when about one- 
half of its mass is heated above 1 X 108 °C. Since, 
further, the temperature ir... .a’yer correspond- 
ing to M,=0.5M is about one-half of the central 
temperature, a reasonable approximation will be 
achieved by changing the model at the point 
where T7.=2X108 °C. Using the fundamental 
equilibrium equations (1) with k=const., we find 
easily that under such conditions the luminosity 
of the star becomes independent of its radius. 
Estimating the points for which T.=2 108 °C, 
. we find for the changes along these parts of the 
tracks II: A log [R/R(O)]=—0.76, —0.73, 
—0.68, and —0.67, and A log [L/L(©)]=+0.76, 
+0.73, +0.7, and +0.7. In Fig. 2 these points 
are indicated by asterisks; beyond them evolu- 
tionary tracks run parallel to the R axis. Thus 
we come to the conclusion that by the end of the 
contractive evolution every massive star comes to 
the state in which its total luminosity becomes con- 
stant.'° These maximum luminosities for the star 
masses 5M(©), 10M(©), 20M(O©), and 40M(O) 
are: log L/L(O) = +4.6, +5.2, +5.6, and +6.2. 

In order to get some idea of the accumulated 
errors in our method which uses three successive 
homology transformations, we can compare the 
obtained luminosity values with the direct 


® Ph. M. Morse, Astrophys. J. 92, 27 (1940). 
10This holds, of course, only for star masses above 
Chandrasekhar’s critical mass: 5.754°=22M(©). For 
lighter stars the electron gas in the interior becomes 
degenerate at a certain stage of contraction, the luminosity 
begins to decrease, and the star evolves into a white dwarf. 








24 G. GAMOW 





homology relation between L and M for k=const, 
models. Using again the first four equations (1) 
with k=const. we get: 


L~ M*s*. (7) 


With the above values of 8 (for 4=1.7) the re- 
lation (7) gives: 


Lio/Ls5=3.2; L20/Li0 = 2.6; Lo/L20=2.5, 


while the ratios of calculated luminosities are 
3.7, 2.7, and 3.7. Taking into account the ap- 
proximate character of the method used,. we can 
consider the agreement quite satisfactory. 

We can also compare our results with the exact 
numerical integrations of k=const. model which 
were carried out by Keen \ (for a complete 
point-source model) and by Henrich” (for a 
point-source model with a convective core). 
For u4=1.7 and the star masses used in the 
present article, Henrich’s calculation leads to 
log [L/L(©)]=4.7, 5.3, 5.8, and 6.3. A very 
good agreement between these exact values and 
those obtained by the method of successive 
homology transformations indicates again that 
our method can be used safely for the calculation 
of evolutionary changes of massive stars. 

It must be noticed here that, using the 
homology transformations based on the sun, we 
have been working all the time with the point- 
source model and that the integrations of 
Henrich, with which we compared our results, 
have also been carried out for that particular 
model. As it was already mentioned above, the 
use of the point-source model is permissible only 
up to a certain point on the evolutionary track 
II where nuclear energy generation gives place 
to purely gravitational contraction. Beyond that 
stage the star must be described by a contractive 
model (polytrope of the index 3), and the transi- 
tion between two different models cannot be in 
general followed by the methods of homology. 


It was, however, shown by Biermann" that for’ 


equal masses, radii, and opacities, the con- 
tractive model leads to luminosities which are 
about twice as high as in the case of point-source 
model (this is due to the more uniform distribu- 


it Ph. C. Keenan, Astrophys. J. 89, 499 (1939). 
#2. R. Henrich, Astrophys. J. (Sept. 1943). 
13. Biermann, Zeits. f. Astrophys. 3, 116 (1931). 
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Fic. 4. Comparison of calculated evolutionary curves with 
the results of observation. 
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tion of the energy sources in the star). Thus the 
later parts of evolutionary tracks must run 
somewhat higher than originally calculated 
(broken lines in Fig. 2), and the maximum lumi- 
nosities become approximately : [log (L/L(©) ]cor 
=5.0, 5.6, 6.1, and 6.6. Similar corrections must 
also be applied to the results of Henrich calcu- 
lations. 

Comparing our results with observational 
material, we find that there actually exist a 
number of stellar objects which seem to fit the 
predicted properties of stars in the late stages 
of their evolution. While the lower part of the 
main sequence is represented by a narrow band 
containing the stars whose M, R, and L values 
are in a good accord with C—N cycle assump- 
tion, the upper part of the sequence is bordered 
by stars with anomalously small radii, and 
luminosities which are seemingly too high for 
their masses. In Fig. 4 we give the comparison 
between the calculated evolutionary tracks and 
the observed data for some massive stars. We 
find that such stars as U Ophiuchi, u Scorpi, 
V Puppis, and both components of 29 Canis 
Majoris and A O Cassiopeiae seem to be normal 
H-consuming stars. On the other hand in Y 
Cygni we probably have a case of rather pro- 
gressed H evolution since its luminosity is 


definitely too high for its mass. In the same 
general region of R—L diagram lies the group 
of stars known as P Cygni stars. Although mass 
values are not available for any member of that 
group, we can conclude that their luminosities 
are abnormally high for their masses from the 
observations which indicate that all these stars 
are ejecting their atmospheres into the sur- 
rounding space; we shall see indeed in the fol- 
lowing sections that such ejection requires rather 
high (L/M) values. A still more typical group 
is formed by the so-called Wolf-Rayet stars 
which are generally similar to P Cygni stars with 
the difference that the atmospheric ejection 
process is considerably more vigorous. W-R stars 
lie considerably farther from the main-sequence 
line than P Cygni stars, and can be thus assumed 
to be in a later stage of their contractive evolu- 
tion. The best available data pertain to one W-R 
star which happens to be the component of a 
binary system (' ~*’—193576), thus permitting 
an estimate of its mass. According to Wilson," 
the mass of the W-R component is 9.7M(©), its 
radius ~3R(Q©), and its total luminosity about 
2X105L(©). Thus the star is about 20 times as 
luminous as could be expected from the ordinary 
mass-luminosity relation. Plotting Wilson’s data 
on the diagram of Fig. 4, we find that the repre- 
sentative point of this star is located very close 
to the calculated evolutionary track for the mass 


10M(©). 
3. SOME REMARKS ABOUT SHELL FORMATION 


In the considerations of the previous section 
it was tacitly assumed that chemical changes 
produced by thermonuclear reactions near the 
stellar center are propagated through the entire 
body of the star sufficiently fast to secure the 
homogeneity of its material in every given stage 
of the evolution. This is equivalent to the as- 
sumption of some convective currents circulating 
between stellar center and the surface with the 
period which is shorter than the rate of stellar 
evolution. It was shown by Eddington" that 
such currents must actually exist within the star, 
being produced by asymmetry resulting from the 
axial rotation. He estimated the period of these 
currents to be of the order of 10'* years which 


4Q. C. Wilson, Astrophys. J. 91, 379 (1940). 
16 A, Eddington, M. N. R. A. S. 90, 54 (1929). 
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Fic. 5. Possible evolutionary track for a star with a 
growing energy-producing shell. Vertical cross indicates 
the place beyond which isothermal and radiative equi- 
librium solutions cannot be fitted together according to 
Chandrasekhar and Schoenberg. 


is definitely too long to secure any mixing-up of 
stellar matter. It was, however, indicated by 
Randers'* that Eddington’s value of circulation 
period represents most probably an enormous 
overestimate since it does not take into account 
the possibility of turbulent motion within the 
star. If, instead of the ordinary coefficient of 
gaseous viscosity, used by Eddington, we take 
the so-called ‘‘eddy diffusivity’’ obtained from 
the observation of the turbulent gas motion on 
the solar surface, the value of the circulation 
period drops from 10'* to only 10 years! It is true, 
of course, that the turbulence inside of the stellar 
body should be expected to be considerably 
smaller than on the surface, so that the above 
value of ten years represents only a lower limit 
of the circulation period. Thus a more detailed 
study is necessary before we will be able to 
decide whether convective currents due to rota- 
tion can secure the equipartition of material in 
the stellar body. Should the period of these 
currents turn out to be shorter than say 10° or 
107 years, the considerations of the previous 
sections will hold. In the opposite case the star 
will develop a ‘‘shell source’’ of energy slowly 
“growing” in size as more and more hydrogen is 
being consumed. The stellar model corresponding 
to this case will consist of: (1) isothermal core 
with zero hydrogen content (u=1.7) deprived of 
any energy sources; (2) thin energy-producing 

'*G. Randers, Astrophys. J. 94, 109 (1941). The author 


is grateful to Dr. Randers for the discussion on energy and 
paraplyen transfer. 


shell with the temperature suitable for C—N 
cycle; (3) outer envelope with the original 
H content (u=1). Such a model was studied by 
Critchfield and the author,’ and more recently 
by Chandrasekhar and Henrich,'* and Chandra. 
sekhar and Schoenberg." The latter investigators 
discovered a very peculiar fact that the fitting of 
isothermal-core solution, with the radiative. 
equilibrium solution of the envelope becomes 
mathematically impossible as soon as the mags 
contained in the core exceeds ten percent of the 
total stellar mass. Since, due to the progressing 
thermonuclear reactions, the amount of the 
dehydrogenized material increases steadily, we 
can ask ourselves what happens when the masg 
of the growing core trespasses the permitted ten 
percent. It was suggested by Chandrasekhar and 
his collaborators that the impossibility of solution 
for larger core masses will result in a non-equi- 
librium transition which could be interpreted ag 
the explosion of the star. This point of view is, 
however, not quite satisfactory since there seem 
to be no other stable states of lower energy into 
which such transition can lead. It appears, in 
fact, that the difficulty with the upper limit of the 
core mass lies essentially not in the general sta- 
bility problem, but rather in the boundary 
conditions usually assumed on the surface of the 
star. The solutions for the envelope used by 
Chandrasekhar and his collaborators in their 
calculations were obtained under the assump- 
tion that both density and temperature vanish 
on the outer boundary of the star, so that the 
impossibility of obtaining the fit with the 
isothermal core in the interior means actually 
that the conditions pyurs=0 and Tyurr=0 can no 
longer be fulfilled simultaneously. This result 
becomes particularly clear if, instead of trying to 
fit two solutions in the middle of the star, we 
carry the integration straight forward from the 
center to the surface. If the mass of the core is 
sufficiently small, it will always be possible to 
choose initial conditions (i.e., p- and 7.) in such 
a way that both density and temperature vanish 
for certain finite value of stellar radius. If, how- 


17 Chas. Critchfield and G. Gamow, Astrophys. J. 8, 
244 (1939). 

18S. Chandrasekhar and L. R. Henrich, Astrophys. J. 
04, 525 (1941). 

19S. Chandrasekhar and M, Schoenberg, Astrophys. J. 
96, 161 (1942), 
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ever, the mass of the core becomes too large, the 
point where the surface conditions can be fulfilled 
moves into the infinity; this is in accord with 
Chandrasekhar’s result that the rate of radius 
changes becomes infinite when the mass of the 
core approaches the critical value of 10 percent. 
Although this would mean mathematically that 
the star expands into infinity, physically it is not 
necessarily so. We know, in fact, that the opacity 
formula of Kramers’ type is valid only for rather 
high temperatures and that for temperatures 
below several thousand degrees the opacity of 
gaseous material drops rather sharply due to the 
absence of sufficient ionization. It follows that 
the ordinary solution of stellar interior should be 
carried out only up to a certain non-vanishing 
temperature corresponding to the photospheric 
layer, and that beyond that limit the transparent 
chromospheric layer should be fitted on.*° This 
will lead to a star model of finite dimensions al- 
though its radius will be abnormally large and 
central condensation of material very high as com- 
pared with the stars of the main sequence. It is at 
once apparent that such a model can be very 
useful for the understanding of red giant and 
supergiant stars, which are located on the right- 
hand side of the upper part of the main sequence. 
After the energy-producing shell grows to such 
an extent that very little material is left in the 
outer envelope, the escape of radiation from the 
generating zone will exceed the energy produc- 
tion by nuclear reactions. Nuclear processes will 
stop for the lack of sufficiently high temperature, 
and the star, with all its hydrogen practically 
gone, will go over into the state of gravitational 
contraction. Finite stages of this contraction will 
be identical with those expected for the stars 
retaining their point-source structure. In Fig. 5 
we give the possible track of the star evolving 
according to the shell model, as compared with 
the normal point-source evolution studied in the 
previous section. 

It is very possible that both modes of evolution 
are taking place in the stellar world: fast rotating 
stars with vigorous internal convection may 
evolve according to the point-source model, 
whereas slowly rotating stars may develop an 
energy-producing shell, and pass on their way 


*° Compare: A. Eddington Internal Constitution of Stars 
(Cambridge University Press, 1926), Chap. IV, Section 67. 
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through the red giant stage. This point of view 
may gain some support from the observational 
fact that the upper part of the main sequence is 
formed almost exclusively by stars possessing 
abnormally high axial rotations. 


4. MASS EJECTION BY CONTRACTING STARS 


As we have already mentioned in previous 
sections, the striking property of massive stars 
in the later stages of their evolution lies in their 
ability to eject their atmospheres into the sur- 
rounding space. The velocity of ejection, esti- 
mated from Doppler broadening of emission 
lines, varies from comparatively low values of 
about 100 km/sec. for P Cygni stars up to 
several thousand km/sec. for the typical W-R 
stars. There can hardly be any doubt that the 
observed continuous ejection is due to the effect 
of radiation pressure which overbalances the 
forces of gravity near the surface of the star. We 
can form some idea about the material density 
in the region where the main acceleration of 
ejected masses takes place by comparing the 
initial momentum of radiation coming from the 
stellar interior with the momentum flux of the 
ejected masses. If L is the total luminosity of the 
star, and R* the mean radius of the layer where 
the primary absorption takes place, the mo- 
mentum flux of radiation will be 


L/4rR™-c. 
If, on the other hand, p* is the mean density in 


this region and V the observed ejection velocity, 
the flux of mechanical momentum is 


p*V-V=p*V?. 
According to the law of conservation of momen- 
tum we can write: 





Sp* V2. (8) 
4nrR*2¢ 


TABLE I. Acceleration of H and He* due to bound-free 
transitions for p=10-". 














T 25,000°C 50,000°C 100,000°C + —_ 200,000°C 
H 
abt 4.5108 2x10" 2x10" 1x10" 
a 2x16°° 3x xis x10 
a *” 
Het 
abf 40 1x10? 3x10" 1x10" 
r ~1 6x 10-6 2x10-* 3x10-" 
abl.» 40 60 60 30 
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Taking, as an example, LZ =2X105L(©) =8 x 10* 
erg/sec., R*=R=3R(O)=2X10" cm.” (the 
data for W-R component of H.D. 193576) and 
V=2 X 108 cm/sec., we obtain for the density in 
the accelerative region : 


p*21.5X10-” g/cm*. 


This value is in good agreement with Aller’s” 
estimate of densities in W-R atmospheres. 

Turning to the particular type of the radiation- 
pressure mechanism which will be of most im- 
portance at these low densities and high tem- 
peratures (25,000 to 200,000 °C) involved, we 
must distinguish between three kinds of absorp- 
tion : (1) absorption due to bound-free transitions 
(b-f) in highly ionized atoms; (2) absorption due 
to the so-called free-free transitions (f-f) of 
electrons in the neighborhood of charged nuclei; 
and finally (3) absorption due to Compton scat- 
tering by free electrons (s).”* 

The absorption cross section for b-f transitions 
in hydrogen-like ion of atomic number Z is: 


642‘m,e" 2 
GP tmeerinnnen, (9) 
3v3h%c* 


where » is the frequency of incident radiation. 
Multiplying this expression by Planck’s distri- 
bution function J,, integrating over all frequen- 
cies 2 v,° (v,9=(2rm.e*/h*)z*), and dividing by 
the mass Amy of the atom in question, we obtain 
for the radial acceleration of the absorbing atom : 


ea) a,»-! } a 
a,f=] r -—dv 
v20 Amy Cc 





2562°m,.e'Z' 
& 3v3h*c'mpA 
*kT-In (1—exp (—hy.°/kT))! 


Zz} 


=7.3-1022-— 
A 





-kT In (1—exp (—/Ay,°/kT))—. (10) 


1Qn the basis of the absorption coefficient in W-R 
atmospheres one arrives at the conclusion that the thick- 
ness of the accelerative layer must be comparable with 
the stellar radius. 

21. H. Aller, Astrophys. J. 97, 135 (1943). 

23 Compare, for example, S. Chandrasekhar, Introduction 
into the Study of Stellar Structure, Chap. VII. 


On the other hand the fraction 7 of all atoms 
of the given element which retain their K ele. 
tron is given by the thermal ionization formula: 


2 





(2xm,) 


ht 





3 
(kT)i exp (—hv,°/kT) 
2 
141 10°(kT)! exp (—hv,°/kT) 


1+2 


e 





(11) 


where JN, is the electron density. Taking V,2108 
electrons/cm* (which corresponds to p*210-® 
g/cm*) we obtain the following values of a, 
and » for H and Het, respectively (see Table J), 
The third and sixth lines of the table represent 
the actual acceleration (in cm/sec.”) expected in 
the stellar atmosphere built entirely from the 
given gas. We notice that the effect on Het jg 
somewhat larger than that on H, but the calcu. 
lations show that beyond helium radiative 
acceleration again decreases rather rapidly. It 
may be mentioned here that in his work on radi- 
ation pressure in hot stars, Gerasimovit?* comes 
to much higher values of expected accelerations, 
This is, however, due to the fact that he assumes 
the densities in accelerative regions to be much 
higher (by a factor of about 10,000) than is 
assumed in the present article; such high den- 
sities would be, in fact, entirely inconsistent with 
the momentum conservation law as expressed by 
our Eq. (8). 

Turning now to f-f transition, we have the 
following formula for the effective cross section; 


4v2-nte®N.Z* 1 
~ 3V3mahc(kT)* v® 


f—f 


(12) 





The total acceleration of stellar matter comes 
out to be: 
me 8v2-rle®N, - c.Z" 
Atotal = = : kT) ——— 
3v3mShc'my ~ A 





OF A) 


; cL" 
=1.8-105(kT)§ + — (13) 


where ¢, is the relative abundance of element Z. 
For an atmosphere consisting of 65 percent of 
Russell mixture of heavier elements the sum in 


* B. P. Gerasimovit, M. N. R. A. S. 94, 737 (1934). 
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(13) is equal to 4. Taking again N.=10", we 
find from (13) that for the temperatures from 
25,000 °C to 200,000 °C, total acceleration due to 
free-free transitions will be only between 0.1 and 
0.6 cm/sec.’. . 

In case of scattering by free electrons the 
cross section is given by the well-known formula: 


87e! 
=6.6X10-** cm? (14) 





a.= 
3m.c* 


and the corresponding acceleration : 
a, ol* T* 


-——=3.8X10-'*—, (15) 
MeN 2c Me 





G,= 


u, being the mass per free electron. Because of 
the strong ionization we can assume yp, to be not 
much larger than 3 or 4. The formula (15) gives 
then that with 7 varying from 25,000 to 200,000 
°C, acceleration a, will increase from 40 to 
160,000 cm/sec.’. 

Thus we come to the conclusion that under the 
conditions existing in the accelerative layers of 
massive contracting stars radiation pressure due to 
scattering by free electrons must play the most 
important role in the process of ejection of stellar 
atmospheres. 

For the ratio of radiation pressure to the force 
of gravity on the stellar surface we have: 


is ( a,L ) eS) 
Qgrav . Mem H8r Rc ’ R? 
(L/L(O)) 


=4x10-*-__—_——._ (16) 
(M/M(©)) 








Thus the condition that radiation pressure 
overbalances gravity may be written in the form: 


(L/L(O))/(M/M(©))>2.5X10°. (17) 


Applying this condition to the normal stars of 
the main sequence, we find that even for most 
luminous stars [L=2X10°L(©); M=40M(O)] 
the left side of (17) is only 5 X 10*. Since, however, 
luminosities of contracting stars reach much 
higher values for the same mass, conditions for 
the ejection become more favorable. Using maxi- 
mum luminosities calculated above for stellar 


masses 10M(©), 20M(©), and 40M(©), we find 





EVOLUTION OF CONTRACTING STARS 


29 





for the ratio of radiation pressure to the force of 
gravity 0.1, 0.2, 0.3, respectively. These values 
are still smaller than unity, and, since the mass 
ejection by W-R stars is a well-established fact, 
the conclusion seems to be necessary that in our 
estimate of the ratio araq/@gray We have either 
underestimated the value of radiation pressure or 
overestimated the acting force of gravity. The 
latter possibility seems quite reasonable since the 
gravity on the stellar surface can be largely 
reduced by the centrifugal force resulting from 
the axial rotation of the star; and it is known, in 
fact, that most stars belonging to the upper part 
of the main sequence rotate very rapidly with 
equatorial velocities easily reaching the values 
as high as 200 km/sec. Such high rotation 
velocities, increased still further by the steady 
contraction of stellar body, will help the radiation 
pressure in the process of atmospheric ejection, 
and will make this process possible already for 
smaller values of star masses. Further observa- 
tional studies in this direction are evidently very 
desirable. 

It must be noticed in the conclusion of this 
section that all the treatment of the dynamical 
state of stellar atmosphere under the action of 
escaping radiation has been carried out only in 
a very schematic form. In a more exact treatment 
it would be necessary to take into account the 
effect of re-emission of the absorbed energy, 
leading to the backflow of radiation through the 
stellar atmosphere, and exercising ‘breaking 
action’’ on the outgoing gas masses. Some inter- 
esting contributions to this problem can be found 
in the works of Gerasimovit. 


5. LATE STAGES OF CONTRACTIVE EVOLUTION 


Since in the later stages of the evolution the 
energy liberation in the star is due exclusively to 
gravitational forces, we can work out an evolu- 
tionary time scale on the basis of the self-evident 
equation :* 


3 
1.3-MG(—-—) =1-a4 (18) 
R_ Ro 


The resulting time changes of stellar radius 
for star masses 10M(©), 20M(©), and 40M(©) 
2°The numerical coefficient 1.3 in the expression of 


potential energy corresponds to the mass distribution in 
the contractive model. 
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are shown in Fig. 6a. We see that the contraction 
of the star becomes quite noticeable a few 
thousand years after the luminosity has reached 
its maximum value. Since the total luminosity 
does not change any more, the surface tempera- 
ture of the star must increase inversely with the 
square root of the radius (Fig. 6b). With the 
increasing surface temperature, the maximum of 
the energy curve will move further and further 
into the ultraviolet, and the visible part of the 
spectrum will be more and more weakened. The 
visual luminosity of the star will be proportional 
to the surface temperature and to the square of 
the stellar radius; since 7yu.s~R-' we have 
Lvis~Rt!. These luminosity changes, expressed 
in terms of stellar magnitudes, are shown graphi- 
cally in Fig. 6c, and we see that a contracting 
star must become visually fainter by two or three 
magnitudes in the period of a few thousand years. 

We can now ask ourselves what becomes of the 
gases which are being continuously ejected from 
the surface of a contracting star. Let us suppose 
first that they continue to move outwards with 
the velocity of several thousand km/sec. observed 
near the surface of the star. According to the 
continuity law 


pr? V =const., (19) 


the density in the expanding envelope will vary 
inversely with the square of the distance from 
the star. Since the density of gases in the ac- 
celerative layer near the surface of the star is 
about 10-” g/cm’, the density of the envelope 
will drop to 10-*° g/cm? at the distance of 104 
stellar radii (or 200 astr. units), and to 10-” 
g/cm® at the distance of 10° radii (2000 astr. 
units).2* The later value represents the average 
size of the celestial objects known as planetary 
nebulae, and consisting of a very hot central star 
(planetary nucleus) surrounded by a luminous 
gaseous envelope. The study of planetary en- 
velopes carried out by Menzel and Aller?’ leads 
to the densities about one hundred times lower 
than the above-given estimate of the expected 
densities for the W-R ejection. Since the lumi- 
nosity of the gaseous envelopes is due entirely 


*¢ Because of the high velocities of ejected gases, it 
would take only 5 years before the envelope will expand 
to that size. 

(94h). H. Menzel and L. H. Aller, Astrophys. J. 93, 195 
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Fic. 6. Time changes of the properties of a contracting 
star which has reached the maximum of its total luminosity, 
() (top) Radius, (b) (middle) surface temperature, (¢) 

bottom) visual magnitude. 


to the re-emission of light coming from the 
central star and must be roughly proportional to 
the amount of matter forming the envelope, we 
should expect W-R envelopes to be relatively 
about one hundred times fainter than those of 
the planetary nebulae. The relative intensity can 
be still further reduced because of the fact that 
W-R stars seem to have somewhat lower surface 
temperatures than the nuclei of planetary 
nebulae, and correspondingly possess less inten- 
sity in the active ultraviolet part of their spec- 
trum. It is difficult to give the exact estimate of 
the expected luminosity of W-R_ envelopes 
without more detailed calculations, but it seems 
that one should be able to observe them at least 
as a faint glow around the star. 

On the other hand, observational studies of the 
known W-R stars do not give even a slight indi- 
cation of the existence of such envelopes. If this 
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conclusion will be confirmed by more careful 
observations and more detailed calculations of 
the expected intensity of W-R envelopes, the 
conclusion will be inevitable that the ejected 
gases do not continue to move with the ofiginal 
velocity, but are probably slowed down by the 
gravitational forces of the star. Such an assump- 
tion seems quite reasonable since the escape 
velocities for W-R stars, as calculated from their 
masses and radii, are of the same order of mag- 
nitude as the observed original ejection-velocities. 
Thus, for example, for M=10M(©) and 
R=3R(O©) we have Vescape= 1200 km/sec. 

If we assume that the gases ejected from the 
typical W-R stars are decelerated by gravity to 
such an extent that their original velocity is 
reduced by a factor of a hundred, the equation 
of continuity (18) will lead to the densities of 
W-R envelopes comparable with those of plane- 
tary nebulae for equal distances. Furthermore, 
this would reduce the expansion velocity to only 
20 km/sec., thus making it of the same order of 
magnitude as the actually observed expansion 
velocities of planetary envelopes. On the other 
hand, the time period necessary to reach the size 
of a planetary nebula will be lengthened to 
several hundred or a thousand years. With the 
further progress of the contraction, visual 
luminosity of the central star will decrease, and 
the gaseous envelope will become more lu- 
minous because of the relative strengthening of 
the ultraviolet part of the spectrum. All in all, 
the aged W-R star will look rather similar to a 
typical planetary nebula, and the suggestion may 
be put forward that the planetary nebulae observed 
at present actually represent the results of the 
evolution of the W-R stars of the past. 


6. ULTIMATE FATE OF THE CONTRACTING STAR 


In considering the ultimate fate of a con- 
tracting star, we must first of all take into con- 
sideration the dissipation of mass due to the 
atmospheric ejection. Using the data given in the 
previous section, we find that the outgoing flow 
of gases (both in case of W-R stars and planetary 
nebulae) carries away about 10** g/year or 10-5 
sun mass per year. Thus far a star of 10M(©) 
substantial mass changes cannot be expected to 
take place in less than say a hundred thousand 
years. On the other hand, the central temperature 
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of the contractive star increases inversely with 
its radius, and as it can be seen from Fig. 6a, will 
increase by a factor of 5 in about 10‘ or 10° years 
after the star has reached the maximum of its 
luminosity. Assuming, as in: previous sections, 
that the maximum luminosity is reached for the 
central temperature of about 2X10*® °C, this 
would correspond to the interior of the star 
heated up to 1X10* °C. 

Now, it has been shown. by the present author 
and Schoenberg** that the temperature of this 
order represents a threshhold for the so-called 
urca-processes by which a large number of 
neutrinos can be produced in the stellar interior. 
Escaping through the body of the star these 
neutrinos carry away the heat energy of the 
central regions at a rate which can be as high as 
10° or even 10" erg/g sec. Thus the formation 
and the escape of neutrinos act as a kind of 
thermostatic arrangement preventing the central 
temperature of the star from rising much above 
10° °C. Under such conditions, “adiabatic” de- 
pendence of gas pressure on gas density in the 
stellar interior will be reduced to simple propor- 
tionality (instead of the usual relation p~p*’*), 
and the star will become radially unstable. This 
can be simply seen from the fact that the increase 
of gas pressure (Pgas~p) produced by a small 
contraction of stellar body will be unable to 
counteract the increase of interior pressure due to 
gravity (Pgrav~p*/*). Thus the contraction will 
continue, and go over into a major collapse. In 
Fig. 4 we have indicated by vertical bars the 
points of contractive-evolution tracks where such 
a collapse is expected to take place. 

The treatment of the dynamical problem of a 
collapsing star is, of course, immensely difficult, 
and we can give here only some preliminary qual- 
itative considerations concerning this process. 

A comparatively simple solution could be ob- 
tained under the assumption of the spherical 
symmetry of the collapse motion. In this case we 
would get a rapidly but uniformly shrinking 
model with the radiation streaming outwards 
along the radii in accordance with the ordinary 
opacity formula. Since the collapse is expected 
towards the end of the evolution, when the 
opacity of stellar matter had reached its mini- 


28 G. Gamow and M. Schoenberg, Phys. Rev. 59, 539 
(1941); G. Gamow, Phys. Rev. 59, 617 (1941). 
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mum value (scattering opacity), and since the 
temperature in the center cannot rise above 
10° °C, the flux of the radiant energy through 
the body of the spherically-symmetrically col- 
lapsing star cannot be substantially larger than 
that existing before the collapse started. Thus 
almost the entire gravitational energy liberated 
in the collapse will have to be removed by 
neutrinos generated in the urca-region, and no 
sudden flare up characterizing stellar explosion 
can be expected. The star would simply shrink 
very rapidly while keeping its original luminosity. 

It must be clear, however, that the spherical- 
symmetrical solution cannot possibly correspond 
to the reality, in the case of a collapse taking 
place within only a few hours. It is, in fact, much 
more likely that the lack of supporting pressure 
in the stellar interior will cause a very irregular 
motion of gas masses, some of them precipitating 
towards the center and the other being ‘‘squeezed 
out”’ and thrown away in the form of vigorous 
prominences. The rotation of the star will exer- 
cise some kind of directing influence on this 
phenomenon since the presence of centrifugal 
force will hinder the inward motions of matter in 
equatorial regions. We may expect that the polar 
regions of a collapsing star will, so to speak, cave 
in, ejecting masses of central and equatorial 
regions. Radiation energy, which is saturating 
the material of stellar interior, will be thus bodily 
carried out and set free in the outer regions of 
the star. A rough estimate of the amount of 
radiative energy which can be liberated in such 
a process and of the velocities of ejected gas 
masses seems to be in satisfactory agreement 
with the observational data on supernovae ex- 


plosions. In particular it is evident that magg 
ejection by supernovae must be of a purely 
mechanical character since the total mechanical] 
momentum of the ejected envelope considerably 
exceeds the total momentum of the emitted radj. 
ation. In fact, according to Baade** and Minkoy. 
ski*® the velocity of gases ejected by the Chinese 
supernova of the year A.D. 1054 is about 1009 
km/sec., whereas the total mass of the gas shelf 
may be as large as 15 sun masses. These numbers 
lead to the value of 3X10 erg-sec./cm for the 
total mechanical momentum of the shell. On the 
other hand the total amount of light energy 
emitted by a typical supernova is of the order of 
magnitude of 10** erg corresponding to the 
momentum of only 10** erg-sec./cm. One can 
argue, of course, that the above-given number 
does not take into account the invisible part of 
supernova spectrum, but it is extremely unlikely 
that the actual luminosity of a supernova could 
exceed by a factor of ten thousand its visible 
luminosity. Thus we can conclude with a large 
degree of certainty that in contrast to the mass 
ejection by a W-R star, nuclei of planetary 
nebulae, and possibly in contrast to the ordinary 
novae, the explosion of supernovae is of more 
mechanical character and is directly connected 
with the dynamical peculiarities of a collapse 
process. 

In conclusion the author considers it to be his 
pleasant duty to express his gratitude to his 
astronomical friends Dr. J. A. Hynek and Dr, 
Ph. C. Keenan for many helpful discussions on 
the subject. 


29 \W. Baade, Astrophys. J. 96, 188 (1942). 
%°R. Minkowski, Astrophys. J. 96, 199 (1942). 
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An Independent Determination of Fixed Points on the High Voltage Scale 


A. O. HANson AND D. L. BENEDICT 
Department of Physics, University of Wisconsin, Madison, Wisconsin 
(Received September 14, 1943) 


Since present methods allow the comparison of high voltages with an accuracy of 0.1 percent, 
it was felt that a more accurate determination of certain fixed points would be desirable. 
This work is an attempt to establish more precisely the voltages of several well-known nuclear 
reactions which have been used as convenient voltage calibration points. The absolute values 
of the proton energies obtained for the various calibration points are: 1.883 Mev for the 
Li(p, ») threshold, 2.058 Mev for the Be(p, ”) threshold, 0.877 Mev for the first strong F(p, y) 
resonance, and 0.4465 for the Li(p, 7) resonance. The absolute values are thought to be accurate 
to within 0.3 percent and the relative values to about 0.1 percent. 





HE high voltage scale used up to the present 

time is based on a determination of the 
voltage at the Li(p, y) resonance by Hafstad 
et al., using a calibrated resistance voltmeter.! 
This measurement, which was also checked by 
Parkinson ef al.,? gave a value of 0.440 Mev for 
this resonance. The consistency of their measure- 
ments indicated a probable error of around two 
percent but additional checks on the high voltage 
scale obtained from the scattering of protons in 
argon® and the C(p, m) threshold‘ seemed to in- 
dicate that the scale was correct to within one 
percent. The values found in this work, however, 
are approximately 1.5 percent higher than 
previously accepted values. 


EXPERIMENTAL METHOD 


The method used in this work consisted of 
measuring the energy of the proton beam at the 
Li(p,) threshold with a calibrated electrostatic 
analyzer. The calibration of the analyzer was 
made by using electron beams of accurately 
measured energies and was checked by calcu- 
lating the deflection constant of the analyzer 
from its geometry. The proton energies corre- 
sponding to the Be(p, 2) threshold and to the 
Li(p, y) and F(p, y) resonances were found by 


'L. R. Hafstad, N. P. Heydenburg, and M. A. Tuve, 
Phys. Rev. 50, 504 (1936). 

*D. B. Parkinson, R. G. Herb, E. J. Bernet, and DBs 
McKibben, Phys. Rev. 53, 642 (1938). 

*N. P. Heydenburg, L. R. Hafstad, and M. A. Tuve, 
Phys. Rev. 56, 1078 (1939). 

“R. O. Haxby, W. E. Shoupp, W. E. Stephens, and 
W. H. Wells, Phys. Rev. 58, 1035 (1940). 


comparing these energies with that at the 
Li(p, m) threshold. 


ANALYZER DESIGN AND CALIBRATION 


An electrostatic analyzer with curved plates 
and a single fixed detector had been used for 
some time in connection with the concentric 
electrode electrostatic generator for measuring 
the energy of the ion beam. This analyzer and 
the auxiliary apparatus are described elsewhere.‘ 
In the present work, since the analyzer was to 
be calibrated with electrons of moderate energy, 


‘it seemed desirable to be able to reverse the de- 
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flector voltage as a check on the effect of stray 
magnetic fields and of possible surface charges. 
The analyzer was, therefore, made to be of the 
parallel plate type with three fixed detectors as 
shown in Fig. 1. One detector was placed on the 
axis of the analyzer in order to establish the 
position of the undeflected beam, and the other 
two were equally spaced above and below it. 
The position of the lower detector was made 
slightly adjustable by means of an eccentric cam. 
This position was adjusted during preliminary 
runs with the electron beam so that reversal of 
the deflector voltage caused the beam to shift 
from the center of the upper detector to the 
center of the lower detector. A leveling screw on 
the detector end of the analyzer provided the 


’ A. O. Hanson, ‘‘Voltage-measuring and -control equip- 
ment for electrostatic generators,” recently submitted to 
The Review of Scientific Instruments. This paper will be 
referred to elsewhere in the text as (I). 
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Fic. 1. Electrostatic deflector. 


adjustment for centering the undeflected beam 
on the center detector. 

For the calibration experiment the analyzer 
was mounted with its axis parallel to the earth’s 
magnetic field. The exterior of the tube contain- 
ing the deflector was wound with wire and a 
current was maintained in the solenoid thus 
formed so as to reduce the magnetic field in the 
interior to a low value. 

The electron gun used to produce the beam of 
electrons is also shown in Fig. 1. An additional 


lens for focusing the beam was provided. This . 


lens was shorted out, however, allowing the full 
voltage to be applied across the first gap since 
concentrating the beam by use of the additional 
lens made the density of the beam less uniform. 

The gun was mounted about 30 cm ahead of 
the slit which defined the flat beam through the 
analyzer. The slit width was about 0.2 mm and 
the beam as seen on the glass behind the de- 
tectors appeared to be a uniform line about 0.5 
mm wide. At first it was found that the beam 
hit the mica washers, which insulated the de- 
tector plates, causing them to charge up and 
distort the beam. This difficulty was overcome 
by inserting a metal shield which allowed the 
beam to strike only the central semicircular 
parts of the detector plates. 

The schematic diagram of the electrical circuit 
for the electron calibration experiment is shown 
in Fig. 2. The voltage used to accelerate the elec- 
trons was obtained from the negative side of the 
high voltage deflector supply described in (I). 


A bank of B batteries was used to furnish the 
deflector voltage. The accelerating voltage and 
the deflector voltage were measured by means 
of good voltage dividers used with the same 
potentiometer. 

The data shown in Table I give the numbers 
upon which the calibration is based. The pro- 
cedure in taking these data was as follows: (1) To 


align the analyzer so that the beam was centered 
on the central detector with no voltage across 
the plates. The correct centering of the beam on 
the detectors was indicated by a zero reading on 


the galvanometer. (2) To set the deflector 
voltage at a given value and adjust the acceler- 
ating voltage until the beam was centered on the 
upper detector and to measure this voltage with 
the potentiometer. (3) To reverse the deflector 
voltage, readjust and measure the accelerating 
voltage when the beam was centered on the lower 
detector. (4) To check the deflector voltage. All 
the potentiometer readings were taken to four 
significant figures and were multiplied by the 
proper constants to get the values in the table 
which are recorded in volts. The values of V,/ Vz 
represent the average accelerating voltage divided 
by the average deflector voltage for each run. 

As shown in Table I the deflector constant 
varies with the energy of the electrons because 
of the relativistic effect. In order to correct for 
this effect it is necessary to examine the equation 
for the deflection of a charged particle by a 
condenser. 

The deflection at the far edge of an ideal con- 
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denser is given by 


d=6&et?/2m= Vuel2/2Smv, (1) 


where & is the electric intensity, Vz the voltage 
across the plates, ¢ the time for the particle to 
travel the length of the plates, / and S the length 
and separation of the plates, and e, m, and v the 
charge, mass, and velocity of the particle. If the 
detector is placed at a distance L from the center 
of the deflector plates, the deflection D at the 


distance L is 


D=Ld/ (1/2) = VaeLl/ Smo’. (2) 


TABLE I. Data on deflection measurements. 

















Va V. 

De- Accelerating Deflecting Ve a ‘in Veo ) 
tector voltage voltage Va Va 2(511 +Vay,) 
Upper 11,749 218.20 
Upper 11,747 
Upper 11,744 
— seer 53.838 53.236 
Lower 11,749 
Upper 11,744 
Lower 11,746 218.14 
Upper 19,520 360.20 
oan ie 54.315 53.217 
Lower 19,534 360.04 
Upper 8,263 153.96 
— oa 53.669 53.223 
Lower 8,263 153.96 
" Neutralizing magnetic field raised 50% 

pper 8,263 153.96 
Lower 8,268 153,96 93-682 53.238 
- 7 neutralizing magnetic field 

pper 8,2 154.00 
Lower 8,277. — «154.00 53-747 53.301 











Fic. 2. Wiring diagram for electron calibration experiment. 








If we use the non-relativistic relation mv’ = 2 Vie, 
where V, is the accelerating voltage applied to 
the particle, we find that the deflector constant 
for a given deflection is simply 


Ve/Va=Ll/2SD. (3) 


Where the velocity of the particle approaches 
the same order of magnitude as that of light we 
must use the transverse mass and the actual 
velocity of the particle in the kinetic energy 
relation. If E is the total energy of the particle 
and E, the rest energy, the kinetic energy is 


V,.e=E—Ep,, 


where 


moc? 
E=————-=me and Eo=m ec? 
(1—v*/c)! 


and m and mp represent the transverse mass and 
the rest mass, respectively. From these relations 


we find 
v me’ EF-E? 
ce me FE? 
F?-E,? 
——— = 2V.4( 1- 
E 


——_.), (4) 
2(Eot+ Vie) 
Equation (3) then becomes 


Ve (: Vie ) 

Va 2(Eot+ V.e) 
The left-hand side of this equation was deter- 
mined experimentally and represents the de- 


and 


mv* 


V.e 


Li 
28D 
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flection constant of the analyzer for particles 
having non-relativistic velocities. The data in 
Table I show that the value obtained for this 
constant is independent of the energy of the 
electrons. Magnetic effects are almost negligible 
and effects due to surface charges are also unim- 
portant since the results agree consistently when 
the deflector voltage is reversed. The value of the 
constant is seen to be 53.22+0.10. The estimated 
uncertainty of 0.2 percent is based on the fact 
that the resistors used in the voltage dividers 
were determined with a precision of 0.1 percent. 


CALCULATION OF THE DEFLECTOR CONSTANT 





After the analyzer was used to determine the 
proton energy at the Li(p, m) threshold, as de- 
scribed later in this paper, it was dismantled and 
the important dimensions were accurately meas- 
ured. These dimensions were: The length of the 
deflector plates, /=6.17 cm. The distance from 
the center of the deflector plates to the center 
detector, L=47.59 cm. The average separation 
of the plates, S=0.994 cm. The variation in the 
separation of these plates was 0.005 cm. The 
distance between the center detector and either 
of the other detectors, D=3.132 cm. 

The deflector constant can be calculated ap- 
proximately from Eq. (3) but in order to get any 
accuracy it is necessary to consider the effect of 
the fields at the ends of the deflector plates. The 
conditions at the end near the entrance slit are 
sufficiently like those considered by Herzog to 
allow reading the correction from the graph 
(Fig. 3) given in his paper.* The conditions at 
this end correspond to a thick aperture having 
a width 0.80 S (S being the plate separation) 
placed at a distance of 0.87 S from the edge of 
the plates. This correction (C,) amounts to 0.31 
cm. The conditions at the other end of the 
plates correspond to a very large exit slit and the 
correction cannot be read directly from Herzog’s 
graph. An independent estimate of this correction 
was made by determining the equipotential lines 
of a model of the analyzer placed in an electro- 
lytic tray and calculating the effect of the field 
at this exit end of the plates. The value for the 
end correction at this end of the plates (C2) 
obtained in this way was 0.50 cm. This is in 


*R. Herzog, Physik. Zeits. 41, 18 (1940). 
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Fic. 3. Neutrons from the Li(p, m) reaction. Runs I, II, 
and III taken with the beam on the upper detector, 
Runs IV and V taken with the beam on the lower detector, 







agreement with the value extrapolated from the 
graph given by Herzog. 

The length of the plates used in calculating the 
deflection constant is therefore 1’=/+(C,+(, 
= 6.98 cm. The distance L is also changed since 
the corrections at the two ends are not the same 
and becomes approximately 


L’=L+32(Ci— C2) = 47.44 cm. 


The deflection constant as calculated from these 
values is 


V,/Va=L'l'/2SD =53.18-+0.30. 












The larger uncertainty given in this case is due 
to the larger probable error in determining the 
end corrections. The value for the deflector 
constant obtained in this way is, however, in 
good agreement with the experimental value. 









THE DETERMINATION OF THE LITHIUM: 
PROTON-NEUTRON THRESHOLD 







For the determination of the threshold voltage 
for the Li(p,m) reaction, the analyzer was 
mounted on one of the outlets of the concentric 
electrode electrostatic generator. Metallic li 
thium was evaporated upon a slitted tantalum 
target and was placed in the tube in front of 
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Fic. 4. Gamma-ray yield from a thick lithium target using 
the diatomic beam. 


the defining slit to the analyzer as shown in 
Fig. 1. A }-inch circular aperture placed ahead 
of the target tube defined the beam before it 
struck the lithium target. A galvanometer was 
used to measure the total proton current entering 
the target tube but approximately one-half of 
this current passed through the slit in the lithium 
target and struck the smaller slit defining the 
beam passing through the analyzer. With this 
arrangement the part of the beam which passed 
through the analyzer should have the same dis- 
tribution in energy as that which hit the lithium 
target. 

A simple ionization chamber filled with boron 
trifluoride placed as shown in Fig. 1 was used to 
detect the neutrons. The data on the neutron 
yields as a function of proton energy were ob- 
tained by setting the deflector voltage at a 
definite value and adjusting the voltage of the 
generator so as to keep the beam centered on the 
upper (or lower) detector. A regulator operating 
from the detector plates served to keep the beam 
centered automatically. This regulator as well 
as the stabilized deflector voltage supply is 
described in (I). It was estimated that the fluc- 
tuations in the energy of the beam were at the 
most 0.3 percent and that the average energy 
of the beam was constant to about 0.1 percent. 
The deflector voltage remained constant to 
better than 0.02 percent. The neutron intensity 
as a function of proton energy is shown in Fig. 3. 
The runs were taken with different beam currents 
and various amounts of paraffin surrounding the 
target tube and the ionization chamber. The 
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first three runs were taken with the beam centered 
on the upper detector and the remaining two 
with the deflector voltage reversed. The last run 
(No. 5) was taken with all paraffin removed. The 
ordinates for the various runs shown in Fig. 3 
were adjusted and displaced for convenience in 
plotting. In all but the last run the intensity at 
the end of the run was as high as the detector 
system could record accurately (about 300 counts 
per second with a total proton current of one-half 
microampere). There was no background count 
below threshold except that due to voltage 
fluctuations. 

It is seen that the extrapolated threshold lies 
at a potentiometer reading of about 1.265. Most 
of the counts recorded at this potentiometer 
setting were observed to occur during voltage 
fluctuations. The counts at the potentiometer 
setting of 1.267, however, came in almost con- 
tinuously indicating that the energy of the 
protons was definitely above the threshold 
energy. The extrapolated threshold is probably 
the most reliable value to take for the threshold 
since it should be independent of small voltage 
fluctuations. The potentiometer reading of 1.265 
corresponds to a deflector voltage of 35.35 kilo- 
volts and if the deflector constant is taken as 
53.22 the energy of the protons at the threshold 
is 1.881 Mev. If we consider the small relativistic 
effect for protons at this energy we find that the 
above value should be raised by about 0.1 percent 
giving a final value of 1.883 Mev. 


COMPARISON WITH OTHER STANDARDS 


Although the relative values of other calibra- 
tion points were considered to be quite reliable 
it was thought that it might be of interest to 
compare them directly by use of the analyzer 
which was at that time permanently on the elec- 
trostatic generator. (This was the first analyzer 
with the 40-cm collimator described in I.) The 
targets used in this work were thick targets of 
lithium and beryllium metal and potassium 
fluoride. These targets were mounted on a single 
target holder in such a way that any of the 
targets could be placed in the beam without 
breaking the vacuum. The threshold energies for 
the Li(p, 2) and Be(p, m) reactions were deter- 
mined from the intercepts of curves similar to 
those shown in Fig. 3. The gamma-ray yields 
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relative values of the calibration points are ip 
good agreement with the previous values. 

It is disturbing to find that the value obtained 
for the Li(p, ) threshold is as much as 1.5 per. 
cent higher than that based on the 0.440-Mey 
Li(p, y) calibration point. It is felt, however, 
that the values obtained in this work should be 
quite reliable since the sensitivity in both the 
calibration experiment and in the threshold work 
was such that the results were consistent to 
better than 0.1 percent. The wire wound re. 
sistors used to measure the voltage of the 40. 







TABLE IT. Resonance energies. 
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Fic. 5. Gamma-ray yield from a thick target of 
potassium fluoride. 







from -a thick target of lithium metal obtained 
with the diatomic beam and from a thick target 
of potassium fluoride with the proton and the 
diatomic beams are shown in Figs. 4 and 5. These 
yields were measured by means of a hydrogen 
filled Geiger-Mueller counter used with a Neher- 
Harper quenching circuit feeding a scale of 64 
recording circuit. The resonance energies were 
taken as that corresponding to the position of the 
maximum slope in the yield curves. These ener- 
gies and the previously accepted values are listed 
in Table II. The values obtained by using the 
diatomic beam are listed at one-half the energy 
indicated by the electrostatic analyzer. The 
values obtained for the fluorine resonance by 
using the proton and the diatomic beams indicate 
that the analyzer gives a linear measurement of 
energy to about 0.1 percent. It is seen that the 














0.440 Mev present in- previous 
Reaction Li resonance work crease work 
Li(p, m) 1.856 Mev 1.883 Mev 15 a 
Be(p, m) 2.028 2.058 1.5 a 
F(p, 7) 0.862 0.877 a 
0.867 1.2 c 
Li(p, vy) 0.440 0.4465 1.5 c 





















@ See reference 4 of text. 

+E. J. Bernet, R. G. Herb, and D. B. Parkinson, Phys. Rev, 54 
398 (1938). 

¢ See reference 1 of text. 







kilovolt supply were rechecked shortly after this 
work and were found to be accurate to within 
0.1 percent. The standard cell used with the 
potentiometer was checked before and after this 
work by the University of Wisconsin Standards 
Laboratory and was found to be reliable. It is 
still possible that there are some systematic 
errors in work of this nature, and it would be 
desirable to check at least one of these points 
by using a somewhat different method. 

We wish to express our appreciation to Dr. 
R. G. Herb under whose guidance this work was 
done and to the Wisconsin Alumni Research 
Foundation for financial support. 
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Interpretation of the Spinning Electron with Bipolar Coordinates* 


Mary HEw ett PAyNeft 
Brown University, Providence, Rhode Island 


(Received October 16, 1943) 


The bipolar coordinates of a point may be loosely 
defined as the distances of the point from each of two 
reference points. Consistent definitions are developed for 
the bipolar coordinates of a point on a sphere with the 
reference points at opposite ends of a diameter, and it is 
shown that the bipolar coordinates so defined transform 
like Pauli spin functions under rotations. The electron is 
visualized as a point mass with charge e moving on the 
sphere with its motion expressed in terms of its bipolar 
coordinates. It is assumed observable only at two dia- 


metrically opposite points of its orbit. This picture allows 
an exact physical interpretation of the spin variable and 
the Pauli spin operators. The magnetic moment associated 
with the motion is eh/2mc; the radius of the electron is 
taken as the radius of the sphere which is A/2mc; the 
wave-length associated with the spinning motion is h/mce. 
Two related functions must be introduced to take care 
of reflections and these together with the two original 
functions transform like Dirac spin functions under 
rotations, reflections, and Lorentz transformations. 





INTRODUCTION 


N this paper an attempt will be made to set 
up a geometrical picture for the Pauli and 
Dirac spin theories.' The proof of the equivalence 
of this picture with the analytical theories will 
be based on the transformation properties of the 
Pauli and Dirac spin functions. It will be con- 
venient to change the usual notation slightly, and 
we shall take the Eulerian angles ¢, 0, and y as 
in Fig. 1. g is a rotation about the z axis; 0, a 
rotation about the new y axis after the ¢ rotation 
has been performed; and y, a rotation about the 
final z axis. This differs from the usual conven- 
tions where @ is taken as a rotation about the x 
axis. The other important change is in the time 
variable. One usually takes the fourth coordinate 
as ict; we shall take it as ct. 

With the notation described in the previous 
paragraph, one finds that the Pauli spin func- 
tions a;(a) and a2(a@) transform under a rotation 
through ¢, 6, and y in that order by a unitary 


matrix U, given by 
—sin }6e(e-») 
. (A) 


e 36e-4 i(e+y) 
cos 3 be} i(e+y) 


sin }6e!'-#) 


a is the spin variable and is allowed to assume 


* Extract from a thesis submitted in partial fulfillment 
of the requirements for the degree of Doctor of Philosophy 
in physics at Brown University, Providence, Rhode Island. 

P Present address: 816 Picher Avenue, Joplin, Missouri. 

1W. Pauli, Zeits. f. Physik 43, 601 (1927); P. A. M. 
Dirac, Proc. Roy. Soc. A117, 610 (1928), A118, 351 (1928) ; 
E. L. Hilland R. Landshoff, Rev. Mod. Phys. 10, 87 (1938); 
B. L. van der Waerden, Gruppentheoretische Methode 
(Julius Springer, Berlin, 1932); L. De Broglie, L’ Electron 
Magnetique (Hermann et C'*, Paris, 1934). 


two values a; and az corresponding to spin in the 
positive and the negative z directions, respec- 
tively. The spin operators corresponding to the 
components of spin angular momentum in the 
x, y, and z directions are, of course, the Pauli 
matrices multiplied by? 3h: 


n(° ’ n(° ong 

s;= 9 S,= 

1 0 iit Te 
1 


’ 0 
$2= 5h 0 es | 


The significance of the*spin¥variable and the 
actual form of the spin ‘operators are not given 
in the Pauli theory. » 

The Dirac theory makes use of four com- 
ponent wave functions. If we denote the four 
components by yi, ~2, ws, and ¥, and investigate 
the transformation properties, we find that 


(2) 





Fic. 1. g=rotation about Z»9; @=rotation about Yj; 


¥ =rotation about 22. 


2h =h/2s; h=Planck’s constant. 
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by the matrix 


y,-(¥% °) * 
o= (4 U,) ( 


Under the space inversion (reflection in the 
origin) and the time reversal they transform by 
the matrices 


0 6 0 6 
Ui=+/ ) and Ur=+( ). (4) 
6 0 —65 0 

X 








Fic. 2. 


respectively, and under a Lorentz transformation 
they transform by 
e“ 0 0 O 
U,= 0 e 0 O (5) 
0 oO e O 
0 0 O 
cosh p= 1/(1—6?)!, 
sinh p= —8/(1—8?)!, B=v/c, (6) 


and v=velocity of the moving system. 


e~ 
where 


PRELIMINARY REMARKS ON BIPOLAR 
COORDINATES 


For the purposes of this discussion we take the 
Cartesian axes as shown in Fig. 2. Note that this 
is a right-handed system. Consider the sphere 
of diameter a with its center at 0, the origin. 
Points on this sphere have coordinates x, y, and 
z, such that 


x?+-y?+2?=a?/4. (7) 


We define the bipolar coordinates of a point on 
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under rotations the four components transform 


PAYNE 





this sphere as the distances of the point from the 
points (0, 0, $a) and (0, 0, —4}a) and call them 
u and v, respectively. The two points from 
which u and v are measured are called the right. 
hand reference point and the left-hand reference 
point, respectively. On occasion we shall call 
them N and S, respectively. We shall confine 
ourselves to points on the sphere. This is suf. 
ficient for the purpose of developing a model of 
the electron, and the definitions of the bipolar 
coordinates become very difficult for points not 
on the sphere. 

Suppose that we have a point P on the sphere, 
up and vp will depend on ¢, @, and ¥ where g, §, 
and wy are the Eulerian angles defined in the 
introduction. @ is the angle made by OP and the 
positive z axis; g is the angle made by the plane 
of OP and the z axis with the xz plane; y is the 
angle of rotation about the line OP (see Fig. 3), 
At first sight one might think that @ and ¢ are 
sufficient to determine the position of a point on 
the sphere uniquely, and it is true that they are, 
However, if we wish to think of the sphere asa 
whole and to specify the position of all the points 
simultaneously, the variable y is also necessary, 
One might think of the y dependence as giving | 
the answer to the question: ‘‘What rotation 
would be necessary to get the points from a 
given initial position to their present position?” 

Now it is evident that with our present def- 
nition of the bipolar coordinates we have no 
means of distinguishing among the various 
points on a circle formed by the intersection of 
a plane perpendicular to the z axis with the 
sphere. The next three sections will be devoted 
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to developing a consistent set of conventions for 
assigning bipolar coordinates to each point on 
the sphere. To facilitate the discussion we shall 


introduce the variables : 
Xo =2. (8) 


xi=x+ty, x1=x—ty, 


THE 6 DEPENDENCE OF BIPOLAR COORDINATES: 


We take a point P on a circle of diameter a 
with center at the origin in the xz plane. From 
plane geometry the magnitudes of u and v are 


given by 


u=asin 30, v=acos 30. (9) 


We shall take (9) as the definition of bipolar 
coordinates in a plane. Note that we can now 
distinguish between points placed symmetrically 
with respect to the z axis (e.g., P and P’ in Fig. 
4). Note also that each point has two sets of 
coordinates—that is, the period of u and v is 4z. 
In other words, the bipolar coordinates of a 
point in a plane (as here defined) are double 
valued. We shall observe this characteristic in 
the g and ¥ dependences as well. In this case we 
have the relations 


X1=x)*=x_,)=uv/a, xo=(v?—u?)/2a. (10) 


The matrix for a rotation through @ about 


the y axis is 
rae - 26 —sin 4" 
sin 36 cos $6 
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(11) 


We shall assume 


x:=uv*/a, x*=u*v/a, 
xo= (|v|*—|u|*)/2a. 


This is consistent with the equations used in the 
previous section. For all points on the sphere we 
have the following equations 


(12) 


|x1|=|x:*|=4asin 6, |u| =asin 48, 
|v| =a cos 38. (13) 
Therefore 


|x1| = |x1*| =|u| |v] =| uo]. (14) 


and Eq. (12) above merely specifies a phase 


factor. 


* The development of bipolar coordinates in a plane was 
first given by W. T. Payne in an unpublished paper. 
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Now if the point P is such that the plane of 
OP and the z axis makes an angle ¢ with the xz 


plane, 
x1 = $a sin 6e**, x,*=}a sin be-‘*, 
xo = 4a cos 8. (15) 

We may take 
(16) 


u=asin }0e*, v=acos }6e%, 


To proceed further we note that the angle ¢ will 
appear positive from the right-hand reference 


X 





v 
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point and negative from the left-hand reference 
point. We take this as sufficient justification to 
assume that the ¢ dependence of u and v shall 
be opposite and set 


a=-— 8. (17) 


Making use of this relation, we obtain 


a—B=2a=¢: a=—B=}¢ (18) 


and 
(19) 


u=asin $6e%", v=acos }6e—#, 


Again we note that half-angles appear. The 
dependence on ¢ distinguishes among the points 
on a circle formed by the intersection of the 
sphere with a plane perpendicular to the z axis. 
No two points have the same coordinates, but 
again each point has two sets of coordinates— 
that is, the bipolar coordinates are double 
valued in ¢ as well as in @. 

The matrix for a rotation through ¢ about the 


z axis is 
ete 60 
én ( a (20) 
0 eil¥ 
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To obtain the y dependence we shall consider 
the bipolar coordinates 


u=asin }6e%", v=acos}Oe—**. (21) 
First we rotate through ¢ about the z axis: 
u’=asin }0, v'=acos 30 (22) 
and then through @ about the y axis: 
u"’=0, v’ =a. (23) 


Now we apply a rotation through —y about the 
z axis. This is the correct rotation since it is 
about the line joining the point with the origin. 
We had to make this line coincide with either 
the y axis or the z axis since these are the only 
two axes about which we know how to take rota- 
tions so far. The rotation through —y gives 


(24) 


u’"’=0, vw’ =ae—®, 


Now we get the point back to the position 
specified by @ and ¢ by rotating through —é@ 
about the y axis: 


u**=a sin }0e~-*¥, v'*=acos 46e—*%; (25) 
and finally through — ¢ about the z axis: 
u*=a sin $6e8*"-), v*=acos }6et), (26) 


Again we note that half-angles appear, so that 
the coordinates are double valued in y and have 
period 4x. No two points on the sphere may have 
the same coordinates, but the y dependence 
makes possible an infinite number of coordinates 
for each point. We may write the general form for 
the bipolar coordinates of a point as follows: 


u=a sin }6e*"-), yv=acos }6e e+), (27) 


GENERAL ROTATION ABOUT THE ORIGIN 


Let-us perform an arbitrary rotation defined 
by the Eulerian angles ¢, 6, and y in that order. 
What will be the matrix for the transformation 
of the bipolar coordinates? Remembering that @ 
is a rotation about the z axis, we first operate on 
u and v with the matrix ®. @ is a rotation about 
the y axis, so we operate next with the matrix 0. 
¥ is a rotation about the final z axis, and so 
finally we operate with the matrix V, which is 


ev 0 
v-( ; a) 


(28) 
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Therefore the matrix U for the complete set of 
rotations is 








U=VO (29) 


and direct substitution shows that this matrix 
is identical with U, in Eq. (1). 







SPIN FUNCTIONS AND THE SPIN VARIABLE 





We shall interpret the spin variable as position 
on the sphere, and the restriction of the spin 
variable to two values implies that only two 
points on the sphere are of importance in dis. 
cussing electron spin. We shall interpret this as 
meaning that the electron, which we may 
visualize as a point mass moving over the 
surface of the sphere, may be observed only when 
it occupies one or the other of the two points 
specified by the spin variable. Now, what is the 
geometrical relationship of these points? We 
recall that if the z component of spin is to be 
observed 














d2(a) =(), 
@2(a2) = 1. 


a;(a@) = 1, 
a;(a@2) =(, 





(30) 






In our bipolar coordinate scheme 





un =0, 
us =d. 


Uy =a, 
vs=0, 





(31) 





If we identify vy and uy with a;(a;) and a2(a), 
and vs and us with a;(a2) and a2(a2), respec 
tively, then a; and ae refer to the right- and 
left-hand reference points, respectively, and we 
may say more generally that a; and ae refer to 
points at opposite ends of the diameter parallel 
to the direction in which spin is being observed. 

In the Pauli spin theory a; and a2 are supposed 
to give information about the direction of spin 
and refer to spin in the positive and negatives 
directions, respectively. At the point J in the 
bipolar coordinate scheme uy is zero, and wy is 
equal to the diameter of the sphere. Now oy is 
the coordinate of N relative to S and may be 
regarded as parallel to the z axis. Similarly us 
may be thought of as antiparallel to the z axis. 
Thus we may take the bipolar coordinates a 
associated with the direction of spin. If the axis 
of the bipolar coordinate picture and the direc 
tion of spin do not coincide, then none of the 
coordinates of the points on the spin axis is zero, 
In this case we think of the v's as associated with 
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the parallel and the w’s with the antiparallel 
components of the spin. 

There is one further implication of the scheme 
we have developed which should be mentioned. 
The electron has been taken as a point charge up 
to the present. This does not mean, however, 
that an observation on the diameter of the elec- 
tron will give the value zero. For in the course of 
time the electron will occupy two points a dis- 
tance a apart, and these two points may them- 
selves rotate in space. This means that on the 
average the electron will appear smeared over 
the sphere and the observation would record a 
as the diameter of the electron. We shall compute 
ain terms of universal constants ina later section. 


THE PAULI SPIN OPERATORS 


To find the analog of the spin operators‘ in the 
bipolar coordinate scheme, we must investigate 
the spin variable somewhat more precisely. So 
far we have used the spin variable to reduce the 
configuration space of the electron to two points 
and to find the geometrical relationship of these 
two points. That is, we think of the electron as 
observable at only two points which are at 
opposite ends of a diameter of the sphere of 
diameter a. Now we have assumed throughout 
that only points on the sphere are of any sig- 
nificance, and if we wish to imagine a path for 
the electron between its observable points, we 
must think of it as being on the surface of the 
sphere. Thus we may think of the mechanism for 
the transfer of the electron from one point in its 
configuration space to the other as a rotation. 

We assume that the electron travels with a 
velocity whose component in any direction is the 
velocity of light c. Then the resultant velocity 
will also be c from the Einstein formula for addi- 
tion of velocities. The angular velocity will then 
be given by 


w=c/(a/2) =c/cAdt=1/At, (32) 


where Af is the time required to traverse the 
radius at velocity c. Now 8, ¢, and y differ by 
at points located at opposite ends of a diameter, 
so that we may think of the change in a, the 
angle in question, in going from one point to the 


‘ The idea of using difference operators was proposed by 
. T. Payne. 


other to be of magnitude xr. Then 


Aa/2At=2/2At=}xw=}nr(Ax/At), (33) 


where x is the angle implicit in the angular 
velocity. The term 2At is used in the denominator 
of (33) since the velocity component along the 
axis is c and the time required, therefore, to 
travel from one point to the other is 2At. 

Now, following the usual quantum mechanical 
procedure, we take the angular momentum 


operator to be 
thA/Ax =i(h/2)x(A/Aa). (34) 


But since Aa is just 7, the angular momentum 
operator reduces to 


3thd (35) 


or a finite rotation. 

There will be three independent rotations 
about the sphere each giving a path joining the 
two reference points which the electron might 
follow. One is a great circle in the xz plane 
(rotation through w about the y axis); another 
is a great circle in the yz plane (rotation through 


mw about the x axis); and the third must contain 
a rotation through wz about the z axis. It may 
be taken as a rotation through $2 about the y 
axis followed by a rotation through x about the 
z axis and finally a rotation through 42 about 
the y axis. 

Let us compute these rotations about the x, 
y, and z axes. The rotation about the x axis is 
given by a rotation through $m about the z axis 
followed by a rotation through about the new 
y axis and then by a rotation through —}x 
about the final g axis: 


QO --—1 
Rae=( } ). 
—2t 0 


where the subscript x means that the rotation is 
to be taken about the x axis, and the subscript 
m means that the angle of rotation is x. For the 
y axis we just rotate through x about the y axis: 


0 -!1 
nin? 
© 


For the z axis we carry out the sequence of 


(36) 


(37) 
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rotations given above: 
¥ iy a, —1 
R,-= zs ) 
1 0O = io 1 0 
* 2 (38) 
“No TF 


Now the operators sz, sy, sz of the Pauli theory 
should be given by 


Se= HAR. 2, Sy=HthR,,., S:=HAR,«, (39) 
as indeed they are. 


SOME DEDUCTIONS FROM BIPOLAR 
COORDINATE THEORY‘ 


We may compute the magnetic moment of the 
electron. If a particle of mass m and charge e 
moves in an orbit, then its magnetic moment po 
is e/2mc times its angular momentum s, if the 
moment is measured in e.m.u. This is proved by 
considering the torque on the system due to an 
external magnetic field. The magnetic forces will 
act on the electron only when it is moving per- 
pendicular to the field. If the electron is moving 
in a circle in the plane of the field, then at two 
diametrically opposite points in its orbit it will 
move perpendicular to the field and at points 
halfway between these two points it will move 
parallel to the field. Elsewhere it will have com- 
ponents of velocity parallel and perpendicular to 
the field. Thus we may say that the magnetic 
forces act effectively over only half of the orbit 
and the factor 2 in the denominator of the 
gyromagnetic ratio takes account of this fact. 
Now in the bipolar coordinate model of the elec- 
tron the orbit is defined only at two diametrically 
opposite points of its orbit, so that the parts of 
the orbit over which the magnetic forces would 
vanish are missing. Hence the factor 2 should be 
omitted here, and the magnetic moment should 
be given by 

vo= (e/mc)s, 


or for some component 


u=eh/2mc. (40) 


The magnetic moment in e.s.u. should be the 


5 The considerations in this section are due to W. T. 
Payne. They were first deduced from a model similar to 
the bipolar coordinate model of the electron, but much 
less complete. 


charge of the electron times its velocity times 
the radius of the orbit. We have taken the 


velocity to be c and hence 
(a/2)=h/2mc. (41) 


This is what we have decided to call the radiys 
of the electron. It is (137/2) times the radiys 
deduced from the electrostatic energy of the 
electron. 

Our spin functions are double valued, and 
hence the wave-length to be associated with the 
spinning motion should be twice the circum. 
ference of the orbit or 


\=4n(a/2) = 24h/mc=h/mce, (42) 


Cuo=ceh/2mc=eca/2: 


which is the Compton wave-length. 


THE RELATIVISTIC PROPERTIES OF BIPOLAR 
COORDINATES 


The relativity transformations include the 
time as well as the usual space variables. This 
means that we shall have to examine the intro. 
duction of time into the picture rather carefully, 
In introducing the time in the previous sections 
we assumed implicitly that the points on the 
sphere of radius $a did not all have the same 
proper time, for a point is of importance only 
when we can imagine the electron passing 
through it, and at a given instant only one point 
may be occupied by the electron. We shall take 
the proper time of the right-hand reference point 
to be zero and that of the left-hand reference 
point to be 2r. This means that the time required 
to travel from the right-hand to the left-hand 
reference point at velocity c is 27 where the 
radius is taken as }a=cr as before. The proper 
times of the points elsewhere on the sphere will 
be between zero and 27. Our choice of ¢=0 at 
N is a little arbitrary; it would do equally well 
to take t=0 at S and ¢=2r at N. This alternative 
choice will be discussed later. 

The bipolar coordinates of a point on a sphere 
will have exactly the same spatial characteristics 
as before, and we need not repeat the discussion 
for rotations. We may, therefore, pass directly 
to Lorentz transformations. It may be show 
immediately that the bipolar coordinates of a 
point do not transform like Dirac spin functions 
under a Lorentz transformation. However, 4 
quantity very closely connected with them does 
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We shall consider relative motion in the 2 
direction only since any other motion may be 
obtained by combining this with the appropriate 
rotation. We have 


x 2’=zcoshyu+ct sinh yp, 
y 


ct’ =z sinh w+ct cosh yp, (43) 


x’ 
y’ 


where cosh yw and sinh uw are as defined in the first 
section. Equation (43) may be written as 


x’ =x (2’+ct’) = (s+ct)e, 


y’=y (—2'+ct')=(—2+0cl)e*. (44) 


These equations hold for all z and ct whether they 
are represented on the sphere or not. In par- 
ticular 

(—2'+cr’) =(—2+cr)e*, 

(2’+cr’) =(s+cr)e 

holds for all z where cr is the radius of the sphere. 
Of course, the point (z, cr) occurs on the sphere 
only at z=0. We have the equations 


z=(|v|?—|u|*)/2a, 
cr=}a=(|v|*+|u|*)/2a 
and hence 


(45) 


(46) 


s+cr=|v|*/a. 
Substituting in (44), we obtain 

u'//a’ =(u//a)e™, 

v'/V/a' =(v/V/a)e™, 


where a’ =2cr’. u/\/a and v/\/a rather than the 
bipolar coordinates transform like a pair of 
Dirac spin functions. These quantities are the 
square roots of the distances between the inter- 
section of the normal from the point on the sphere 
to the z axis and the reference points. 

The transformation of a point on the sphere is 
shown in Fig. 5. For a positive velocity of the 
primed system with respect to the unprimed 
system, e*>1, so that BC is shortened and AC 
lengthened in the transformation. 


—2+cr=|u|?/a, (47) 


(48) 


REFLECTIONS AND FOUR COMPONENT 
SPIN FUNCTION 


It is easily verified that under a space inversion 
(49) 


u’=v*, vo = —u*. 
Let us change the notation slightly and take 


uU=U, V=?}, u’ = Ue, v' =02, (50) 


where now 4, and 2, are the bipolar coordinates 
of some point on the sphere and uw and 22 are the 
bipolar coordinates of the diametrically opposite 
point. If now we call u,, v1, u2, and v2 the four 
functions to be associated with a given point, we 
may obtain linear transformations for reflections 
in space. The matrix for the transformation 
under an inversion is easily seen to be identical 
with that given in the first section, Eq. (4). Now 
it will be recalled that the bipolar coordinates 
are double valued in all of their variables and, 
therefore, we may associate two laps with any 
rotations. The fact that the 6’s appear with the 
same sign in U; means that we are to think of 
the two points being in the same lap in the path 
of the electron. 

Let us now consider the time reversal. We 
have so far assumed that ty=0 and tgs=2r. A 
time reversal means that we take ty=0 and 
ts= —2r. Now if ts= —2r, we must think of the 
electron as moving from S to N instead of from 
N to S as previously. Evidently, then, NV and S 
interchange their roles as in the space inversion. 
In addition, however, one of the two points at 
opposite ends of a diameter changes laps because 
the motion is taking place at essentially 27 
seconds earlier. If originally the electron was 
moving through P on its way to S from N and 
now it is moving through P’ on its way to V 
from S, 27 seconds later it will be at its original 
position at its original time. But as it goes 
through N it must change laps. On the other 
hand, if the electron were originally at P’ and 
the time reversal were performed it would find 
itself at P 27 seconds earlier, but it would be on 
the same lap since it passes only through S and 
not through JN before arriving at the original 
position. Now the bipolar coordinates of a given 
point differ only by sign in the two laps, and 
hence the transformation is to be given by Ur 
in Eq. (4). 

Evidently uz and v2 transform just like u, and 
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v; under rotations, and it remains to investigate 
how they will transform under a Lorentz trans- 
formation. “2 and v2 will be associated with a 
point whose z coordinate is the negative of the 
point associated with u, and 2, and hence if we 
change the sign of a in Eq. (47) and replace u 
and v by wu: and v, respectively, we shall obtain 
the correct transformation. This gives 


ta! //a' = (u2/r/a)e™, 


v2'/V/a’ = (02/V/a)e™. -_ 
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The Lorentz transformation for the set of func. 
tions “1, 21, us, and v2 is thus U, in Eq. (5), and 
we have complete agreement with the Dirac 
theory as far as transformation properties go. 

In conclusion the author would like to express 
her deep appreciation of the many helpful and 
encouraging discussions with Professor Leop 
Brillouin, under whose direction the work wag 
carried out, and with Dr. William T. Payne, to 
whom she is also indebted for the use of his yp. 
published work on the subject. 
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Adsorption of Thorium on Tantalum* 


CHARLES J. 


GALLAGHER 


Department of Physics, University of Notre Dame, Notre Dame, Indiana 
(Received September 30, 1943) 


It is known that thorium adsorbed on the surface of tungsten or molybdenum results in a 
decrease in the work function and also a decrease in the value of A, the emission constant. 
The behavior of thorium on tantalum was investigated, and similar effects were noted. Thorium 
was evaporated from thoriated tungsten onto a well-outgassed tantalum ribbon, and the 
thermionic emission at 1390°K was studied as a function of the time of evaporation. The 
work function of the ribbon changed from 4.07 volts to 2.52 volts, while A decreased by a 
factor of thirty. The ribbon was then held at 1550°K, and the emission from the two sides 
was studied separately. The emission from the activated side decreased, while that of the 
clean side increased until both reached the same final level after about twenty hours. This 
indicates that the adsorbed thorium had migrated around the edges of the ribbon and dis- 
tributed itself uniformly over the emitting surface. 


INTRODUCTION 


HE thermionic properties of metals are 
highly sensitive to the presence of adsorbed 
surface layers. For example, thorium adsorbed 
on tungsten or molybdenum is known to produce 
very large increases in the thermionic emission. 
Similar results are reported for other electro- 
positive elements deposited on these metals.' It 
is also known that these adsorbed atoms are able 
to migrate over the underlying surface in the 
direction of decreasing concentration at a rate 
dependent on the temperature. Evidence for this 
migration effect has been obtained by measuring 


* Submitted to the Committee on Graduate Study of 
the University of Notre Dame in partial fulfillment of the 
requirements for the degree of Doctor of Philosophy. 

1A. L. Reimann, Thermionic Emission (John Wiley & 
Sons, 1934), pp. 99, 134. 


the emission from different parts of the surface! 
and by taking photographs of projection tube 
patterns.*-® 

It was the object of the present experiments 
to determine whether thorium would produce 
similar effects if adsorbed on a surface of tanta. 
lum. A tube similar to that used by Brattain 
and Becker? was designed to permit separate 
measurement of the thermionic emission from 
either side of a tantalum ribbon. Thoriated 
tungsten filaments were mounted so that thorium 





os... H. Brattain and J. A. Becker, Phys. Rev. 43, 428 
(1933). 

+E. Briiche and H. Mahl, Zeits. f. tech. Physik, No. 12, 
623 (1935); also, No. 8, 262 (1936). 

*M. Benjamin and R. O. Jenkins, Proc. Roy. So 
A176, 262 (1940); also, A180, 225 (1942). 
; a J. Ahearn and J. A. Becker, Phys. Rev. 54, 48 
1938). 
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could be evaporated onto either side without 
being deposited on the other. Thus, either side 
of the ribbon could be activated independently, 
and its thermionic properties studied for different 
amounts of adsorbed thorium. The migration 
effect could be investigated by activating the 
one side and comparing its emission with that 
of the other side, holding the ribbon at a suitable 
temperature. 

The relation of the saturated emission current 
density J to the absolute temperature T is given 
by the familiar Richardson equation 


I=AT? exp (ge/kT), (1) 


where ¢ is the work function of the surface, A 
is the emission constant, and e and k are the 
electronic charge and Boltzmann’s constant, 
respectively. 

Theoretically, A is a universal constant with 
a value of 120 amp. per cm? per deg.’, while ¢ is 
a constant for a given surface. However, experi- 
ment shows that A is also dependent on the 
nature of the surface, particularly if impurities 
are present. In general, a deposit of an electro- 
positive film on a metallic surface produces 
decreases in both A and g¢.' This increases the 
thermionic emission since the exponential factor 
in Eq. (1) predominates. 

Both of these constants can be obtained for a 
surface if the temperature versus emission char- 
acteristic is known. If we take the common 
logarithm of each side of Eq. (1) and rearrange 
terms, we obtain 


log T2/I = ye/2.3kT —log A. (2) 





Fic. 1. Electrode arrangement. A, B’—Half-ring col- 
lectors; length—1.27 cm, diameter—1.91 cm. C, C’— 
Shields; length—2.54 cm, width—1.27 cm. D, D’—Guard 
rings; length—2.54 cm, diameter—1.91 cm. W, W’— 
Tungsten filaments; length—6.00 cm, diameter—0.076 
mm. Ta—Tantalum ribbon; length—6.00 cm, width— 
0.159 cm, thickness—0.00293 cm. 
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Fic. 2. Richardson plots. 


By plotting the experimental values of log (T?/J) 
against (1/7), we obtain a straight line with 
slope (ve/2.3k) and intercept (—log A) on the 
log (T?/I) axis. The value of ¢ in volts is given 
by the relation 


y = 0.862 X 10 Xslope. 


THE EXPERIMENTAL TUBE 


The general features of the tube are shown in 
the diagram of the electrode arrangement (Fig. 
1). The emission from opposite sides of the 
tantalum ribbon Ta was collected by the semi- 
cylindrical plates A, B. Electrostatic shielding 
between the two interelectrode spaces was pro- 
vided by the flat plates C, C’. Cylindrical guard 
rings D, D’ were mounted coaxially with the 
collectors to insure a uniform field along the 
length of the ribbon. The parts of the collecting 
system were all of preglowed nickel and were 
connected by tabs of outgassed Fernico to a 
supporting structure of Pyrex and No. 704 glass. 

The structure was mounted in a Pyrex blank, 
which was first evacuated and baked for twenty- 
four hours at 500°C. The whole assembly was 
then baked on the pumps for eighteen hours at 
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Fic. 3. Activation curves. 






450°C. This was followed by vigorous heat 
treatment of the metal parts, carried out over a 
period of fifteen hours, after which the tube was 
sealed off. 

A typical heating schedule followed for out- 
gassing is given next. The tungsten filaments 
were glowed at 2000°K for twelve hours, flashed 
to 3000°K for two minutes, and dropped to 
1700°K until after the seal-off. At the same time, 
the tantalum was flashed to 2350°K for eight 
periods of about five minutes each and baked at 
1750°K between flashes. During each such flash, 
the collector system was bombarded by electrons 
from the ribbon. A barium-aluminum getter was 
exploded in a side tube an hour before the seal-off. 
All plates and filaments were kept glowing while 
the tube was removed from the vacuum system. 


















MEASUREMENT OF EMISSION AND 
TEMPERATURE 







The tantalum emission current was read on a 
sensitive galvanometer with a range extended by 
| a suitable shunt from 5X10-' ampere to 
1.25X10-* ampere. All emission measurements 
were made with a potential of 100 volts between 
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the ribbon and the elements of the collector 
system. The filaments were heated by alternating 
current fed through isolation transformers to 
permit grounding the galvanometer circuit. 

The temperature of the ribbon was obtained 
with a Leeds and Northrup disappearing filament 
pyrometer. The brightness temperatures so found 
were converted to true temperatures by use of 
the temperature scale of Malter and Langmuir.* 
Over the range of temperatures used in the 
emission studies (1100°-1700°K), the pyrometer 
settings could be reproduced to within four 
degrees. 


TANTALUM CLEAN-UP 


Three attempts were made to evaporate 
thorium onto the tantalum ribbon. In each case, 
before evaporation was begun, the tantalum was 
subjected to a prolonged and vigorous heat 
treatment. This consisted of five-minute flashes 
to 2350°K, followed by baking periods of several 
hours at 1750°K. The thermionic emission was 
measured at intervals, and the treatment was 
continued for an average of thirty hours until 
the emission reached a stable value. 

The Richardson lines (Fig. 2—IA, IIA, IIIA) 
obtained at the end of each clean-up period gave 
the following values for ¢: 


I[A—4.52 volts, 
ILA—4.08 volts, 
IIITA—4.04 volts. 


The last of these should be the most reliable 
since it was obtained after the third period of 
outgassing. It compares favorably with the value 
of 4.07 volts reported by Dushman.’ In each 
case, identical values of emission were found for 
both sides at all temperatures. This agreement 
is a good check on the symmetry of the structure. 


ACTIVATION OF TANTALUM 


After each clean-up period, thorium was 
allowed to evaporate from one of the tungsten 
filaments and condense on the near side of the 
ribbon. In order to prevent the ribbon from 


taking up gas given off by the tungsten, a 


6 L. Malter and D. Langmuir, Phys. Rev. 55, 743 (1939). 
7S. Dushman, H. N. Rowe, J. Ewald, and C. A. Kidner, 
Phys. Rev. 25, 338 (1925). 
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special heating schedule was used for evapora- 
tion, similar to those used by Coomes® and 
Daniel® in depositing thorium on tungsten. It 
included these three steps: 

1. The tantalum was set at 1600°K. 

2. The tungsten was flashed for a minute or 
more to about 2900°K and then dropped to an 
evaporating temperature of about 2500°K. 

3. The ribbon was cooled to room temperature 
and held there for deposition of the thorium. 
The initial temperature of 1600°K was chosen 
for the ribbon because preliminary studies in a 
simple projection tube indicated that tantalum 
outgasses very readily in the range 1600°K to 
1800°K. Thus, if held at 1600°K, the ribbon 
would not be likely to take up gases released in 
flashing the tungsten. 

The results for each of the activation runs are 
shown in the curves in Fig. 3. These are loga- 


rithmic plots of the emission current density in ° 


amperes per cm? at 1390°K versus time of 
evaporation of the thorium. Although the runs 
vary among themselves in the rate of change of 
emission and in the final value reached, yet each 
shows a definite increase in the emission due to 
the presence of the thorium. 

The abrupt termination of the first run (Fig. 3, 
I) was due to the burning out of the tungsten 
filament after ten hours of evaporation. The 
general appearance of the curve suggests that a 
much greater increase in emission would have 
resulted had it been possible to extend the time 
of activation. The second run (Fig. 3, II) was 
carried out for twenty hours and the emission 
increased by about one order of magnitude. This 
small change was probably caused by the forma- 
tion of large crystals in the tungsten. It is 
known that large grain filaments are very diffi- 
cult to thoriate.5 

In the third run (Fig. 3, III), a steady increase 
in the emission was noted over a period of 
seventy hours, with a final value four orders of 
magnitude above the original level. The filament 
burned out after the last reading was taken, so 
that it is not certain whether the leveling-off 
near the end of the run represents the condition 
of optimum thoriation. 


SE. A. Coomes, Phys. Rev. 55, 519 (1939). 
* J. A. Daniel, Phys. Rev. 61, 657 (1942). 
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The sharp dip in the curve after seventeen 
hours was probably caused by the migration of 
some of the thorium around the edges of the 
ribbon to the other side. The temperature of the 
ribbon was inadvertently left at 1434°K for 
about two hours, during which the decrease 
occurred. A check on the emission from the other 


_side revealed a sharp increase, which is sup- 


porting evidence for the view that migration was 
responsible for the effect. 

At the end of each run, Richardson plots were 
made for the surface and values of the work 
function calculated. The lines in Fig. 2 give the 
following values for ¢: 


IB—3.60 volts, 
IILB—3.80 volts, 
III B—2.52 volts. 


Each of these represents.a considerable drop 
from the value obtained for the clean surface. 
The last value is probably very nearly equal to 
that for optimum thoriation. 

The value of A was computed for this last 
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Fic. 4. Surface migration curves. 
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TABLE I, Thermionic constants. 











Mo Ta 
¢: (Clean) 4.54 4.15 4.07 volts 
¢2 (Thoriated) 2.63 2.59 2.52 volts 
Ai/A2 20-33 37 31 








run for both the initial and final stages: 


IILA—A =15.5 amp./cm*/deg.? 
IIIB—A = 0.5 amp./cm?/deg.?. 


These cannot represent the true values for the 
two surfaces since there is considerable uncer- 
tainty in assigning an exact value to the effective 
emitting area of the ribbon. The shields draw 
some of the emission current, so that the effective 
area would be less than that calculated from 
the dimensions. This would increase the values 
of A given above. However, since the area is the 
same for both cases, the ratio of these numbers 
is independent of the actual area. The ratio of 
the initial A to the final A isseen to be thirty—one. 
It is of interest to compare this ratio with the 
corresponding ratios obtained for thorium on 
tungsten and thorium on molybdenum. For 
tungsten, the value is about twenty; for mo- 
lybdenum, a value of thirty-seven is obtained. 
When the large experimental uncertainties are 
considered, these ratios are not far apart. The 
results indicate that the nature of the surface 
film of thorium is the same in each case. To 
emphasize the similarities for these three cases 
of thoriated metals, the values of the work 
functions for the six different surfaces and the 
values of the ratios of the A’s for the three 
transitions are listed in Table I. It should be 
noted that the values for thoriated tantalum do 
not necessarily represent the condition for opti- 
mum activity although the appearance of the 
activation curve suggests that the maximum 
emission point had almost been reached. 
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SURFACE MIGRATION 


The effects of surface migration on the thermj. 
onic activities of the two sides of the ribbon are 
shown in Fig. 4. These curves were obtained 
after the second attempt to activate the ribbon. 
They show how the emission at 1390°K changed 
with time when the tantalum temperature was 
held at 1550°K between readings. The upper 
curve (I) represents the decline of the emission 
from the originally activated face, while the 
lower (II) shows the growth of the activity of 
the opposite side. The two sides reach the same 
final level in about twenty-four hours. 

These curves indicate that some of the thorium 
had traveled around the edges of the ribbon and 
distributed itself over the opposite side. The 
final limiting value of the emission corresponds 
to a uniform concentration of thorium over the 
entire emitting surface. 


CONCLUSION 


These experiments show definitely that the 
thermionic emission from tantalum increases 
considerably when sufficient thorium is adsorbed 
on the surface. They indicate that the thorium 
film is of the same nature as the films on thoriated 
tungsten and thoriated molybdenum. A further 
similarity is noted in the phenomenon of surface 
migration. Thorium atoms adsorbed on the 
surface of any of these metals can migrate over 
the surface at temperatures of about 1600°K. 

It would be of interest to extend the activation 
curves for tantalum to determine the work 
function for optimum thoriation more precisely, 
It would also be of value to obtain more migra- 
tion curves at various temperatures. 

The author gratefully acknowledges the aid of 
Dr. E. A. Coomes, who suggested this problem, 
and Dr. E. Guth, who contributed many valuable 
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The Hydrogen Content of White Dwarf Stars 


GEORGE CHERTOCK 
760 Princeton Place, N. W., Washington 10, D. C. 


(Received November 10, 1943) 


A model for the white dwarf star, Sirius B, is considered in which the whole star is a com- 
pletely degenerate gas; the central core containing 11 percent of the mass and extending to 
21 percent of the radius has no hydrogen, and the remainder of the star has 45 percent hydrogen. 
This model is mechanically stable and satisfies the boundary values for mass and radius, but 
the energy production in the outer layers due to the reaction H+H=D-+e* far exceeds the 


observed luminosity. 





I. INTRODUCTION 


HITE dwarf stars are now believed to 

represent a final stage in the evolution of 

a star, a stage in which all major energy sources 

are exhausted and gravitational contraction has 

compressed the interior of the star to the state 
of a degenerate electron gas. 

The theoretical equilibrium equations for the 
interior of a degenerate gas star gives the radius 
as a function of u., the mean molecular weight 
per electron.' But at the high densities present 
in degenerate gas stars wu, will depend only on 
the hydrogen content X. For at these densities 
(~105 g/cm’), pressure ionization must be com- 
plete, and for all elements but hydrogen u.=2, 
while for hydrogen u.=1. Hence 


For the star Sirius B, if uw. is assumed constant 
throughout the star, then u, must equal 1.43 and 
X =0.40 in order to satisfy the observed values 
of mass and radius.’ 

But as shown by Wildhack,* and Marshak,' 
the core of Sirius B cannot have any appreciable 
hydrogen content, for the energy production 
due to the reaction H+H=D-+e?*, even at the 
zero-point energy of completely degenerate star, 


1S. Chandrasekhar, Introduction to the Study of Stellar 
Structure (University of Chicago Press, 1939), Chapter 11. 

* Chandrasekhar gives ue=1.32, X =0.52, although he 
uses the same observational data and the same mass-radius 
relation. His value seems to be a numerical error. 

*W. A. Wildhack, Phys. Rev. 57, 81 (1940). 

*R. E. Marshak, Astrophys. J. 92, 321 (1940). 
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would exceed 10 ergs/g sec. while the observed 
energy production is only 6X10~— erg/g sec. 

This paper investigates the possibility that 
there is no hydrogen at all in the core of Sirius B 
but that it is present in the outer layers. We 
calculated the size of the core and the hydrogen 
content in the outer layers which will make the 
star mechanically stable. 

This proposed composite model for Sirius B is 
somewhat plausible because of its supposed mode 
of evolution from the main sequence to the 
degenerate state. In the main sequence the 
energy source is a thermonuclear reaction in 
which H is changed to He with C and N as 
catalysts.® The reaction presumably starts at the 
center, and then proceeds outward, thus forming 
an energy-producing shell of gradually increasing 
diameter. Inside the shell the matter is isothermal 
without hydrogen; outside the shell, there is 
hydrogen, and the distribution of density and 
temperature corresponds to that in the ordinary 
point source model. 

But it was shown by Schénberg and Chandra- 
sekhar® that such a model is at all possible only 
if the core contains not more than 10 percent of 
the star’s mass. This means that after 10 percent 
of the star’s hydrogen is exhausted, the star 
must either continue to evolve through non- 
equilibrium configurations, or somehow cool 
down, stop the thermonuclear reaction, and 
continue to evolve by gravitational contraction. 
In either event it seems plausible that the 
degenerate star which eventually results has a 
core without hydrogen. 

5H. A. Bethe, Phys. Rev. 55, 434 (1939). 


*M. Schénberg and S, Chandrasekhar, Astrophys. J. 
96, 161 (1942), 
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Fic. 1. Mass vs. radius for core and envelope of Sirius B. 











II. CALCULATIONS OF A COMPOSITE MODEL 


The equilibrium equations for a degenerate gas 
star can be written! 


. 
r? dri. dr a 2A ” 
arm'‘c® 





A= = 6.01 X 10”? c.g.s. units, 
3h* 


8armicr HT 


3h* 


= 9.82 105 c.g.s. units, 


where the mass interior to radius 7 is 





—8A dy 
M(r)= fm, (3) 
GBu, dr 
the density at r is 
p=Bu.(y?—1)!, (4) 


and the pressure at r is 
P=A[y(2y?—5)(y?—1)! 
+3 sinh-! (y?—1)#]. (5) 


These equations were solved numerically for 
the structure of the outer, hydrogen-bearing 
layers of this proposed model. They were started 
under the boundary conditions P=0 and 
M=4M.», at the observed radius of Sirius B, and 
then were numerically integrated inwards by 
the method of successive approximations. Three 
separate calculations were made for the values 
of log M, log R, u., and X shown in Table I. 
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Fic. 2. Pressure vs. radius for core and envelope of Sirius B, 


The values of X were so taken because, as 
shown by Strémgren, the mean hydrogen content 
of main sequence stars varies from 25 to 45 
percent.’ The results of the integration are 
shown by the solid line curves in Figs. 1 and 2, 

For the hydrogenless core the equations must 
evidently satisfy the conditions that at r=0, 
dy/dr=0, and P=P(yo), where yo is the center 
point value of y. Chandrasekhar has numerically 
solved these equations for ten different values of 
yo by integrating outwards from the center. His 
results are given for u.=1, but they can be 
easily adjusted for u4.=2 (no hydrogen) and the 
accepted mass and y results are shown by the 
dashed line curves in Figs. 1 and 2 for three 
different values of yo. 

If now a particular core solution and a par- 
ticular envelope solution intersected at the same 
value of the radius in both Figs. 1 and 2, then 
this radius could be taken as the interface be- 
tween the non-hydrogen core and the hydrogen- 
bearing envelope, and the composite star would 
be mechanically stable; for the equation of 
state would be satisfied within the star, and the 


TABLE I. Boundary values used in integration for 
outer layers. 











No. log R/R(O) log M/M(O) Me xX 
1 —1.71* —0.01* 1.38 0.45 
2 —1.71* —0.01* 1.40 0.42 
3 —1.71* —0.02** 1.48 0.35 








* G. P. Kuiper, Astrophys. J. 88, 472 (1938). 
** H. N. Russell and C. Moore, The Masses of the Stars (University 
of Chicago Press, 1940). 


7B. Strémgren, Zeits. f. Astrophys. 7, 222 (1933). 





ma 
cor 
sur 
neg 


solt 
pro 
this 
of 7 


lear 
simi 
of a 


the 

they 
relat 
to 1 
Emc 


are § 








HYDROGEN CONTENT OF STARS 


mass and pressure would satisfy the boundary 
conditions at the center, interface, and outer 
surface. The effect of temperature is assumed 
negligible and is neglected. 

An inspection of the figures indicates that the 
solution for the envelope with yu,.=1.38 will 
probably intersect a core solution with 1/0? 
=0.45 at r=2.9X10* cm. The characteristics of 
this composite model are given in the first row 
of Table II. 

When the present work was completed, we 
learned that Chandrasekhar and Henrich made 
similar calculations for the same composite model 
of a degenerate star.* Their methods were the 
same as ours except that instead of integrating 
the exact differential equation for the envelope, 
they used the P=kp** approximation (non- 
relativistic). This approximation enabled them 
to use the F and M solutions to the Lane 
Emden equation with »=3.° 

Their results which most nearly fit Sirius B 
are given in the second row of Table II. 


Ill. CONCLUSIONS 


The above calculations have shown that it is 
possible to obtain a valid model of Sirius B with 
the central 11 percent of mass completely de- 
prived of hydrogen. It seems, however, that such 
possibility does not remove the paradox of high 
energy production since it is easy to show, that 
even in the outer parts of the model the tempera- 
ture will be high enough to produce rapid 
transformation of hydrogen. In fact, the decrease 
of temperature from the center towards the 


§ Private letter. See also Pub. Am. Astronom. Soc. 10, 
306 (1943). 
* See reference 1, p. 420. 


TaBLE II. Characteristics of composite model. 








Interface Core | Envelope 


log M(r)/ 
R/R(O) | Ms r/R 


—1.71 | .11 214 
—0.70 | 115 .22 


Star 


1/ye Me 


0.45 | 1.38 
0.50 | 1.33 


log 
M/M(O) 


—0.01 
0.035 





0.96 
1.28 

















surface of the degenerate model is so slow that 
even in the outer parts of the stellar body the 
conditions for the nuclear reaction will still be 
satisfied. According to Marshak‘ the temperature 
within 0.001R from the surface of Sirius B 
reaches 5X10®°C and the density reaches 10°. 
Now the rate of energy production for the 
proton-deuteron reaction is given by 


3.38 X 108 
€=4.6X10"%p7T— exp (- ——) (6) 


and at this density and temperature e=8.5 
X10~ erg/g sec. which is already in excess of 
6X10-*, the average energy production. Since 
the rate of energy production goes up expo- 
nentially with temperature and since the temper- 
ature will increase toward the interior, the 
integrated luminosity for the whole star will far 
exceed that observed. 

Thus it seems that the conclusion about the 
complete absence of hydrogen in the entire body 
of the star is quite unavoidable. The composite 
model although mechanically stable is still 
unsatisfactory on this score. 

I wish to express my gratitude to Professor 
G. Gamow for advice and assistance throughout 
this work, and also to Dr. S. Chandrasekhar for 
informing me of his calculations on the composite 
model. 
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The application of methods of geometrical optics for the 
investigation of electron trajectories requires the knowledge 
of the refractive index u of the equivalent optical problem. 
A general formula for » in terms of the electrostatic 
potential V and of the magnetic vector potential A is 
known. If A0 this formula contains a unit vector s 
parallel to the electron velocity, which is not known in 
general. It is therefore important that s be eliminated. 
This can be done only if a suitable integral of the equations 
of motion of the electron is known. If V and A are sym- 
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The Refractive Index in Electron Optics 
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metric about an axis, 8 may be eliminated from 4 by 
means of a momentum integral (Glaser). This elimination 
is generalized in the present paper to include cylindrical 
fields (not necessarily symmetrical) as well as some more 
general fields with a transversal field component. The 
result is reached by interpreting the magnetic scalar 
potential U as a velocity potential of an ideal fluid ang 
by making use of the corresponding stream function w tg 
obtain an expression for yu. 





















© te treat at the same time both the case of 
constant and relativistic mass we assume 
that the mass m of the electron is a function of 
the magnitude of its velocity v. The equation of 
motion in a charge free space may be written in 
any of the following three forms:! 


d(mv)/dt=F, (1) 
5 f ‘Lat=0, (2) 


Pi 
sf pds=0, where p=L/v. (3) 
Po 


The meaning of the symbols is obvious. F is the 
Lorentz force acting on the electron (which is 
here taken as a point mass and a point charge) ; 
L, the corresponding Lagrange function (or elec- 
trokinetic potential) : 


L= f modo— V+A-v, (4) 


where eE= —grad V; eH=c rot A; E and H are 
the intensities of the electrostatic and magnetic 
field, respectively; e the electron charge (which 
is to be taken with its sign); and c the velocity 
of the light in vacuum. Equation (4) is a formula 


1 Concerning the Hamilton’s principle for a mass which 
is a function of its velocity, see C. Somigliana, Rend. 
Accad. Lincei 13, 804-809 (1931) or Memorie Scelte 
(Lattes, Torino, 1936), pp. 623-629. 
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of Schwarzschild? and Levi-Civita* in a slight 
generalization of the writer*® for the case of 
variable mass. Equation (1) has an energy 
integral :* 5 







mst— { modo-+ V=K=constant, (5) 





so that by (3) and (4): 
p=mv+A-s—(K/v), 


where s=v/v. We use cylindrical coordinates 
(p, ¢, 2) and put v= gy—eg+ez, where é is either 
0 or 1. Let U be the magnetic scalar potential 
defined by the equation rot A= —grad U. We 
limit ourselves to the case in which U is a fune-. 
tion of p and » only, that is U= U(p, v), whereas 
the electrostatic potential is: 


V=g(p, vy) +9q(ee+2), 


where g is also a function of p and v only andg 
may be any constant. If e=1, the variable» 
coincides with z, and U, g become potentials 
symmetric about the z axis. If e=0, the variable 
v becomes ¢, and U and g are plane potentials. 
It is seen in this way that the quantities ¢ and? 
are introduced only for the purpose of a concise 
treatment of these two types of fields. Since 














2 K. Schwarzschild, Géttingen Nachrichten Math.-Phys 
Kl. 3, 126-131 (1903). Also F. Zerner, Handbuch da 
Physik, Vol. 12, p. 182. 

3T. Levi-Civita, Archiv f. Math. Naturvid. 31, 147 

1911). 
‘ 4 oe Comptes rendus 208, 638-640 (1939). 

5 I, Opatowski, ; ath. Phys. 20, 418-424 (1941). In 
Eqs. (5), p. 419 of this paper the condition uH should 
be added and the constant p should be taken equal to 
zero throughout the paper. 


























u by 
ation 
lrical 
more 
Calar 


and 
Ww to 











REFRACTIVE INDEX 


U(p, v) is a harmonic function it may be looked 
upon as a velocity potential of an ideal fluid 
motion. This is a well-known analogy between 
magnetism and hydrodynamics.® One may easily 
check that the stream function w(p, v) cor- 
responding to the velocity potential U(», v) satis- 
fies the equations :°~® 


dw/dp=p (ep? +1—€)dU/d», 
dw/dv= — pdU/dp. 


Various methods for the calculation of w are 
known. If U is expanded in series of the type 
U=R(p)N(v), which includes the classical ex- 
pansions in cylindrical, spherical, trigonometric, 
and other functions, w is obtainable from U by 
a simple operation of differentiation.” *'° 

Since>*® A=w grad (s—ez—eg), Eq. (4) may 
be changed into: 


L= f modv— V —(eg+ez—Z)w, (6) 


so that by (3) and (5) we have:" 
p=mv— (eo+ez—Z)vw—v'K. (7) 


Equation (1) has the following momentum in- 
tegral :° 


m(ep?o+ez— 2) =w+q(1—2e)t+M, 
(M=constant), (8) 


where ¢ means the time. The left-hand side of 
(8) represents a component of the momentum 
of the electron; it equals mv-p? grad ¢ in the 
axially symmetric case (e=1) and —mv-grad z 


See, for instance, M. Abraham, Enc. d. math. Wiss., 
Vol. IV, 2, I, pp. 16, 17. 

7H. Lamb, Hydrodynamics (Cambridge University 
Press, London, 1932), sixth edition, pp. 63, 125-128. 

SE. Trefftz in Frank-Mises, Differential- und Integral- 
gleichungen der Mechanik und Physik (Vieweg, Braun- 
schweig, 1935), Vol. 2, second edition, pp. 399-402, 446. 
oo — E. H. Love, Enc. d. math. Wiss., Vol. IV, 2, I, pp. 

, 101. 

107. Opatowski, Atti Accad. Sci. Torino 68, 135-146 
(1932-1933) ; Phil. Mag. 15, 849-850 (1933) ; Rend. Accad. 
Lincei 17, 1049-1054 (1933) and 18, 18-25 (1933); Mem. 
Accad. Lincei 6, 327-351 (1936). 

4 Equations (4) and (6) are not the only possible forms 
of the Lagrange function. In fact, the magnetic field does 
not change if we add to A the gradient of an arbitrary 
function ¥, because rot grad ¥=0. Therefore, since v- 
grad y =dy/dt = y we may add y on the right-hand side of 


(4) and (6) without changing the physical meaning of L. 
We get in this way from (6), if we take for instance 
¥=(ep+es—z)w, 

L= f{mvrdy— V+(eg+es—2)wv. 
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in the cylindrical case (e=0). Since «¢ is either 
0 or 1 the time derivatives in (7) may be elimi- 
nated by means of (8). This gives 


p=mv—wP(mv)"— Ko, 


where ©=(w—gt+M)/p* in the axially sym- 
metric case (e=1, y=z) and @=w+gqt+WM in the 
cylindrical case (e=0, v=g). It is seen in this 
way that the presence of the transversal field 
component (represented by the constant q) intro- 
duces time into the expression of the refractive 
index.” 

Since, by (2), we can take L+X instead of L, 
and therefore 4+(K/v) instead of u, we may 
write also 

w= mv— (mv)—dw. (9) 


The above formulae for the refractive index hold 
independently of any relationship between m and 
v. The velocity v which appears here explicitly, 
or through m, may be eliminated by means of 
the energy integral (5), which in the case of 
constant mass is :5 


°=2(K—V)/m (10) 
and for the relativistic mass becomes: 


v=c(l—moem) with m=(K—V)/ce. (11) 


If the mass is treated as constant, we may take 
as the refractive index the expression (9) divided 
by (2m). The elimination of v by means of (10) 
gives then 


u=(K—V)!—w/[2m(K— V)*). 
For relativistic mass we obtain from (9) and (11): 
u=c(m? —m¢)*—wb/[c(m* —m¢)*] 


with m=(K—V)/e. 
Taking into account the relation 


A=w grad (z—ez— eg) 


it is not difficult to check that the above equa- 
tions for u are equivalent to formulas which have 
been known in particular cases." 


2 This is not a first example of a time dependent re- 
fractive medium. The human eye is another example, cf. 
J. Picht, Ann. d. Physik 15, 927 (1932). 

13 See, for instance, P. Frank in Frank- Mises, reference 8, 
pp. 35, 36, 42, 43; W. Glaser, Zeits. f. Physik 80, 451-454 
(1933) and 81, 647-649 (1933); also in H. Busch-E. 
Briiche, Beitrige zur Elektronenoptik (Barth, Leipzig, 
1937), p. 29; E. Briiche-O. Scherzer, Elektronenoptik 
(Springer, Berlin, 1934), pp. 20-22, 38. 
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Comments on the Anisotropy of Metals 


H. S. CoLEMAN AND H. L. YEAGLEY 
Department of Physics, The Pennsylvania State College, State College, Pennsylvania 


(Received November 12, 1943) 


URING the past ten years we have been 

greatly concerned with such properties of 
metal systems as might be related to inter- 
diffusion in the solid state. It is the purpose of 
this paper to report a number of what might 
have been previously regarded as unrelated facts. 
Inasmuch as our interest has been primarily from 
an experimental point of view, it is our hope that 
these facts may be of interest to the theoretical 
physicist. 

Recently Coleman and Yeagley' described a 
new method of observing the interdiffusion of 
metals in the solid state at temperatures far 
below the melting points of the metals. This 
method consists of forming a bimetallic film on 
glass by vaporization methods and noting the 
change in optical properties as the metals inter- 
diffuse. Diffusion coefficients as determined from 
data obtained by this method were several orders 
of magnitude higher than those determined by 
older methods. A number of possible explanations 
have been considered in an attempt to explain 
this discrepancy. Among these were considered 
the possibility that the diffusion laws may be 
different for films, that the crystal structure of 
the films may not be normal, that the contact 
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Fic. 1. Crystal planes of face-centered cubic system (Cu). 


1H. S. Coleman and H. L. Yeagley, Trans. A. S. M. 31, 
No. 1 (1943). 
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between the metals may be more intimate than 
in older methods, and that there may be some 
contamination during the formation of the 
bimetallic film that would accelerate the inter- 
diffusion by catalytic action. All of these except 
the question of crystal structure were eventually 
eliminated. Special studies in which x-ray diffrac- 
tion was used were then made on the films, by 
Coleman, Burr, Davey, and Yeagley.? Their 
results indicate that a large number of metals 
condensed on glass at room temperature have 
the same crystal structure as do the metals in 
bulk except that certain crystal planes are 
preferentially oriented with respect to the glass 
surface. It is found that the (111) planes for 
face-centered cubic metals are oriented parallel 
to the glass surface on which they are deposited. 
In the case of body-centered cubic metals, the 
(110) planes are oriented parallel to the glass 
surface on which they are deposited. Metals 
having the hexagonal close-packed crystal struc- 
ture are found to have their (00.1) planes 
oriented parallel to the glass surface on which 
they are deposited. Thus the anomaly in rates of 
interdiffusion of films appears to be a matter of 
anisotropy. 

Further consideration of this matter indicates 


100 001 


Fic. 2. Crystal planes of body-centered cubic system 
(8-brass). 


2Coleman, Burr, Davey, and Yeagley, to be published 
in Trans. A. S. M. (January, 1944). 
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that there may be some general physical principle 
involved since it is found that similar anisotropy 
exists in a number of other cases. For example, 
it is found that the anisotropy of magnetization 
is related to these same planes ;* also that in the 
case of plastic deformation, it is these planes 
that slip.‘ 

Since it was evident that the above crystalline 
planes have similar properties, it seemed ad- 
visable to see if these planes have similar 
structure. Large atomic models were carefully 
constructed, the face-centered cubic model hav- 
ing 172 atoms, the body-centered cubic model 
having 91 atoms, and the hexagonal model 
having 81 atoms. These models were mounted on 
special stands so that they could be viewed 
conveniently from any angle. When the three 
models were oriented, having the face-centered 
cubic (111) planes, the body-centered cubic (110) 
planes, and the hexagonal close-packed (00.1) 
planes in a horizontal position, it was quite 
apparent that the atoms were arranged in very 
similar patterns. It was also evident that these 
planes have the greatest atomic density. Sections 
of these planes are shown in Figs. 1-3. It is 
evident that the special planes mentioned are 


3F. Bitter, Introduction to Ferromagnetism (McGraw- 
Hill Book Company, Inc., New York, 1939), Chapter VI. 

‘F, Seitz, The Physics of Metals (McGraw-Hill Book 
Company, Inc., New York, 1943), page 75. 
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Fic. 3. Crystal planes of hexagonal close-packed system 
(Zn). 


made up of equilateral and isosceles triangles, 
forming hexagonal figures. 


CONCLUSIONS 


The following conclusions are suggested by 
the above facts: 

(1) The anisotropy in the diffusion rate of one 
metal into another indicates that orientation 
may be the primary factor in controlling this rate. 

(2) There is an intimate relation between the 
crystal structure, plastic deformation, ferro- 
magnetism, growth nucleation of metals, and 
intermetallic diffusion in the solid state. 
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—— publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. The closing date for this department is the 
third of the month. Because of the late closing date for the sec- 
tion no proof can be shown to authors. The Board of Editors 
does not hold itself responsible for the opinions expressed by 
the correspondents. Communications should not in general ex- 
ceed 600 words in length. 





‘The Structure of 25461 of Mercury (198) 


Jacos H. WIENS 
Department of Physics, University of California, Berkeley, California 
December 7, 1943 


INCE the first announcement of the production of the 

mercury isotope of mass 198,! it has been possible to 
produce much larger quantities of the single isotope. The 
process is essentially the same as described before except 
that a one-ounce gold sample was bombarded by stray 
neutrons from the 60” cyclotron for a period of 10 months. 
The discharge tube was 4 cm long and had an inside 
diameter of 0.3 mm, and the mercury discharge was not 
diminished in intensity in 30 minutes of operation. A blast 
of air was used to cool the tube. 

Two Fabry-Perot etalon spectrograms of the line \5461, 
one greatly overexposed, are shown (Fig. 1). No structure 
is visible in the spectrogram of mercury (198) although the 
structure is clearly visible in the mixed mercury. From the 
relative position of the line due to mercury (198) and the 
main component of the mixed mercury, it will be possible 
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Fic. 1. Fabry-Perot etalon spectrograms of the line 45461 of mer- 
cury (198) and mixed mercury. A, normal exposure. B, mercury (198) 
greatly overexposed. ° 


to verify Schiiler’s* assignments of the components of the 
line 45461. 

Silvered etalons with a separation of 9.015 mm were 
used, and the exposure time for the mercury (198) in A 
was 30 seconds, and in B was 2 minutes, while the time for 
the mixed mercury was kept constant at 20 seconds. The 
slit was moved in order that the plate need not be disturbed 
between exposures, and the center of the etalon pattern 
falls exactly between the two patterns. 

A complete report covering the production of the mer- 
cury isotope mass 198 will soon be published and this 
report will include the work being done at the present on 
a mercury source to be used by the Bureau of Standards 
in determining the wave-length of a new primary standard 
of wave-length.? 

1 J. H. Wiens and L. W. Alvarez, Phys. Rev. 58, 11, 1005 (1940). 


2 H. Schiiler and J. E. Keyston, Zeits. f. Physik 72, 430 (1930). 
3W. F. Meggers, Rev. Mod. Phys. 14, 2-3, 59 (1942). 





Proceedings of the American Physical Society 


MEETING AT EVANSTON, ILLINOIS, NOVEMBER 12-13, 1943 


HE 257th meeting of the American Physical 

Society occurred at Northwestern Univer- 
sity, Evanston, Illinois, on November 12 and 13, 
1943. This deviation from our traditional policy 
(of meeting at Thanksgiving-time at the Univer- 
sity of Chicago) proved a singularly happy one: 
travelling conditions were probably better than 
at the holiday week end, and our host institution 
threw open to us its magnificent new Tech- 
nological Institute for our sessions. There were 
slightly over two hundred registrations; we do 
not know how great was the fraction of these 
which was due to the Symposium hereinafter 


58 


mentioned, but are confident that it was notable. 
Professor B. J. Spence was head of the Local 
Committee. 

An outstanding feature of this meeting was 
the Symposium on the Physics of Rubber and 
Other High Polymers, consisting of sixteen in- 
vited papers of which the speakers and the titles 
are listed near the end of these proceedings. The 
credit for this notable event belongs almost 
entirely to Warren F. Busse, to whom authority 
was delegated to arrange it. It was an innovation 
in the policy of our Society, and the welcome 
which was given it suggests that there will be 
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ample cause for similar symposia hereafter. A 
group of those in attendance petitioned the 
Council for a new Division of the Society, and 
appointed some of its members to advise the 
Council in this matter. 

There were thirty contributed ten-minute 
papers, the abstracts of which are printed here- 
inafter. 

The dinner of the Society was held on Friday 
evening, November 12, at the Georgian Hotel, 
Vice President Dempster presiding. The speakers 
were President F. B. Snyder of Northwestern 
University welcoming us; Henry Crew com- 
memorating our late Past Presidents J. S. Ames 
and W. F. Magie; Paul Klopsteg speaking of 
the new mansion of the American Institute of 
Physics; Warren F. Busse and R. Bowling Barnes 
speaking on topics pertinent to the Symposium. 

The Council met on Friday, its opening action 
being the election to membership of no fewer 
than 149 candidates, probably a record in the 
history of the Society. Petitions were received 
for the formation of a new Division of the Society 
including within its scope rubber, textiles, plas- 
tics, and other high polymers. The request was 
granted, and the President was empowered to 
appoint an Organizing Committee composed of 
members of the Council and of sponsors of the 
new Division; this Committee shall define the 
scope, choose the title, prepare by-laws, and 
propose initial officers of the Division, submitting 
these recommendations to the next meeting of 
the Council. 

Reports reaching the office of the Society 
attest the loss by death of the following Fellows 
and Members: John Almy, J. S. Ames (former 
President), W. E. McElfresh, A. P. R. Wadlund, 
Irving Weiner. The report of the death of C. M. 
Kilby, quoted in a previous issue of these Pro- 
ceedings, was fortunately fallacious; apology is 
expressed to Professor Kilby for having given it 
currency. 

Elected to Membership: Alexander, Samuel N.; 
Alpher, Ralph Asher; Anderson, R. T.; Atherton, 
Ralph; Baldwin, Thomas Watson; Banca, M. 
Charles; Barrett, E. Otto; Bennon, Meyer; 
Benson, Allan Ingvald; Birks, Laverne S., Jr.; 
Boettner, Edward A.; Breed, Henry Eltinge; 
Brinster, John F.; Broyles, Arthur A.; Bruck, 
George G.; Burkig, Jack W.; Bush, Robert R.; 
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Cazenavette, Lionel L.; Chako, Nicholas; Col- 
lins, E. Hobart; Conant, Floyd S.; Coogan, 
Charles Henry, Jr.; Culp, Eugene H.; Davis, 
Luther, Jr.; DeLany, B. P.; Drukey, Donald L.; 
Duerig, William H.; Duroux, Dorothy M.; Egy, 
Willard L.; Engelkemeir, Donald W.; Evans, 
John; Fagan, Henry D.; Falk, Charles E.; Fein- 
berg, Theodore; Ferguson, John P.; Fox, Ken- 
neth Russell; Franklin, Philip J.; Fuchs, Robert 
A.; Gates, Halbert F.; Gernsback, Hugo; Ginc- 
berg, Arthur S.; Goodman, Bernard; Gregory, 
Christopher; Griffith, Gordon Lamar; Grimm, 
Henry H.; Gooze, Irving; Gran, W. H.; Grosse, 
Aristid V.; Gruner, Wayne; Hale, Prentice D.; 
Hall, George L.; Hartig, Lucy Acquaviva (Mrs.); 
Harvalik, Zaboj Vincent; Hazen, Wayne E.; 
Helmick, Paul S.; Hennessey, James J., S.J.; 
Hess, David C., Jr.; Hinch, William H.; Hoff- 
man, John W.; Hurd, Frank W.; Huston, Nor- 
man Earl; Hyde, Walter Lewis; Jacobs, Julian S.; 
Jaffey, Arthur H.; Johnson, Carlton F.; Johnson, 
Charles F., Jr.; Katz, Robert; Kelsey, Robert H.; 
Kelso, John Morris; Kessler, Karl G.; Kitchen, 
Sumner W.; Klein, Martin Leroy; Kolstad, 
George A.; Kubitschek, Herbert E.; Kuppen- 
heim, Hans F.; Lawrence, Jack; Leland, Wallace 
T.; Lengyel, Bela A.; Levinthal, Elliott Charles; 
Lewis, Margaret N.; Long, Edward Crawford, 
Jr.; Lufey, Carroll W.; Marton, Claire (Mrs.); 
Marvin, Marguerite; Marvin, Philip Roger; 
Mason, C. LeRoy; Mates, James; McLaren, Ian 
H.; McNeil, Edward; Menius, A. C., Jr.; Meyer, 
Charles F.; Meyer, Herbert M.; Meyerhof, 
Walter E.; Miller, Victor A.; Mooers, Calvin 
Northrup; Mutter, Walter E.; Orden, Alexander; 
Orr, Rollin S.; Oster, Gerald; Packard, Karle 
Sanborn, Jr.; Page, George V.; Pease, Robert L.; 
Pepiot, Harold M.; Perry, George L.; Pochapsky, 
T. E.; Prentky, Peter; Purcell, Edward M.; Reed, 
Allan W.; Rifenbergh, C. Morgan; Roberts, 
James Herbert; Roberts, Louis W.; Rothschild, 
Robert Friend; Rouvina, James; Rueger, Lauren 
J.; Samuelson, Wallace H. E.; Saxon, David; 
Schlegel, Richard, Jr.; Schneerer, William F.; 
Scholl, Allen W.; Sheppard, C. Bradford; Smart, 
J. Samuel; Smith, Karl F.; Smith, Lester C.; 
Smith, Lyle Willard; Snyder, Ralph; Sorensen, 
Earl G.; Staub, Hans; Stechco, Allen C.; Sturm, 
William James; Sullivan, William A. (Rev.); 
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Sweet, Alex L.; Sydoriak, Stephen G.; Uchiya- L.; Wittig, William L.; Wollman, Seymour H.; 
mada, Hilda Castle; Van Driest, Edward Regi- Yalow, Abraham Aaron; Zilsel, Paul R. 

nald; Viscardi, John E.; Wada, Walter; Wang, 
Theodore Joseph; Wasserman, Robert Harry; 
Watts, Richard John; Welton, Theodore A.; 
West, S. F.; Wheeler, Samuel Crane, Jr.; Willig, ° 
Frank J.; Winslow, George H.; Wiser, Herbert 


Kari K. Darrow, Secretary 
American Physical Society, 
Columbia University, 
New York 27, New York 


SYMPOSIUM ON THE PuysiIcs OF RUBBER AND OTHER HIGH POLYMERS 


(in order of presentation) 


H. M. James, Purdue University: Present Status of Theory of Elasticity of Rubber. 

L. A. Woop anp F. L. Rotu, National Bureau of Standards: Stress-Temperature Relations in Natural Rubber. 

H. A. BRAENDLE AND W. B. WIEGAND, Columbian Carbon Company: GR-S an Elastically-Inverted Polymer. 

J. H. Ditton, I. B. PRETTYMAN, AND G. L. HALL, Firestone Tire and Rubber Company: Hysteretic and Elastic Properties 
of Rubberlike Materials under Dynamic Shear Stresses. 

M. Mooney, W. C. WoLSTENHOLME, AND D. S. VILLARD, United States Rubber Company: Drift and Relaxation of Rubber. 

P. DEBYE, Cornell University: Light Scattering in Solutions. 

R. B. BARNES, American Cyanamid Company: Physical Methods of Analysis of Synthetic and Natural Rubber. 

H. A. RoBINSON, Armstrong Cork Company: Some Physical Properties of Commercial Polymerized Methyl Methacrylate. 

R. H. KELSEY AND J. H. DILLon, Firestone Tire and Rubber Company; H. A. SCHULTZ AND R. C. Bryant, B. F. Goodrich 








Company: A New Extrusion Plastometer. 


I, FANKUCHEN AND H. Mark, Brooklyn Polytechnic Institute: Crystallization of Long-Chain Compounds. 
S. D. GEHMAN AND J. E. FreLD, Goodyear Tire and Rubber Company: X-Ray Observations on the Structure of Rubber and 


the Size and Shape of Rubber Crystallites. 


A. V. ToBoLsky, Princeton University; 1. B. PRETTYMAN AND J. H. DILton, Firestone Tire and Rubber Company: Stress 


Relaxation of Natural and Synthetic Rubber Stocks. 


H. S. Sack, Cornell University: Hysteresis Losses in High Polymers. 
P. H. Emmett, Johns Hopkins University: Measurement of the Surface Area of Carbon Black by Low Temperature 


Absorption. 


R. S. SPENCER AND R. F. Boyer, The Dow Chemical Company: Second-Order Transition in High Polymers. 


ABSTRACTS OF CONTRIBUTED PAPERS 


Al. K-Electron Capture in Radioactive Silver. L. K. 
Hurst* AND M. L. Poot, Ohio State University.—The fail- 
ure to find x-rays! associated with the nuclear K-electron 
capture postulated? as characteristic of Ag’ led to a re- 
investigation of the radioactivities induced in silver. 
Bombardments of Rh, Pd, and Ag with a, p, m, and d have 
given several long periods in the silver fraction. Of these 
the 8.2 day substance is best produced, free from the other 
periods, by the new reaction Rh’ a, n Ag!®, Results have 
thereby been obtained which gave definite evidence for 
K-electron capture. X-rays of 0.59 angstrom have been 
observed which are close to the Ka x-rays of Pd. The 
excited Pd! then emits a gamma-ray of 1.07 Mev. Elec- 
trons of 0.74 Mev of relatively low intensity have been 
observed and presumably are Compton recoils due to the 
gamma-ray. No beta-rays were observed. Electrons of very 
low energy associated with the fluorescence yield were 


present approximately in the amount expected. The 
ionization ratio of x-rays to gamma-rays indicates that 
for each K-electron captured a gamma-ray is emitted. 

* Now in the U. S. Army. 


1L. Alvarez, Phys. Rev. 54, 495 (1938). 
2M. L. Pool, Phys. Rev. 53, 116 (1938). 


A2. Artificial Radioactivity of Columbium and Zir- 
conium. WILLIAM N. Mogurn* anp M. L. Poo., Ohio 
State University Energy measurements have been made 
on the characteristic radiations emitted by radioactive 
columbium and zirconium. The assignments of the various 
periods already have been made through the use of neu- 
tron bombardments.' Deuteron and proton bombardments 
give very strong activities and lead to several new reac- 
tions. The isomer Cb® emits a 0.94-Mev gamma-ray, 
0.11-Mev electrons and 0.78-angstrom x-rays. The latter 
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two observations indicate a very highly internally con- 
verted gamma-ray. The long period in zirconium emits 
gamma-rays of 0.94 Mev and beta-rays of 0.57 and 0.29 
Mev. Proton bombardment of zirconium gave Cb® with 
gamma-rays of 1.0 Mev, beta-rays of 0.59 Mev and x-rays 
characteristic of zirconium. 


* in the U. S. Arm 
fy tt et al., Phys. "Rev. 57, 1179 (1940). 


A3. Beta-Ray Spectra of the Artificially Radioactive 
Isotopes, Pb®* and Bi®’. ARTHUR L. Lutz,* M. L. Poo, 
AnD J. D. Kursatov, Ohio State University —The beta- 
rays emitted by the artificially radioactive isotopes Pb** 
(52 hours) and Bi*’ (6.5 days) were studied by means of a 
newly built beta-ray spectrograph. Exposures of up to 
several half-life lengths were made on x-ray film. The beta- 
rays from the 52-hour (Pb**) activity were found to form a 
line spectrum consisting of the energies, 185, 255, 267, and 
308 kev. The first three were shown to be due to K, L, 
and M electrons, respectively, which were emitted by the 
internal conversion of a 270-kev gamma-ray. The last line 
is believed to be due to Compton recoil electrons expelled 
by a 470-kev gamma-ray. The presence of both gamma- 
rays is shown by lead absorption measurements. The 6.5 
day (Bi*") activity was found to have a more complicated 
line spectrum. Reproductions of the spectrograph films 
will be shown, and the energy values correlated with those 
obtained by absorption measurements of the beta- and 
gamma-rays. Evidence of decay by K-electron capture in 
both cases will also be presented. 


* Now at Wittenberg College, Springfield, Ohio. 


A4. Isolation of Pure Radioactive Cerium. J. D. Kursa- 
Tov AND M. L. Poot, Ohio State University—The long 
life (140 days) radioactive isotope of cerium obtained by 
activation of lanthanum with deuterons has been isolated 
without addition of inactive cerium. The separation of 
cerium in minute concentrations from activated lanthanum 
consists of two processes of coagulation at different pH. 
The first is the coagulation of tetravalent cerium at pH 
3.1-3.3 in the presence of MnO,~, NO;~ and Na* while La 
remains in the solution.! The La contaminating the Ce 
portion after a single separation has been obtained as low 
as 0.005 percent La. It has been established that the 
separation of Ce in minute quantities is independent of its 
concentration, contrary to the concepts of conventional 
chemistry, but is affected by the ions in the solution. The 
Ce is freed from Mn by a second coagulation at a high pH 
after NO;~ is replaced by Cl,~ in solution. The isolation 
of radioactive titanium and ionium,’ light and heavy 
tetravalent elements, has been reported previously. This 
bridges the gap between. 


Quill and H. L. Wilson, unpublished thesis, Ohio State 
Univer, 1937. 
Allen, Pool, Kurbatov, and Quill, Phys. Rev. 60, 425 (1941). 
Hoh Karjavina, and Samojlo, Comptes Rendus Acad. Sci. 
USSR 69 (1930), 
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AS. Investigation of High Energy Gamma-Ray of Na* 
by the Photo-Neutron Method. M. Go_pnaser, G. S. 
KLAIBER, AND G. SCHARFF-GOLDHABER, University of 
Illinois.—For gamma-rays with an energy hy above the 
threshold for photo-disintegration of Be® (1.62 Mev) or D* 
(2.18 Mev) Ay can be determined by the “photo-neutron 
method” previously described.? It consists in determining 
the energy of the photoneutrons from Be or D,0O irradiated 
by the gamma-rays under investigation, with a hydrogen- 
filled ionization chamber. Na™ has a gamma-ray of suffi- 
cient energy to disintegrate both Be and D. For several 
reasons it seemed of interest to apply the photo-neutron 
method for a redetermination of this gamma-ray energy. 
A source of Na™ of several hundred millicuries initial 
intensity, produced by bombarding sodium with 10.5-Mev 
deuterons from the University of Illinois cyclotron, was 
used in this investigation. The photo-neutron energies 
found were (0.34+0.03) Mev for Na*+D, and (1.13+0.05) 
Mev for Na*+Be. This leads to a mean value for the 
gamma-ray energy hy=(2.87+0.05) Mev, in satisfactory 
agreement with two recent measurements, based on an 
analysis of the Compton electron spectrum produced by 
the gamma-rays of Na™, which yielded values of (2.762-0.6) 
Mev? and (2.94+0.06) Mev‘. 

1 G. Scharff-Goldhaber, Phys. Rev. 59, 937 (1941). 


2G. S. Raber and G. Scharff-Goldhaber, Phys. -” 61, 733 (1942). 
*L. G. Elliott, M. Deutsch, and A. Roberts, Phys. Rev. 63, 386 


ae Cx. Mandeville, Phys. Rev. 63, 387 (1943). 
A6. Energy of Gamma-Rays from Na**-Mg"™. W. E. 
OGLE AND P. GERALD KRUGER, University of Illinois.— 
Gamma-rays from a source of 250+ 200 mC of radiosodium 
were used to produce-+telectron pairs in the gas (air) of a 
cloud chamber. A lead collimator 40 inches long, 4 inch 
high, and 1 inch wide provided adequate collimation so 
that it was easy to apply rigid rules to the selection of pairs. 
Only those showing no evidence of scattering and occurring 
within a 20° cone in the forward direction were accepted 
as good pairs. Six thousand pictures yielded 31 pairs, 
including two “triplets.” Twenty-seven of these pairs 
occur between 2.41 and 2.89 Mev and give an average 
gamma-ray energy of 2.69 Mev. The other 4 pairs indicate 
two higher energy gamma-rays at 3.22 and 3.61 Mev. 
There is some evidence that the 2.69-Mev group may be 
complex. This point is being investigated further. All of 
these energy values are subject to correction when the 
calibration of the mutual inductance used to measure the 
magnetic field is checked by the National Bureau of 
Standards. The gamma-ray energies calculated from the 
two “‘triplets’’ (pairs formed in the field of an electron) are 
2.67 and 2.66 Mev. The accuracy and resolving power of 
this method of measuring gamma-ray energies will be 
discussed. The results of other investigations! indicate the 
existence of gamma-ray energies of 1.38+.03 and 2.85+.09 
Mev. These values are compatible with the present results, 
since the pair production at 1.38 Mev is very small and the 
gamma-rays above 3 Mev are very weak. 

1 Elliott, Deutsch, and Roberts, rove, | Rev. 63, gos 985); C. E. 
Mandeville, Phys. Rev. 63, 387 (1943); M. Goldhaber Klaiber 


and G. Scharff-Goldhaber. Abstract No AS, to be 
meeting. 
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A7. A Cyclotron Which Allows the Accelerated Particles 

to Emerge in a Direction Parallel to the Dee Interface. 
P. GERALD KRUGER, J. R. RicHARDsON, G. K. GROET- 
ZINGER, E. M. Lyman, W. E. Ocie, M. E. Ngtson, G. 
SCHWARZ, J. B. GREENE, R. T. ANDERSON, N. C. CoLBEy, 
R. W. Lee, C. E. McCLetian, D. Scac, Ltoyp Smits, 
F. K. TALLMADGE, University of Illinois.—The construc- 
tion of the dees, deflection channel and vacuum chamber 
of a cyclotron will be described. The design allows the 
accelerated ions to pass through an electrostatic deflect- 
ing field for one-quarter of a revolution and then to pro- 
ceed through the magnetic fringing field for another quarter 
of a revolution before they leave the cyclotron vacuum 
chamber. This provides that the ions emerge in a direction 
which is parallel to the dee interface and perpendicular 
to the principal azimuthal plane through the long axis 
of the cyclotron magnet. There are several advantages to 
this arrangement: (1) Relatively low deflecting voltages 
are necessary. (2) The emerging ions at 10 Mev (100 we or 
more) are satisfactorily mono-energetic. (3) Because of the 
favorable direction of emergence of the ions, it is easy to 
“pipe” the beam to a point many feet from the cyclotron. 
This allows adequate shielding to be placed between the 
cyclotron and other experimental apparatus, such as cloud 
chambers or counting apparatus. (4) It is convenient to 
place magnetic analyzing apparatus in the beam path near 
the cyclotron (but outside of the cyclotron magnetic field) 
and thus attain small currents of very monoenergetic ions 
for charged particle scattering and similar experiments. 


A8. The Refractive Index in Electron Optics. I. Opa- 
TOWSKI, Armour Research Foundation, Illinois Institute of 
Technology.—The application of methods of geometrical 
optics for the investigation of electron trajectories requires 
the knowledge of the refractive index yu of the equivalent 
optical problem. A general formula for yu in terms of the 
electrostatic potential V and of the magnetic vector poten- 
tial A is known. If A+0 this formula contains a unit vec- 
tor s parallel to the electron velocity, which is not known in 
general. It is, therefore, important that s be eliminated. 
This can be done only if a suitable integral of the equations 
of motion of the electron is known. If V and A are sym- 
metric about an axis, 8 may be eliminated from yu by means 
of a momentum integral (Glaser). This elimination is 
generalized in the present paper to include cylindrical 
fields (not necessarily symmetrical) as well as some more 
general fields with a transversal field component. The re- 
sult is reached by interpreting the magnetic scalar poten- 
tial as a velocity potential of an ideal fluid and making use 
of the corresponding stream function to obtain an expres- 
sion for w. Both relativistic and constant mass are con- 
sidered. 


A9. Influence of Pressure on Intermetallic Diffusion. 
F, J. RADAVICH AND R. SMOLUCHOWSKI, General Electric 
Company.—The influence of elastic deformation of crystal 
lattice on diffusion is of considerable theoretical interest 
because it might help to estimate the various factors in- 
volved in the mechanism of intermetallic diffusion. On 
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the other hand, it has also a bearing upon many technical 
problems in which diffusion occurs in strained material. 
The experimental difficulty lies in obtaining high pressures 
at temperatures which are high enough to produce well- 
measurable concentration gradients. The experiments were 
performed on the Al—Cu system. Two bars, one of pure 
aluminum and the other of a 9 atomic percent (4 weight 
percent) Al—Cu alloy were welded together in vacuum. 
The welded sample was placed in a mold and subjected to 
high pressures at 500°C for a week. The analysis was made 
spectrochemically and the diffusion coefficients calculated 
from the concentration vs. distance plots. There was no 
measurable change of the diffusion rate observed at the 
pressure of 1500 kg/cm*. However, at 7000 kg/cm*, a 
decrease of the rate of diffusion was found corresponding to 
a decrease of the diffusion coefficient by about 30 percent 
at a copper concentration of 2 atomic percent. 


Al0. A Series Generator as a Negative Resistance. G. 
H. Fett, University of Illinois.—A series generator in the 
operating range below maximum voltage possesses both 
d.c. negative resistance (—V/J), and incremental nega- 
tive resistance (—dV/dJ). It is shown that the incre- 
mental negative resistance is equal to Kn—R, where n is 
the speed, K the induction factor per unit of field current, 
and R is the armature circuit resistance. Since the series 
generator has the incremental resistance property, it will 
support oscillations if the tuned circuit has the proper 
constants. !t is shown that a generator load consisting of a 
separately excited motor can be represented by a condenser 
of sufficient size (in the order of farads) that oscillations 
at mechanically feasible frequencies of 0.2 to 1.5 cycles per 
second can be obtained. It is shown analytically that the 
magnitude of the capacitance varies directly with the 
moment of inertia of the motor armature and inversely 
with the square root of the induction factor. It is shown 
experimentally that the frequency of oscillation of the 
circuit depends upon the motor induction factor and mo- 
ment of inertia, and slightly upon the circuit resistance. 
The amplitude of oscillation is controlled by the circuit 
resistance. 


All. New Experiments About the Magnetic Current. 
FELIX EHRENHAFT, New York City.—Experimentally we 
have established the existence of magnetic currents in our 
discovery of magnetic ions, magnetolysis of water, and the 
electric field around magnetic currents.'! The following is a 
continuation of these experiments. In the darkfield we 
observe that microscopic testbodies (e.g. chromium, 
manganese, nickel) move in spiral paths or parts of spirals 
around a completely constant vertical homogeneous mag- 
netic field in air of atmospheric pressure, created in a 
magnetic condenser (magnet poles of soft Swedish iron 
apart 0.5-1.5 mm, diameter 6-12 mm). In dilute copper 
sulfate solution blue light scattering particles of colloidal 
size become visible in the darkfield between the aforesaid 
magnetodes. Each single colloid particle spirals up or 
down around the constant homogeneous vertical magnetic 
field ending its movement on the upper or lower magnetode. 
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These are the well known copper-particles that precipitate 
on iron.? The described experiments in gases and liquids 
prove clearly the existence of bodies charged simul- 
taneously electrically and magnetically. Because of their 
electric charges they move around the magnetic field 
along the electric force lines closed in circles. Because of 
their magnetic charges, north or south, they move up or 
downwards. The helix is the resultant path.’ 


1 Felix Ehrenhaft, Phys. Rev. 60, 169 (1941); 61, 733 (1942); 63, 216, 


1943); 64, 43 (1943). ’ : 
461 comm see H. V. Juepter, La lumiere electrique 10, 468 (1883). 
3 The experiments can be seen at C. Zeiss, Inc., New York City. 


Bl. Counter Corrections at High Counting Rates. 
TRUMAN P. KouMAN, University of Wisconsin.—A method 
has been developed for determining counter losses at high 
rates which does not involve assumptions regarding counter 
mechanism. Several pairs of sources are used, the members 
of each pair being approximately equal. For each pair the 
sources are counted separately and together, and a first 
approximation to the fractional correction for the combined 
source count is calculated on the assumption of no losses 
in the separate source counts. The corrections are plotted 
against recorded rates and the best smooth curve drawn. 
Second approximations to the combined source corrections 
are made by using this curve to correct the separate source 
counts, and a better curve obtained. The process of 
successive approximations is repeated until the curve 
converges to a stationary position. With a Geiger-Miiller 
coanter and Neher-Harper circuit, the apparent resolving 
time decreased markedly with increasing rate, as expected 
from theory. Simplified assumptions about the mechanism 
yielded a theoretical expression for the average deadtime 
tw as a function of true rate N: ry =tmin+(tmax—Tmin)p/ 
(o+N). tmax is the deadtime following an isolated pulse, 
Tmin is the deadtime of the smallest recordable pulse, and 
p characterizes the rate of recovery of the field strength 
around the counter wire. This equation fits the observations 
satisfactorily, 


B2. Recording Negative Ions and Electrons in a Gas 
with a Geiger-Miiller Counter. R. T. K. Murray, Poly- 
technic Institute of Brooklyn.—Free electrons and negative 
ions present in a gas have been recorded by allowing them 
to enter a counter through a window consisting of a small 
open aperture. The pressure outside the counter could be 
adjusted above or below the optimum counter pressure by 
suitable pumping. Negative ions once in the counter suffer 
electron detachment at the high values of X/p and are 
therefore recorded. Electrons and negative ions can be 
introduced into the counter by diffusion, by convection in 
the gas flow and by electrical convection under an applied 
electric filed. For free electrons the electric drift velocity 
for weak fields is still large compared to the jet velocity 
which has an upper limit slightly greater than the velocity 
of sound thereby making it possible to convect free elec- 
trons into the counter against an outwardly directed gas 
jet at the aperture. Electrons from a thermionic source in 
electronegative oxygen at 4-mm pressure have been re- 
corded. This is possible since the average distance for elec- 
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tron attachment is large compared to the distance between 
the source and the aperture. It is planned to develop the 
apparatus further for the efficient detection of space- 
localized low energy ionizing particles. 


B3. Study of Cosmic-Ray Air Showers Using the Method 
of Coincident Bursts in Two Unshielded Ionization Cham- 
bers. L. G. Lewis anp E. W. Lewis, University of Chicago. 
—Simultaneous bursts of ionization were observed in two 
unshielded thin walled ionization chambers. These cham- 
bers, which were 35 cm in diameter, were filled with argon 
to a pressure of 100 cm Hg. The signals, which the bursts 
produced on the grids of two FP-54 tubes, were amplified 
by d.c. amplifiers employing inverse feedback. The out- 
put of each amplifier was connected to a sensitive string 
galvanometer. Both galvanometer deflections were re- 
corded simultaneously and continuously on the same 
photographic film. This method of observation makes it 
possible to plot size-frequency distribution curves for 
single bursts in each chamber, as well as a family of curves 
giving the rate of burst coincidences as a function of separa- 
tion of the chambers. This last family of curves, with size 
of burst as a parameter, yields information concerning the 
lateral density distribution of extensive cosmic ray showers. 
Measurements have been made near Echo Lake, Colorado 
(elevation 10,650 ft.) and further investigations are con- 
templated at Idaho Springs, Colorado (elevation 7556 ft.), 
and at Chicago. Analysis of these coincidence curves for 
different altitudes should yield the energy spectrum for 
the shower producing radiation. 


B4. The Origin of Large Cosmic-Ray Bursts. R. E. 
Lapp, University of Chicago——Large cosmic-ray bursts 
under different thicknesses of iron have been investigated 
with the simultaneous use of an ionization chamber and 
Geiger- Mueller counters. From an analysis of the transition 
curve obtained for iron, three types of bursts are dis- 
tinguished: (1) é-bursts observed in an unshielded cham- 
ber; (2) tm-bursts produced in a 12-cm iron shield; (3) 
t,-bursts originating under a shield equivalent to at least 
20 radiation units. The ¢o bursts were found to be coincident - 
with the simultaneous triggering of the G-M counter ar- 
rangements which were designed to detect extensive atmos- 
pheric showers (A showers). Thus a ¢) burst represents a 
high density of particles in an A shower (the so-called core). 
On the other hand, only about 5 percent of the ¢ bursts 
were found to be coincident with A showers; this is re- 
garded as further evidence that the great majority of t, 
bursts at sea level are produced by mesotrons. Analysis of 
the burst data obtained with a 12-cm iron shield shows that 
the ¢,, bursts are a complex of t and ¢, bursts and in addi- 
tion it seems that some of the ¢,, bursts originate neither 
from the cores of A showers nor from mesotrons. Experi- 
ments are in progress to determine the origin of this fraction 
of the ¢,, bursts. 


BS. An Analysis of Large Cosmic-Ray Bursts. MARcEL 
SCHEIN AND R. E. Lapp, University of Chicago——Through 
the cooperation of the Carnegie Institution over five years 
of burst data obtained at Cheltenham, Teoloyucan, and 
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Huancayo were made available to us. These data have been 
analyzed and size-frequency burst distribution curves were 
obtained for 72, 2285, and 3350 meters elevation. Because 
of the great amount of data available these curves are more 
accurate than those! published earlier; this is especially 
true for the largest size bursts. The sea-level data are dis- 
cussed with reference to the theoretical calculations of 
Christy and Kusaka,* and a discussion is given of the possi- 
bility of determining the spin of the mesotron' from the 
sea-level size-frequency distribution which is presented. 
The origin of bursts under thick shields at high altitude is 
also considered. 


1M. Schein and P. S. Gill, Rev. Mod. Phys. 11, 267 (1939). 
?R. F. Christy and S. Kusaka, Phys. Rev. 59, 414 (1941). 


B6. Production of Mesotrons by Non-Ionizing Radia- 
tion at Altitudes of 10,618 ft. and 14,230 ft. JuLrus TaBIN, 
University of Chicago.——The production of mesotrons in 
various materials was studied by means of an apparatus 
consisting of 46 Geiger-Mueller counters arranged in 12 
three-, four- and fivefold coincidence sets. Each individual 
pulse actuated a neon light which was recorded photo- 
graphically. The producing layers consisted of different 
thicknesses of paraffin, aluminum, iron, and lead. Special 
care was taken to eliminate the influence of air showers. It 
was found that at the above-given altitudes the rate of 
production of both single and multiple mesotrons by non- 
ionizing rays is a very small effect and did not exceed under 
optimum conditions more than 2.5 percent of the intensity 
of the penetrating component. The maximum production 
of single mesotrons was found to exist in paraffin. The pene- 
trating power of the produced mesotrons was determined 
in lead and it was found that a large fraction of these 
mesotrons was stepped by a thickness of 11 cm. Experi- 
ments are in progress to investigate the nature of the 
mesotron producing radiation. 


B7. Secondary Radiation Produced by Cosmic-Ray 
Primaries in Paraffin. Mario IoNA JR., MARCEL SCHEIN, 
AND JuLius TaBIN, University of Chicago.—In recent 
stratosphere balloon experiments counter coincidences have 
been observed between the passage of vertical penetrating 
particles and multiple secondaries emerging below layers 
of paraffin. From electromagnetic theory no appreciable 
shower production is to be expected in the thicknesses of 
paraffin used. However, the experiments show that a large 
proportion of the single penetrating rays passing through 5 
or 10 cm of paraffin in the stratosphere emerge as several 
particles. The relative number of these secondaries in- 
creases rapidly with altitude in the upper atmosphere. This 
strongly indicates that this effect is caused by a direct 
interaction of the primary cosmic radiation with the nuclei 
of paraffin. The majority of the multiple secondaries are 
ejected under a very small angle below the paraffin. 


B8. Electrons in Equilibrium with Mesotrons. H. E. 
STANTON, University of Chicago (Introduced by Marcel 
Schein).—The problem of the origin of the soft component 
of cosmic rays in the atmosphere which arises from the 
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mesotron component is considered in such a way as to 
enable a calculation of the energy spectrum of the electrons, 
The theory is based on the cross sections derived by Bethe 
and Heitler and the simple theory of mesotron decay. It jg 
assumed that the mesotron has spin zero and decays into 
a single neutrino and electron with emission equally prob. 
able in all directions in the reference frame of the mesotron, 
A modified set of diffusion equations is set up and forma] 
solutions found with the boundary conditions for no in. 
coming electrons or photons. By the use of approximate 
methods, it is possible to calculate the spectrum of the 
electron component once the mesotron spectrum is known, 
and these results can be compared with experimental values 
obtained at the lower altitudes. 




















B9. Non-Shower Producing Particles in A Showers. P. 
AUGER AND AN. RocoziNnski, University of Chicago.— 
Penetrating particles in extensive atmospheric showers (A 
showers) were investigated by means of Geiger-Mueller 
counters. It was found that single particles emerging below 
lead screens up to 12.5 cm thicknesses, were associated 
with particles which traversed counters placed in horizontal 
distances up to 15 meters from those counters which were 
located below the lead screen. A careful study of the ob. 
served phenomena shows that part of these events could 
be attributed to knock-on processes arising from collisions 
between single mesotrons and electrons. 
















C1. On Some Physical Properties in the Problem of 
Rubber. GeorGe ANToNoFF, Fordham University.—In the 
course of his work on surface tension, the author evolved 
some methods which can be applied to the study of the 
problem of rubber. Two will be mentioned here. There 
exists parallelism between rubber mixed with pigments and 
paints which are a mixture of oil and pigments. Pigments 
which increase the breaking stress of films of oil will also 
reinforce the rubber and vice versa. In pastes, as a function 
of concentration, an abrupt change in properties is observed 
when fluid paste becomes plastic. Similarly, in rubber, 
there is an abnormal zone where the properties change as 
soon as a certain concentration of pigment is attained, 
Properties of solutions of rubber vary with time. Similar 
changes can be observed in liquids leading to formation of 
rubber. Some fraction from turpentine has been heated. 
Its properties changed with time. Suddenly the transparent 
liquid changed into a brown rubber-like mass. Analogous 
but small changes can be observed in ordinary liquids. They 
can be made apparent by surface tension method described 
by the author.! 


1G. Antonoff, J. Phys. Chem. 46, 492 (1942). 




















C2. A Proposed Substitute for Instrument Oil. S. BLoom- 
ENTHAL, Automatic Electric Laboratories, Chicago, Illinois. 
—Porpoise jaw oil is superior to any known petroleum 
derivative for lubrication of mechanical instruments that 
must operate for long periods of time without attention, 
Because of the present scarcity of the first mentioned oil, 
a substitute lubricant was looked for among certain esters 
and oils that are fluid at —30°F. Help in this search was 
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obtained from the Laboratory of the Distillation Products 
Company, Rochester, New York. Tenacious low vapor 
pressure films on metal may explain the exceptional 
lasting qualities of porpoise jaw oil. CeH«(COOC,Hy)s, 
CsH1e(COOCsH)2 and Distillation Products Company 
Synthetic Lubricating Oil have low vapor pressures and are 
fluid at —30°F. Far better results were obtained with these 
substances in the lubrication of automatic telephone dials 
than with light petroleum or fish liver oils, and it is believed 
that a satisfactory substitute for porpoise jaw oil has been 


found. 


C3. The Energy of Dissociation of Hg:. J. G. WINANS 
AND M. P. HeE11z,* University of Wisconsin.—D for Hg: has 
been determined by a modification of the method used by 
Koernicke.! The absorption intensity of the 2345 Hg band 
was measured photographically at two different tempera- 
tures with constant density. With the equation of Gibson 
and Heitler, D was calculated from measurements of two 
temperatures and the ratio of absorption intensities at 
these temperatures. Constant vapor density was obtained 
by using an absorption cell containing a drop of mercury 
small enough to be completely evaporated at the low- 
est temperature measured. The final value of D was 
0.060+0.003 volt. This agrees within experimental error 
with the value obtained by Koernicke. 


* Now at Cedarburg, Wisconsin. 
1 E. Koernicke, Zeits. f. Physik 33, 219 (1925). 


C4. Raman Spectra of 2,5-Dimethyl-1,5-hexadiene and 
2,5-Dimethyl-2,4-hexadiene. Forrest F,. CLEVELAND, [lli- 
nois Institute of Technology.—Continuing work on olefinic 
hydrocarbons, Raman frequencies, relative intensities, and 
depolarization factors have been obtained for (A) 2,5- 
dimethyl-1,5-hexadiene and (B) 2,5-dimethyl-2,4-hexa- 
diene. The relative intensities and depolarization factors 
were obtained by use of a Gaertner microdensitometer. The 
two molecules are isomers differing only in the location of 
the two double bonds contained in each molecule. In A the 
double bonds are at the ends of the chain and the molecule 
thus contains two XYC=CHgz groups. Mean values ob- 
served in the previous work on non-cyclic mono-olefins for 
the five frequencies characteristic of this group are 1299, 
1412, 1649, 2977, and 3075 cm. Corresponding frequencies 
for A are 1303, 1414, 1651, 2977, and 3075 cm™; the diolefin 
A thus has the same frequencies as the mono-olefins. In B 
the double bonds are symmetrically located in the middle 
of the molecule; the molecule thus contains the group 
XYC=CHZ. Dipentene, the only olefin investigated previ- 
ously in this laboratory that contained this group, had the 
characteristic frequencies 1246, 1682, 3012, and 3050 cm™. 
Corresponding frequencies for B are 1227, 1663, 3012, and 
3040 cm=. The lowering of the 1682 dipentene frequency 
to 1663 in B is due to the conjugation of the double bonds. 


CS. Polymer Spectra of a Cyanine Dye. RicHarp W. 
Mattoon, University of Chicago (Introduced by J. Franck). 
—In dilute aqueous solutions 1,1’-diethyl-2,2’-cyanine 
chloride has two broad absorption bands with maxima at 
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4900 and 5230A, while at concentrations >10~* molar a 
new, intense, and exceedingly narrow absorption band ap- 
pears at 5730A (G. Scheibe, E. Jelley, and others). Accord- 
ing to Scheibe this new band belongs to a system of long 
chains of molecules formed by polymerization. In the 
present investigation the occurrence and behavior of this 
“P band” were studied for thin dye films deposited on glass 
from solutions containing water. The P band shifts about 
50A to >A, becomes broader and weaker, and finally dis- 
appears as the water is removed by pumping and/or by 
increasing the temperature. This behavior is reversible if 
water vapor is readmitted. The position of the P band 
enables a measurement of the humidity of the atmosphere 
surrounding the film. A second, weaker absorption band at 
5430A behaves in a manner similar to that of the P band in 
its dependence upon water vapor pressure and temperature. 
Measurements of fluorescence polarization, etc., taken from 
the literature and the present results are used for an inter- 
pretation of these bands. 


C6. Isotope Shift in Magnesium. L. G. MUNDIE AND 
K. W. MEISSNER, Purdue University —The previous in- 
vestigations concerning the isotope shift of magnesium 
have been repeated and extended by means of an improved 
atomic beam source and a Perot-Fabry interferometer. 
The values of Meissner have been checked to within experi- 
mental error in most cases. In addition, several new lines 
have been resolved and their isotope shift determined. The 
31P—m'"D series has been resolved from m=3 out as 
far as the member m=11. In this series the values of 
Av(=vse— a4) seem to converge to a value of about 0.068 
cm~ rather than the value 0.056 cm obtained by Meiss- 
ner. The isotopic shifts of the members m=5 and m=6 of 
of the series 3'P—m'S have been measured. The results 
obtained are discussed in connection with the theoretical 
work of Vinti, who assumed the splitting to be solely a mass 
effect. The experimental and theoretical values of Av agree 
in sign, but the former are usually larger by a factor of 
about two. 


C7. Forbidden Lines of Bismuth Bi I. S. Mrozowsk1, 
University of Chicago.—A group of forbidden lines corre- 
sponding to transitions between energy levels belonging to 
the electronic configuration 6p* was observed in the spec- 
trum of neutral bismuth. The experimental arrangement 
for excitation of the spectrum was similar to the one used 
by the author in a work on forbidden lines of lead (con- 
figuration 6p*) and also in a recent work on the line spec- 
trum of bismuth.' The photographs of 6 lines and their 
hyperfine structures were obtained with the 30-foot grating 
spectrograph. The lines observed are: 43014, 44597, 44615, 
45640, 46476, and 48755. The two last ones are especially 
interesting, since the structure of the first one gives a con- 
firmation of the value of the hyperfine structure separation 
reported by the author! for the ground level *5%;_ (a value 
almost 3 times as large as the one accepted in the past). 
The last line is situated quite far in the infra-red and its 
photograph can be obtained only with considerable diffi- 
culty. The hyperfine structure of all the lines is showing the 
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predominantly magnetic dipole character of the forbidden 
radiation with a possible only slight admixture of electric 
quadrupole radiation, since no components corresponding 
to a change AF= +2 were observed. 

1S. Mrozowski, Phys. Rev. 62, 526 (1942). 


C8. A New Adsorption Isotherm and Polymolecular 
Films on Solids. WILLIAM D. HARKINS AND GEORGE JURA, 
University of Chicago and the Universal Oil Products Com- 
pany.—Films adsorbed from vapors on the surfaces of 
solids exhibit phases of different types, among which gase- 
ous and condensed phases are particularly important. The 
equation for a condensed monolayer, either liquid or solid, 
on an aqueous subphases may be written 

x=b—aad (1) 
in which =z is the film pressure, and o the molecular area. 
By a method devised by Bangham it is possible to calculate 
the film pressure (x) for a film on a solid, and it is found 
that at the lower film pressures at which Eq. (1) is valid, 
the film is monomolecular, but that, as the film pressure 
increases the film becomes polymolecular without any 
change in the value of the slope (a). By the use of the 
Gibbs adsorption Eq. (1) was changed into 

log (p/po) = B— (A/V?) (2) 
in which po is the gas pressure at saturation, A and B are 
constants for the condensed film, V is the volume of gas 
adsorbed (STP) when the pressure is p. Equation (2) is 
found to give a linear relation over a very much greater 
range of pressure than any earlier adsorption isotherm. The 
area of a very finely divided crystalline solid may be de- 
termined by an absolute method described earlier. 


C9. An Adsorption Method for the Determination of the 
Area of a Solid. GEORGE JURA AND WILLIAM D. HaARKINs, 
University of Chicago and the Universal Oil Products Com- 
pany.—An adsorption isotherm developed by the writers 
has the form 


log (p/po) = B—(A/V*) 
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in which B and A are constants and V is the volume (STP) 
of gas adsorbed at the pressure p, and fp is the saturation 
pressure. It is found that the areas of the surfaces of the 
solids (per gram) are proportional to the square roots of 
the slopes A, or 

==kAt 


where = is the area per gram. An absolute method for the 
determination of the area of a finely divided crystalline 
solid consists in dividing the heat of immersion of the solid, 
covered by a duplex film of the liquid, in the liquid itself, 
The amount of heat divided by the surface energy of the 
liquid gives the area of the solid, and thus the value of the 
constant k is obtained. Thus (1) becomes 


==4.06 Ai 


The equation of Brunauer, Emmett, and Teller gives the 
number of molecules in the first monolayer. If the area of 
the nitrogen molecule is assumed to be 14.8A?, the areas 
obtained by our method agree with theirs within plus or 
minus 9 percent. 


C10. The Stability of a Mercury Arc as a Function of 
Cathode Pool Area. Paut L. COPELAND AND CARROLL 
Lurcy, Jilinois Institute of Technology.—The work of Cope- 
land and Sparing' on the increase in the stability of a 
mercury arc with the cathode spot anchored suggested a 
systematic statistical study of arc life as a function of the 
size of the cathode pool. An experimental low pressure 
mercury arc was started repeatedly and the measured times 
of extinction were recorded. The results established a corre- 
lation between the surface area of the cathode pool and 
the average life of the arc. When the cathode pool was large 
the life was relatively short, and the average life increased 
rapidly as the area of the cathode pool was decreased. 


1 Copeland and Sparing, Phys. Rev, 59, 115 (1941). 





